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FOREWORD 


Oils  report  summarizes  research  carried  out  under 
contract  AF  49(638)- 544  at  the  Graduate  School  of  Aerospace 
Engineering,  Cornell  University,  Ithaca,  Rev  York*  This  contract 
has  been  monitored  by  the  Mechanics  Division  of  the  Office  of 
Scientific  Research,  USAF. 

The  vork  reported  here  was  begun  before  the  extended 
visit  of  Professor  D.  A.  Spence  in  1959-80,  was  pursued  under 
his  direction  during  that  academic  year  and  more  recently  has  been 
continued  in  his  absence.  As  the  reader  will  discover,  the  conclusions 
reached  are  in  some  degree  critical  of  the  solutions  provided  by 
Prof.  Spence  for  the  integro- differential  equations  set  up  by  the 
author,  John  C.  Erickson,  Jr.  Since  Prof.  Spence's  return  to 
England,  where  he  is  employed  at  the  Royal  Aircraft  Establishment, 
he  has  continued  his  Interest  in  the  investigation  and  sincere 
thanks  must  go  to  him  for  his  many  comments. 

This  is  the  last  Technical  Note  to  be  prepared  under 
the  above-mentioned  contract,  which  has  now  been  replaced  by  a 


grant. 


ABSTRACT 


A  linearized  model  for  the  Incompressible,  inviscid, 
ir rotational,  and  unsteady  flow  about  a  thin  airfoil  with  jet- flap 
is  formulated. .  The  unsteady  problems  considered  are  the  transient 
and  oscillatory  deflection  of  the  Jet,  plunging  and  pitching  of 
the  airfoil,  deflection  of  a  "blown- flap, "  and  also  the  penetration 
of  a  sharp-edged  gust.  Justification  is  given  for  representation 
of  the  jet,  in  the  limit  of  high  speed,  small  thickness,  and  constant 
momentum- flux  strength,  by  a  vortex  sheet,  across  which  there  is  a 
pressure  difference  proportional  to  the  momentum- flux  strength  and 
inversely  proportional  to  the  local  radius  of  curvature  of  the  jet. 
The  dynamic  and  kinematic  interaction  of  the  main  stream  with  the 
vortex  sheets  representing  the  airfoil  and  jet  are  shown  to  be 
described  by  a  coupled  set  of  equations  consisting  of  a  third-order 
partial  differential  equation  and  a  singular  Integral  equation, 
along  with  appropriate  boundary  conditions.  The  properties  of  these 
equations  and  their  relationship  to  classical  unsteady  thin- airfoil 
theory  curd  steady  jet- flap  theory  are  discussed. 

For  small  momentum- flux  strength,  and  for  either  small-time 
after  initiation  of  transient  mcftlon  or  high  frequency  of 
oscillation,  a  transformation  is  made  which  leads  to  a  simplified 
form  of  the  governing  equations.  All  types  of  airfoil  motion  are 
reduced  in  this  approximation  to  a  single  set  of  equations,  whose 
solution  is  found  in  terms  of  the  solution  for  jet  deflection. 
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Spence  (196IB)  attempted,  by  making  a  farther  approximation, 
a  solution  for  mall -time  after  instantaneous  Jet  deflection.  Ibis 
solution  is  found  to  be  incorrect,  and  it  is  shown  that  no  solution 
can  be  found  in  the  sense  of  this  further  approximation.  Spence’s 
related  solution  for  high^ frequency  oscillatory  Jet  deflection  is 
also  found  to  be  Incorrect,  and  a  tentative,  corrected  solution  is 
proposed  and  discussed. 
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INTRODUCTION 


The  jet- flap  principle  has  heen  extensively  studied,  both 
expe -H  mental 1 y  and  theoretically,  in  recent  years  since  the  pioneer¬ 
ing  work  of  Davidson  (1956),  Stratford  (1956),  Malavard,  Polsson- 
Quinton  and  Jousserandot  (1956),  and  Helmbold  (1955)*  The  principle 
is  briefly  this:  a  thin,  high-speed  jet  of  air  is  ejected  at  or 
near  the  trailing  edge  of  an  airfoil.  Besides  the  direct- reaction 
lift,  which  acts  upon  the  internal  jet  ducting,  additional  lift  is 
obtained  due  to  the  effect  of  the  curved  Jet  on  the  airfoil  external 
pressure  distribution.  Furthermore,  this  modified  external  pressure 
distribution  accounts  for  recovery  of  very  nearly  the  total  thrust 
of  the  Jet,  independent  of  jet  deflection  angle,  the  so-called 
"thrust  hypothesis,"  cf.  Yen  (i960).  In  the  Jet- flap,  then, 
has  been  found  a  premising  means  of  integrating  the  lift  and 
propulsion  of  an  airfoil.  A  thorough  review  of  the  jet- flap 
literature  has  been  given  recently  by  Korbacher  and  Srldhar  (i960). 

This  integration  has  generated  interest  for  application 
not  only  to  airplane  wings  for  STQL  performance,  but  also  to  airfoil 
applications,  e.g.,  helicopter- rotor  blades,  cf.  Richards  and  Jones 
(1956)  and  Dorand  ( 1 959) >  Jet-engine  compressor  blading,  cf.  Clark 
and  Grdway  (1 959) >  Brocher  ( 1 961 )  and  Paul on  (1959) >  and  very 
recently  to  hydrofoils,  cf.,  Ho  ( 1 961 ) . 

In  view  of  the  proposed  applications  of  the  Jet- flap, 
it  would  seem  desirable  to  extend  the  analysis  to  unsteady  problems. 
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For  if  the  jet- flap  is  used  as  a  controlling  device,  e.g.,  for 
cyclic  control  of  helicopter- rotor  blades,  or  for  aircraft  in  slow 
flight,  the  lift  response  to  time- dependent  Jet  deflections  is  of 
great  importance.  In  flutter  stability  calculations,  the  lift  and 
pitching  moments  of  an  oscillating  Jet- flapped  airfoil  are  required. 

Gbe  motivation  for  the  study  of  unsteady  jet- flapped  airfoil  theory 
therefore  arises  from  the  same  considerations  which  motivated  the 
classical  theory  of  unsteady  airfoil  motion. 

To  study  the  steady- state  lifting  properties  of  jet- flapped 
airfoils,  a  model  has  been  formulated  independently  by  Malavard  ( 1 957 ) , 
Helmbold  (1955)  and  Spence  (1956)  -  which  will  be  referred  to 
subsequently  as  I.  In  this  model,  the  non- homogeneous  flow  of  the 
jet  embedded  in  the  main- stream  is  treated  by  representing  the  jet 
by  a  vortex  sheet,  across  which  there  is  a  pressure  difference 
proportional  to  the  jet  momentum  flux  and  inversely  proportional 
to  the  jet  radius  of  curvature.  This  model  has  been  shown  by  Spence 
in  I  to  be  a  good  approximation  when  the  Jet  velocity  is  very  much 
greater  than  the  free- stream  velocity,  and  when  the  Jet  is  sufficiently 
thin.  With  this  as  the  model  for  the  jet,  the  problem  can  be  linearized 
in  typical  thin- airfoil  fashion.  The  linearized,  two-dimensional 
problem  has  been  solved  in  a  rheoelectrlc  analogy  by  Malavard  ( 1 957) > 
numerically  by  Spence  in  I  and  Spence  (1958);  and  finally 
analytically  by  Spence  (196IA)  -  to  be  referred  to  as  II.  The 
corresponding  problem  for  supercavitating  hydrofoils  has  been 
solved  numerically,  following  I,  by  Ho  (1 961 ) •  Three-dimensional 
theories  have  also  been  put  forth,  but  are  beyond  the  scope  of  this 
present  research. 
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In  the  present  research,  the  unsteady  two* dimensional 
lifting  problem  of  a  Jet- flapped  thin  airfoil  in  incompressible 
flow  is  formulated  as  an  extension  of  the  Malavard-Helmbold- Spence 
model.  The  problems  considered  are  those  where  (i)  the  airfoil 
is  performing  some  time- dependent  motion  normal  to  a  mean,  steady 
position,  or  ii)  the  Jet- deflect! on  angle  is  time- dependent. 

The  Jet  momentum  flux  at  its  exit  from  the  airfoil  is  assumed  , 
constant,  independent  of  time.  Furthermore,  the  Jet  is  assumed 
to  be  fully  developed  in  length  prior  to  the  onset  of  the  unsteady 
motion,  i.e.,  the  time- dependent  motion  is  superimposed  upon  a 
fully  developed  steady- state  configuration. 

In  Chapter  1  the  non- homogeneous  flow  problem  of  a  Jet 
and  mainstream  of  different  total  pressures  is  treated  to  find 
a  model  for  the  dynamical  interaction  of  the  Jet  and  mainstream. 

The  limit  of  a  very  thin,  very  high  speed,  constant-momentum- flux  Jet 
is  taken,  and  the  consequent  representation  of  the  Jet  by  a  vortex 
sheet  is  Justified.  This  vortex  sheet  is  characterized,  exactly 
as  in  the  steady  problem,  by  its  support  of  a  pressure  difference 
proportional  to  the  Jet  momentum  flux  and  inversely  proportional 
to  the  instantaneous  local  radius  of  curvature  of  the  Jet. 

The  resultant  flow  problem  is  linearized  by  the 
assumptions  of  classical  thin- airfoil  theory  in  Chapter  2,  the 
airfoil  also  being  represented  by  a  vortex  sheet  with  the  appropriate 
dcwnwash  boundary  condition.  Consideration  of  the  pressure 
difference  across  the  vortex  sheet  representing  the  Jet  leads 
to  a  third- order  partial  differential  equation  relating  the  Jet 
vortex  strength  and  the  Jet  ordinate,  or  alternately,  the  same- 
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order  equation  relating  the  Jet  vortex  strength  and  the  dovnvash 
on  the  Jet.  The  Integral  equation  for  the  dovnvash  on  the  airfoil 
and  Jet,  calculated  by  the  Blot-Savart  Lav,  or  alternately  considered, 
the  mixed  boundary- value  problem  in  the  main- stream  velocity 
perturbations,  is  solved,  resulting  in  a  singular  integral 
equation  relating  the  Jet  vortex  strength  and  the  dovnvash  on  the 
Jet.  Boundary  conditions  on  the  Jet  ordinate  and  slope  at  the 
trailing  edge  are  specified.  A  proof  due  to  Spence  (I96IB)  -  to 
be  referred  to  as  III  -  is  given  that  the  potential  difference 
across  the  Jet,  at  a  given  instant  of  time  after  initiation  of  the 
unsteady  motion,  vanishes  if  a  point  sufficiently  far  dovnstream 
is  considered.  This  condition  is  shewn  to  be  stated  alternately 
in  the  form  of  the  Vagner  integral  condition  of  the  classical 
theory  of  unsteady  thin  airfoils  without  Jets.  Vith  this  condition 
the  formulation  is  completed  and  seme  properties  of  the  equations 
of  the  problem  are  discussed. 

Expressions  for  the  lift  and  pitching-moment  coefficients 
are  derived  in  Chapter  3  by  relating  the  pressure  distribution  on 
the  airfoil  to  the  vortex  distribution  representing  the  Jet.  These 
expressions  are  compared  to  their  counterparts  in  the  classical 
unsteady  thin- airfoil  theory. 

Detailed  equations  for  the  application  of  this  model  to 
five  fundamental  problems  of  Hat- plate  airfoils  are  given  in 
Chapter  b,  and  are  applicable  for  time- dependence  of  both  transient 
and  steady- state  oscillatory  nature.  These  problems  are; 

i)  An  airfoil  aligned  with  the  free- stream  direction  and 
having  a  time- dependent  Jet  deflection  angle  at  the  trailing  edge, 
the  "Jet- deflection"  or  Singular- bloving"  problem. 


ii)  An  airfoil  performing  a  purely  plunging  motion  about 
a  mean  position  aligned  with  the  free  stream,  the  jet  always  re¬ 
maining  tangential  to  the  trailing  edge. 

iii)  An  airfoil  performing  a  pitching  motion  about  some 
axis  in  its  plane,  the  Jet  being  always  tangential  to  the  trailing 
edge. 

iv)  An  airfoil  having  time- dependent  deflection  of  a 
mechanical  flap,  over  which  a  jet  is  blown  from  the  hinge  point  and 
leaves  tangentially  at  the  trailing  edge,  the  so-called 
"blown  flap"  or  "Jet-augmented  flap"  problem. 

v)  Entrance  of  the  airfoil,  previously  aligned  with  the 
free  stream,  into  a  sharp- edged  gust  of  constant  upwash  amplitude. 
The  jet  remains  tangential  at  the  trailing  edge,  and  the  relative 
speed  between  the  airfoil  and  gust  is  arbitrary,  but  constant. 

These  problems,  being  linear,  may  be  superimposed  in  any  desired 
fashion. 

The  reduction  of  the  equations  of  Chapters  3  and  4  to  the 
case  of  airfoils  without  jets  is  discussed  in  Chapter  5*  The 
classical  unsteady  theory  of  thin- airfoils  is  then  outlined, 
bringing  out  certain  features  of  the  flow  pattern,  e.g.,  the 
downwash  distribution  behind  the  airfoil,  which  have  particular 
importance  in  understanding  the  extension  of  this  theory  to  the 
Jet- flap  case.  The  steady  limit  of  the  equations  is  also  found, 
and  Spence's  steady-state  solutions  of  I,  (1958);  and  II  discussed. 

In  Chapter  6,  the  "boundary- layer"  nature  of  the  equations 
for  small  values  of  the  jet  momentum  and  small  times  after 
initiation  of  transient  motion  or  high-frequency  steady- state 
oscillations  sure  investigated  and  the  equations  are  -given  in  terms 


of  new  "boundary- layer"  coordinates.  All  the  airfoil  motion 
problems  -  ii)  to  v)  above  -  are  shown  to  reduce  to  the  same 
equations for  a  first  approximation  in  these  coordinates.  The 
relation  of  the  first  approximation  solution  of  this  airfoil 
motion  problem  to  that  of  the  Jet- deflection  problem  is  shown. 

The  further  small-time,  or  high-frequency,  approximation 
of  neglecting  X-  derivatives  with  respect  to  -t  -  derivatives 

is  then  examined  closely  in  Chapter  7>  Errors  made  in  III  in  the 
solution  of  the  jet-deflection  problem  in  terms  of  the  Jet  ordinate 
and  jet  vortex  distribution  for  small  times  are  pointed  out. 
Attempts  to  correct  these  errors  by  considering  the  downwash  on  the 
Jet  and  the  Jet  vortex  distribution  as  unknowns  and  a  similar 
approach  to  the  problem  of  airfoil  motion  demonstrate  the  failure 
of  this  approximation  to  give  valid  solutions.  However,  for 
high-frequency  steady- state  oscillations,  tentative  solutions  are 
proposed  for  both  the  jet- deflection  and  airfoil- motion  problems 
to  replace  the  erroneous  ones  given  in  III. 

Chapter  8  briefly  gives  seme  conclusions  of  this  research 
and  points  out  seme  areas  of,  and  approaches  to,  the  problem  of  the 
unsteady  motion  of  Jet- flapped  airfoils  where  further  work  is 
necessary.  These  are  felt  worthy  of  further  research. 


CHAPTER  1  -  DYNAMIC  INTERACTION  OF  THE  JET  AND  MAIN  STREAM 


1 . 1  Assumptions 

For  the  purposes  of  this  section,  no  assumptions 
about  the  airfoil  need  be  made,  except  that  it  has  a  Jet  emerging 
at  tne  trailing  edge,  and  that  it  has  seme  means  of  causing 
unsteady  motion  of  the  jet.  This  means  might  be  its  own  motion 
or  the  motion  of  the  jet  ducting  within  the  airfoil. 

The  jet  is  assumed  to  be  fully  developed  in  length, 
i.e.,  it  extends  infinitely  far  downstream  at  any  instant  its 
motion  is  being  considered.  The  momentum  flux  of  the  jet  at  the 
airfoil  trailing  edge  is  taken  constant  in  time.  The  flow  in  the 
Jet  is  assumed  inviscid,  incompressible  and  irrotational . 

It  is  also  assumed  that  the  local  velocity,  /U”  ,  of 

the  jet  is  very  much  greater  than  the  local  velocity,  LL  ,  of 

the  main  stream  in  the  vicinity  of  the  Jet.  This  velocity,  U.  , 
is  composed  of  the  undisturbed  free- stream  velocity  at  infinity 
upstream,  Uo  ,  plus  perturbations  due  to  the  interaction  between 
the  two  flows,  including  the  velocity  of  downward  translation  of 
the  Jet  boundaries.  For  the  practical  Jet- flap  applications  which 
have  been  proposed,  1/'  »  ,  so  the  perturbations, 

in  particular  those  normal  to  U.  ,  i.e.,  due  to  the  downward 
translation  of  the  Jet  boundaries,  must  be  very  small  compared 
to  V'  .  For  V  »  U.  ,  it  follows  that  the  instantaneous 
streamlines  of  the  flow  in  the  Jet  at  the  boundaries  between  the  Jet 
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and  main  stream  are  substantially  parallel  to  the  instantaneous 
shape  of  these  boundaries.  That  is,  with  reference  to  Figure  1, 
the  jet  velocities,  'Ui  and  Vx.  ,  at  the  boundaries  are  assumed 
to  be  parallel  to  the  boundaries.  Furthermore,  if  the  jet  is 
assumed  thin,  significant  mean  properties  of  the  flow  variables 
at  any  position  along  the  Jet  may  be  defined,  and  the  streamlines 
of  the  Jet  flow,  in  addition  to  being  parallel  to  the  boundaries 
at  the  boundaries,  will  be  all  parallel  to  each  other,  hence 
concentric  at  any  position  along  the  Jet.  This  property  leads  to 
a  great  simplification  in  the  equations  governing  the  problem,  as 
will  be  seen  in  the  next  section.  To  clarify  this  point,  consider 
the  excluded  case,  (JL  -  O(V)  ,  due,  say,  to  large  downward 
velocities  of  the  jet  boundaries.  The  instantaneous  streamline 
pattern,  even  for  a  thin  Jet,  would  intersect  the  Jet  boundaries  at 
an  appreciable  angle,  as  shown  below.  Such  a  geometry  must  be 
excluded  in  order  to  formulate 
a  tractable  model. 

The  restriction,  /1T  »  U.  ,  will  be  assumed 
to  be  met  in  all  types  of  airfoil  and  jet  motions  to  be  considered 
in  Section  2.2.  In  order  to  actually  attempt  solutions  of  these 
problems  in  Chapter  7>  discontinuous  motions  of  the  airfoil  and 
jet,  as  represented  by  the  unit-step- function  (equations  2- 1 T) > 
will  be  considered.  These  motions  will  always  be  treated  in  the 
sense  that  they  are  mathematical  idealizations  of  continuous 
motions  with  /VT>  » 


10 


The  main- stream  flow  is  also  assumed  inviscid,  incompressible, 
and  irrotational.  The  total  pressure  in  the  Jet  is,  for  dT  »U.  } 
greater  than  that  in  the  main  stream.  This  and  the  continuity  of 
static  pressure  across  each  boundary  between  the  flows  require, 
in  the  absence  of  viscous  mixing  (since  both  flows  are  assumed 
inviscid),  a  vortex  sheet  at  each  boundary  to  satisfy  the  velocity 
discontinuity  there. 

In  the  main  stream,  the  integrated  form  of  Euler's  equation 
is  the  so-called  "unsteady  Bernoulli  equation,"  which  is,  here, 

jf>  + 

(i-0 

where  ^  ,  <f>  ,  and  UU  are  the  local  static  pressure, 

velocity  potential  and  velocity,  respectively;  po  is  the  (constant) 
density.  f>o  and  U0  are  the  static  pressure  and  velocity, 
evaluated  at  the  undisturbed  conditions  at  infinity,  hence 
independent  of  time. 

In  the  Jet,  similarly, 

*>  p,  ^  , 

(1-2) 

where  fa  ,  4>y  ,  and  /V"  sure  as  above;  is  the  (constant) 

density,  in  general  different  frcm  p*  •  )».  and  Vo  sure  the 

static  pressure  and  velocity  evaluated  at  the  undisturbed  conditions 
at  Infinity  downstream,*  also  independent  of  time.  The  total 

*  See  Section  2.6  for  remarks  about  this  assumption  of  undisturbed 
conditions  in  the  Jet  at  infinity  downstream. 
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pressures  in  the  tvo  flours,  i.e.,  the  right-hand  sides  of  (l-l) 
and  (1-2),  are  in  general  different. 

The  condition  of  i r rotati  onali ty  in  the  jet  can  he 

written 


w  +  2: 

r 


o, 


where  Y  is  the  radius  of  curvature  of  the  streamlines, 
immediately  integrates  to 


(’-3) 

This 


f'lT  =  constant. 


1 .2  Pressure  Difference  Across  the  Jet 

Consider  an  incremental  element  of  jet  as  in  Figure  1 , 
described  by  its  position,  ^  ,  downstream  from  the  trailing 

edge  at  an  instant  of  time,  *£  .  The  subscripts  J  and  Z. 

refer  to  the  upper  and  lower  boundaries,  respectively,  of  the  Jet. 
Having  assumed  concentric  streamlines  within  the  jet,  as  discussed 
in  the  preceding  section,  the  element  may  be  described  further  by 
a  radius  of  curvature,  ,  to  the  centerline  of  the  Jet,  a 

thickness,  S(St4 r)  ,  such  that  Rfc,  i.)  ^  are  the 

upper  and  lower  boundaries,  respectively,  and  an  Incremental  angle, 
,  (positive  counter-clockwise)  subtended  by  the 
Jet  element.  The  pressure  difference  across  the  Jet  is,  frcm  (1-2), 

h-lt.i)  =  -fr 

_  .p,  lS£JtigKll£  fats.t)  -  VL/s,*)] 


0-5) 
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By  the  arguments  of  Section  1.1,  the  radii  of  curvature 
are  normal  to  the  streamlines  Inside  the  jet,  so  they  are  equi- 
potentials  of  the  Jet  flow.  Therefore, 

-  brjs*)  =  O  . 


The  mean  velocity  of  the  jet  may  be  defined  as 


(1-6) 


Vfs.i'i  =  ' 


(1-7) 


The  irrotational  condition  within  the  jet,  (1-4),  becomes 


[RM- 

(1-8) 

and  with  (1-7),  solving  (1-8)  for  the  velocity  difference, 

•V i t$fb)  -  lS,&)  ,  gives 


-vus,fi  = 

(1-9) 

Substituting  (1-6),  (1-7),  and  (1-9)  into  (1-5)  gives 
the  pressure  difference  across  the  jet  as 


f>,(r.t)  -  fh.  Is,*.) 


R/s.4) 

_  jferi 
Ris.t)  , 


(1-10) 
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where  the  Jet  momentum  flux,  ,  is  defined  as 

JM  « 

(1-11) 

Therefore  it  has  been  found  that  the  pressure  difference 
across  the  Jet  is,  as  in  Spence ' s  steady  formulation  of  I,  proportional 
to  the  Jet  momentum  flux  and  inversely  proportional  to  the  radius 
of  curvature  of  the  Jet.  Here  the  pressure  difference  is  a  function 
of  time  as  well  as  position  along  the  Jet.  The  essential  simplification 
of  the  vanishing  of  the  velocity  potential  difference  across  the  Jet 
has  resulted  from  the  assumption  of  Section  1 . 1  that  AT  »  U. 

1 . 3  Velocity  Difference  Across  the  Jet  in  the  Main  Stream. 

By  treating  the  flow  in  the  main  stream  and  relating 
it  to  that  in  the  Jet  by  the  condition  of  constant  static  pressure 
across  the  boundaries,  the  obIdp stream  velocity  difference  across 
the  Jet  can  be  found.  This  velocity  difference  is  required  to 
determine  the  strength  of  the  vortex  distribution  necessary  to 
represent  the  effects  of  the  Jet  on  the  main  stream. 

Considering  again  the  Jet  element  of  Figure  1,  equation 
(1-1)  evaluated  in  the  main  stream  across  the  Jet  is 

rn-fo  [a, ft*)  -  Uaft fl] 

”  f8  ^  £<&  fci:)-  <&.&*]]  # 


(1-12) 


IV 


The  mean  velocity  of  the  main  stream  across  the  Jet 
may  be  defined  as 

{J(stD  X  +4il$£) 


Substituting  ( 1 - 1 3)  into  (1-12)  and  solving  for 
the  velocity  difference,  tL(S/4)-  {XxlS,-L)  ,  gives 


(1-13) 


ULjs.ti-LLjSt) 


fo  U(s.i) 


Ul>',h)  &t  “  Qd5'^ ' 

(1-14) 


Since  the  static  pressure  is  continuous  across  the 
boundaries  between  the  flows,  the  pressure  difference  found 
considering  the  jet  alone,  (1-10),  may  be  substituted  into  (1-1 4) , 
giving 

u.fs.-t)-  a*/s.4)  *  . 


(1-15) 

This  result  differs  significantly  from  the  steady  case 
of  I  by  the  time  derivative  of  the  velocity- potential  difference. 

1 .4  Velocity  Induced  in  the  Main  Stream  by  the  Jet 

considering  the  elemental  vortices  comprising  the  Jet 
boundaries,  the  velocity  induced  in  the  main  stream  by  the  Jet 
can  be  calculated  using  the  Biot-Savart  Lav.  The  strength  of  these 
vertices  is  related  to  the  main- stream  velocity  at  the  Jet  boundaries. 
For  convenience,  the  complex  velocity  at  a  field  point,  ^  , 

in  the  complex  plane,  due  to  the  element  of  jet  at  £  will 

be  treated.  The  coordinates  are  as  follows: 


15 


The  velocity  induced  by  the  pair  of  elemental  vortices,  eLP,fe,4) 
and  tilfWSit)  >  is 


_  UnwbM +  o/i!- 1 

anrlys)  "  -fnvj-j)4  . 

(1-17) 

ais  nay  be  recognized  as  the  increment  of  velocity  induced 
by  a  vortex  of  strength 

JJYS/fe)  *  fil/Vs.i)  +d/ 9 

(1-18) 


and  a  doublet,  directed  downstream,  of  strength 
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* 


(1-19) 

both  of  which  lie  on  the  Jet  center  line,  plus  higher- order 
singularities  whose  strengths  are  higher- order  in  Sls.i)  . 

Define  a  vortex  strength  per  unit  length  along  the 
center  line,  Xx(S,t)  ,  and  a  doublet  strength  per  unit  length 
along  the  center  line,  yUj  ls,h)  ,  by 

JlFIS,1)  =  klS,i)  A8lS,S 

r 

d jk  1st)  *  } hist)  Mst)  Mist) J  # 


(1-20) 

The  strengths  of  the  elemental  Vortices  iP.IsM) 
and  dPxist)  may  be  written  in  terms  of  the  velocity  they 

induce  in  the  main  stream  at  the  boundaries,  namely. 


Mist)* 

d  list)*  -  +  £^£]  Mist) 


(1-21) 

Rewriting  (l - 18)  and  (1-19)  using  (1-20)  and  (1-21 )  gives 

fUst)M(s.i)  =  U,fc,4)[k/s,t>-  -Uxlstfakti  4  ^^]A«S,4) 


and 

Wi)  WISH)  =  -  jii, A9 


or 
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MA  •  [*«H  -«**)]- 

Slsti)  4 pl.fr, 4) 


Substituting  (1-1 3)  and  (1-15)  into  (1-22)  gives 

**■«•  -  uks 


(1-22) 


and 


(1-23) 

The  influence  of  the  jet  on  the  main  stream  can  thus 
be  calculated  by  replacing  the  jet  by  vortex  and  doublet  distributions 
along  the  Jet  center  line,  their  strengths  being  given  by  (1-23). 

In  order  to  make  precise  the  neglect  of  the  higher- order  singularities 

$ 

in  ( 1 - 1 7) >  and  to  simplify  the  strengths  of  the  singularities  in 
(1-23),  the  limiting  case  of  a  thin,  high- speed  jet  will  be  treated 
in  the  next  section. 

1 .5  Limiting  Case  of  Thin,  High-Speed  Jet 

$y  order  of  magnitude  considerations,  the  limiting  case 
of  a  thin,  high-speed  jet  of  constant  momentum  flux,  3*  ,  will 

be  deduced. 

The  (constant)  mass  flux  of  the  Incompressible  jet  is 

defined  as 


m  =  f, Vis&SM  , 


(1-24) 
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so  the  jet  momentum  flux  may  be  written  using  (1-1 1)  and  (1-24)  as 


J  (S,h)  *  # 


(1-25) 

Writing  the  subscript  (  )#  to  refer  to  conditions  at 
infinity  downstream  where  all  the  quantities  are  assumed  constant 
in  time,  (1-25)  may  be  written  between  seme  jet  position,  5  , 

and  infinity  as 


J(s,t)  _  ,  .  V(s,«-V4 


.  + 

Vo  > 


(1-26) 

where  &WS,4)  »  V/S,t)  -  Vo  ,  and  similarly  AUlS,±)  a  UlS^-Vo 

are  the  perturbations  of  UlS,4;)  and  WS,t)  from  the  conditions 
at  infinity.  It  is  important  to  note  here  that  these  perturbations 
are  not  necessarily  small  for  the  following  analysis.  They  are 
written  this  way  for  convenience. 

Equations  (1-1)  and  (1-2)  may  be  written,  between  an 
arbitrary  field  point  in  or  near  the  jet  and  the  quiescent  conditions 
at  infinity,  as 


-  N  +4£-  [vVs,0  -  \Sj  +  f r  ^  =  O  " 

. 

pW-)h  +4[oYc.«-Li]  +  f.  ^  *  o 


0-27) 
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For  an  order- of- magnitude  analysis  It  la  sufficient  to 
consider  the  pressure,  pi*, 4:)  ,  in  each  of  these  equations 

as  the  same,  say  an  average  pressure  in  the  neighborhood  of  the 
jet.  Eliminating  —  f),  from  this  pair  of  equations 

cud  solving  for  A  Ws, t)  ,  gives 


*  2  iV,+iV/S,t)  2  ZV.tiVa.t) 


(1-28) 


To  estimate  the  order  of  the  time  derivatives  of  the 


velocity  potentials,  write 


a. 

~*ft  TT~  . 


and 

**  a*.  . 

In  the  main  stream,  a  characteristic  time- dependent 
would  be  the  product  of  the  characteristic  perturbation  speed  2  VIS,-L) 
and  a  characteristic  length,  say  the  chord,  C  .  Likewise  a 
characteristic  time,  A  t  ,  would  be  the  chord,  C  ,  divided  by 
the  free- stream  speed,  u.  .  Thus, 


or  C&UIS& 

VUo 


UAU&fc) . 


(1-29) 


In  the  Jet,  a  similar  analysis  would  use  the  jet  speeds, 


&V(S,i)  and  ytS.-k)  f  and  the  characteristic  Jet  dimension, 
$  ^S,4)  .  Thus, 


“7rft  -fisqfir 


V*AWs,tf 


(1-30) 
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Substituting  (1-29)  and  (1-30)  into  (1-26)  gives 


(1-32) 


*  In  the  case  of  a  Jet  or  airfoil  deflection  having  a  unit- step- 
function  time  dependence,  f*^*  ’  aa^L  c°nse4uently  2  U/J.-fe)  , 

would  have  in  the  first  Instant  the  infinity  of  a  Dirac  delta  function, 
in  clear  violation  of  the  assumption  AU^WO(U,).  However,  as 
discussed  in  Section  1-1,  unit- step  functions  are  considered  because 
of  their  mathematical  convenience,  and  may  be  considered  as  the  generali¬ 
zation  of  a  physically  realistic  deflection  of  finite  rate,  where 

OMo)  This  is  the  same  consideration  made  in 

classical  linearized  unsteady  airfoil  theory,  where  the  infinities 
Introduced  by  the  derivatives  of  the  step  function  clearly  violate 
the  small-perturbation  assumption,  unless  understood  in  the  above  sense. 


L 
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Substituting  (1-32)  into  (1-26)  gives,  finally, 

“-'♦"S'. 


Jfri) 

*Jo 


(1-33) 


In  most  practical  Jet- flap  applications,  the  Jet  speed, 
Vo  ,  is  very  much  greater  them  the  free- stream  speed,  U» 


or 


ik 

Vo 


>0). 


Consider  now  the  limit  as  vanishes,  such  that  the 

Jet  momentum  flux  remains  finite.  In  this  limit  Vls.i)  must 
become  infinite  although  the  flow  is  still  considered  incompressible. 
In  detail,  from  (1-33)> 


37s,6)  *  J©  »  =  T  * 


and  as 


and 


£  (S,b)  vanishes  , 

Vts.i)  ~  [£/s,*)]’/x 

VV1  a.  £  £/S,fc)] 


0-34) 


(1-35) 

The  important  relations  from  the  earlier  analysis  may  be 
written  in  terms  of  this  approximation,  treating 
as  the  small  parameter  tending  to  zero.  Equations  (1-10),  (1-1 5), 
and  (1-23)  become 


J 

KM)  , 


(1-36) 

*"T*  a 


(1-37) 
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y,  1st)  S  — ^  » 

(1-38) 


yUT^»*  Ott). 


(1-39) 

Die  jet  vortex  strength  is  zeroeth  order  in  the  jet  thickness, 
iCst4s)  ,  and  becomes,  to  this  order,  equal  to  the  main-stream 
velocity  difference  across  the  jet.  The  doublet  strength  is  of  first 
order  in  jet  thickness  and  will  thus  bemeglected  in  the  lifting 
problem.  Therefore  the  lifting  problem  is  independent  of  thickness 
in  this  approximation.  Furthermore,  the  neglect  of  the  higher-order 
singularities  in  (1-17)  is  justified  by  this  limiting  analysis,  as 
they  are  higher- order  in  the  jet  thickness. 

Therefore,  the  non- homogeneous  flow  problem  of  a  thin, 
high-speed  jet  embedded  in  the  main  stream  has  been  approximated, 
representing  the  Jet  by  a  vortex  distribution  along  the  Jet  center 
line.  This  vortex  distribution,  whose  strength  is  given  by 
(1-38),  interacts  dynamically  with  the  main  stream  in  the  same 
manner  as  the  Jet,  in  the  limit  of  a  vanishingly  thin,  constant- 
momentum- flux  jet. 


CHAFFER  2  -  DERIVATION  OF  TRE  BASIC  EQUATIONS  OF  THE  LINEARIZED  UNSTEADY 

PROBLEM 


2.1  Linearization  of  the  Problem 

No  treatment  of  the  airfoil  vas  required  in  Chapter  1 , 
except  indication  that  It  vas  the  source  of  the  Jet  and  the  cause 
of  its  unsteady  motion.  In  principle,  then,  a  theory  with  no 
restrictions  on  the  airfoil  or  its  motion  could  be  developed  using 


the  model  of  Chapter  1  as  a  representation  of  a  thin,  high-speed  jet- 
flap,  so  long  as  none  of  the  assumptions  of  Section  1 . 1  were  violated. 
Practically,  it  is  desirable  to  make  the  thin- airfoil  assumptions 
and  linearize  the  problem. 

Assume  that  the  airfoil  is  thin  and  symmetrical  about 
a  mean  camber  line,  and  has  chord,  C  .  Choose  %■-  and 
coordinates,  so  that  %  is  the  free- stream  direction  and  is 
positive  downwards,  for  convenience.  If  )//%,€.)  is  the  ordinate 
of  the  mean  camber  line  for  O  %^-C.  ,  and  is  the 

ordinate  of  the  Jet  center  line  for  oo  ,  and  if  , 

0  £  %<  c  ,  is  the  thickness  of  the  airfoil  (the  Jet  having 


been  assumed  of  zero  thickness  in  Chapter  1 ) ,  the  problem  may  be 
linearized  by  assuming  >  ~ * 

and  ^  to  be  small  quantities  with  respect  to  unity.  The 

velocities  induced  due  to  these  quantities  are  then  small  compared 
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to  the  free- stre«»  velocity  at  infinity,  U#  ,#  and  the  squares 
of  these  quantities  nay  be  neglected.  Consistent  with  this,  the 
airfoil 


and  Jet  are  considered  to  lie  along  the  axis,  regardless 

of  their  actual  shape,  and  the  unsteady  kinematic  boundary  condition 
that  they  be  instantaneously  streamlines  is  to  be  satisfied  on  the 
•TL-  axis,  too. 

As  is  well  known  for  linearized  flow,  cf . .  Robinson  and 
Laurmann  0956);  ppe.  129  and  170,  the  thickness  and  lifting  problems 
decouple  and  may  be  considered  separately.  For  the  lifting 
problem,  the  airfoil  may  be  represented  by  a  vortex  distribution,  sub 
the  Jet  has  been.  3he  above  kinematic  boundary  condition  that  the 
airfoil  and  Je  surfaces  be  streamlines  is  satisfied  if  t he  dowmash, 

UT  (%,  -t)  ,  is  given  by  the  linearized  convective  derivative  of  the 

airfoil  or  Jet  ordinate, 

imma- 

(2-1) 

as  given  in  Robinson  and  Laurmann  (1956),  p.  3. 

*  Here,  as  in  Section  1 .3,  infinities  implied  by  the  step- function 
time  variation  must  be  considered  as  the  limit  of  finite-rate 
processes  which  do  not  violate  the  isn  perturbation  approximation. 
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Consistent  with  the  linearization;  the  local  mean  main* stream 
speed  across  the  jet,  0iS,4s}  ,  given  by  (1-13)>  nay  be  taken  to 

be  Uo  everywhere.  The  %r  and  -  components  of  the 

local  stream  velocity  may  be  written  as  Uo  +  Uj%t4rS  and 
UTiJLf-k)  respectively;  where  U}UT  are  to  denote 

the  velocity  perturbations  for  the  remainder  of  this  report.  Ho 
confusion  with  their  previous  usage  should  occur. 

Define  ^  of  a  quantity  to  mean  the  difference  in 
its  value  Just  below  the  axis  from  Its  value  Just  above;  i.e., 


&(  k*(  )r~-(  )r-. 

The  coordinate  S  ,  denoting  distance  along  the  Jet  from  the 
trailing  edge  linearizes  to  %-  C 

The  curvature  of  the  Jet  is  given  geometrically  by 

&  "  -0  +  ,  c^°°’ 

which  is;  linearized; 


i  _  _ 

Rfs,«  >%*• 


c^x.^00 ' 


If  the  Jet- momentum- flux  coefficient  is  defined  as 


(2-3) 


Ct  * 


T 


(2-4) 


the  pressure  difference  across  the  Jet;  (1-36);  and  the  Jet  vortex 
strength;  (1-38);  may  be  written  in  linearized  form  as 
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(2-5) 

and 


*/*,«  =-i  U.tC»  - -|j  ^ 

(2-6) 

where  the  subscript,  (  )j  on  the  vortex  strength  has  been 
emitted  for  simplicity. 

From  the  definition  of  the  velocity  potential,  and  the 
strength  of  the  vortex  distribution  in  terms  of  the  free- stream 
velocity  perturbation, 

%  Ht%jO+,4s)  ,  cu.%jcoo  „ 

(2-7) 

The  jet  vortex  strength,  (2-6),  written  using  (2-7),  gives  an 
equation  in  the  potential  difference  across  the  jet,  namely 

=  iggs  +Uo&gtt>  U?CCT 

(2-8) 

Taking  the  linearized  convective  derivative  of  this  and  using  (2-1) 
gives 

£**££>,  ggjgi  -£u:ur  ^ 

(2-9) 
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Differentiating  (2-8)  with  respect  to  X  and  using  (2-7) 


gives 

j£4L)U40. 


Fran  (2-9)  or  (2-10),  it  also  follows  that 


(2-10) 


tgg  +o:v*p,.  iutcc*  Ttegu  3 

(2-11) 

The  unsteady  Bernoulli  equation,  (1-1 ),  evaluated  across 
the  X*  axis  is,  with  (2-7), 

+  u  =  -AVp «  ,  o^~». 

(2-12) 

If  A  1^/%,-tr)  given  by  the  analysis  of  the  Jet,  (2-5),  is 
substituted  directly  into  (2-12),  (2-8)  seems  to  follow  immediately, 

Just  from  the  linearization  of  the  problem.  The  analysis  of  Sections 
1.3  to  1.5  is  necessary,  however,  to  Justify  the  representation  of  the 
Jet  by  a  vortex  sheet  as  being  a  valid  approximation  to  the  Jet  flow, 
i.e.,  that  the  Jet  thichness  is  really  a  higher*  order  effect  for  thin, 
high-speed  Jets. 

2.2  Nature  of  the  Motion  and  its  Time  Dependence 

The  unsteady  motions  to  be  treated  are  of  the  following 

types: 

i)  The  airfoil  is  fixed  and  the  Jet  deflection  angle, 

~CLh)  ,  relative  to  the  slope  of  the  airfoil  at  the  trailing 


SB 


edge  varies  time- dependent!/ .  Oils  motion  Is  given  by 


t  (h)*  t0  m) , 


(2-13) 


where  X#  is  the  small  amplitude  of  the  deflection.  Each 
point  on  the  jet  boundary  moves  normally  to  the  %-  axis  only. 

11 )  The  airfoil  Is  performing  some  small-amplitude 
time- dependent  motion  about  a  mean  position,  e.g. ,  plunging, 
pitching,  or  rotation  of  a  mechanical  flap  with  the  Jet  deflection 
angle  at  the  trailing  edge  fixed  relative  to  the  airfoil.  Such 
motion  can  be  characterized  by 


V(%)  F(&  ,  O-Mtcc  . 


(2-14) 

Each  point  on  the  airfoil  and  jet  boundary  moves  normally  to  the 
■%-  axis. 

ill)  The  airfoil  and  jet  deflection  angle  at  the  trailing 
edge  remain  fixed  and  a  sharp- edged  gust  passes  over  the  airfoil 
and  jet  from  the  leading  edge.  Such  a  gust  can  be  represented  by  the 
downwash  distribution 


W lr,h\  *  F(i-  *£)  # 


(2-15) 

where  Wg/%)  Is  the  small  amplitude  and  U»/X  is  the 
speed  relative  to  the  airfoil.  Again  the  Jet  boundary  moves  normally 
to  the  %-  axis  only.  Despite  the  vortex  sheet  required  at  the 
edge  of  the  gust,  this  problem  has  been  treated  in  classical  unsteady- 
airfoil  theory  as  if  it  were  a  potential  flour  with  an  Imposed, 
cf.  von  Harman  and  Sears  (193$)  and  Miles  (195^)* 
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The  Y 1%)  ,  \Mj/li)  ,  and  F/t”  functions 

are  assumed  to  be  of  a  form  such  that  in  any  case  the  downvaab  on 
the  airfoil  may  be  represented  by  a  power  series  in  %  .  Such 

a  power  series  is  equivalent  to  the  familiar  CHauert  Series  for 
the  downwash  on  the  airfoil,  Glauert  (1943),  hence  there  is  little 
loss  of  generality  in  the  choice  of  functions  which  can  be  considered. 

Two  types  of  time  dependence,  Ftzr)  ,  are  of 
particular  interest.  First  are  the  transient  problems  where 


=  f /«  1  ft) . 


(2-16) 


The  function  P/'t)  is  any  non-dimensional  function  of  time, 
and  d/t)  is  the  unit- step  function,  defined  by 


-lrt-1 


To  ,* 
*1  •  .* 


iO 


(2-17) 

The  fundamental  transient  problem  is  that  of  response  to  the  unit- 
step  function,  i.e.,  p/t)s  I  ,  since  the  response  to  all 

other  p/'t)  may  be  found  from  it  by  Duhamel  superposition.  Second 
are  the  cases  of  steady- state  oscillations,  where 

F/t)*  e.wt  , 

(2-18) 

The  motion  is  assumed  to  hare  begun  at  t  ■=  -  oo  f  yith  (R/w) 
being  the  circular  frequency  of  the  oscillations,  and  &/<**)  being 
the  damping  factor.  For  ,  the  oscillations  diverge 

exponentially  with  time,  but  began  with  zero  amplitude  at  t*-  •»  , 

and  are  thus  meaningful  for  t  .  The  9/uJ)> O  ,  i.e., 


9 


exponentially  duped,  oscillations  oust  be  excluded  since  they 
arise  from  the  physically  unrealistic  initial  condition  of 
infinite  amplitude  at  -t  *  -  ®°  .  Exponentially  damped 

oscillations  beginning  at  finite  amplitude  at  ■£  *  o  could  be 
treated  as  a  transient  problem,  with  f74)  *  §(u>)  >  o 

For  a  consideration  of  these  damping  considerations  In  classical 
unsteady- airfoil  theory,  see  the  paper  by  Luke  and  Dengler 
(1951)  and  the  subsequent  notes  by  Van  de  Vooren  (1952),  Lai tone 
(1952),  W.  P.  Jones  (1952)  and  Dengler,  Ooland  and  Luke  (1952). 
2.3  Downwash  Conditions  on  the  Airfoil  and  Jet 


Hfce  downwash  due  to  the  vortex  distributions  representing 
both  the  airfoil  and  Jet  may  be  found  by  using  the  Biot-Savart  Lav 
for  the  velocity  induced  by  an  incremental  vortex  element,  and 
integrating  over  the  distributions,  giving 

c 


urtx.tu 


VH.iUi 

i-% 


0*.  %  *■<*> 


(2-1 


This  1b,  in  fact,  two  equations,  one  for  o  *-X*-C  where  the 
first  integral  exists  in  the  sense  of  Cauchy's  Principal  Value,  the 
second  for  where  the  second  Integral  exists  In  the 

same  sense.  Whereae  (2-8)  to  (2-11)  express  the  dynamic  interaction 
of  the  Jet  and  main  stream,  (2-1 9)  expresses  a  kinematic  interaction. 

The  two  equations  (2- 1 9)  and  equation  (2- 1 l )  sure  three 
equations  in  the  four  functions  ur/1,4)  ,  O 

ur/t-,4)  ,  }  i%tb)  , 

0*.%*C  ;  X(%,k\  ,  C*-V-4.00 

Alternately,  combinations  of  and  )i ,  urty.,±) 

and  ,  or  >7%.  4)  and  might 


and 
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toe  considered,  also  giving  three  equation*  in  four  unknown  functions. 
In  the  cases  of  practical  interest;  the  airfoil  shape;  //%.,£) c 
and  from  (2-1)  then;  the  downwash,  Ur a*-K-c  ,  is  prescribed, 
and  along  with  the  boundary  conditions  to  toe  discussed  in  Sections 
2.5  and  2.6 ,  the  equations  (2-19)  and  (2-11)  are  then  three  equations 
in  the  remaining  three  unknown  functions. 

The  downwash  equations  (2-1 9)  may  toe  inverted,  using  the 
prescribed  airfoil  shape,  Yfx.fc'l  (  0  ,  in  a  way  to 

mathematically  decouple  the  airfoil  vortex  distribution,  , 

OC'f.e.c  ,  from  the  functions  describing  the  Jet.  That  is, 

a  form  of  the  equations  may  toe  found  where  the  Jet  ordinate  or 
downwash  and  vortex  distribution  or  potential  difference  may  toe 
solved  for  independently  of  the  airfoil  vortex  distribution.  Upon 
solving  for  the  former,  the  latter  may  then  toe  evaluated. 

There  are  several  alternate  methods  of  inverting  (2-19), 
all  of  which  have  a  similar  mathematical  basis,  although  the 
individual  techniques  are  different.  Spence  used  a  "null- transform" 
technique  (icf.  Section  D12  of  Heaslet  and  Lomax  (1954)  )  in  I. 

In  II  he  used  a  generalization  of  a  result  first  given  toy  Carleman 
(1922)  for  the  Inversion  of  a  certain  singular  Integral  equation. 

For  the  present  report,  the  technique  put  forth  toy  Cheng  and  Rott 
(1954)  is  used.  This  is  based  an  the  solution  of  the  equivalent 
mixed  boundary- value  problem  in  the  velocity  perturbations. 

The  downwash  equations  (2-19)  will  not  toe  considered 
directly,  tout  the  corresponding  mixed  boundary- value  problem  in  the 
complex  velocity  function,  is  treated. 

In  this  form,  boundary  values  are  given  in  terms  of 

and  UT on  the  alternate  segments  of  the 
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%-  axis  corresponding  to  the  airfoil,  Jet  and  the  region 
upstream  of  the  leading  edge.  From  the  strength  of  the  vortex 
distribution  representing  the  jet,  (2-7),  and  from  the  downwash 
required  on  the  airfoil,  the  boundary* value  problem  here  Is  as  follows: 

tLfco*,t)»o  {  U.  hi,  <>■*•,  fcW  _ 


LU^ot.h)  =  O  t  —  «o  «c  o  n 
Ur/^o+. b)  *  &  ,  O  <oc*  C 

UilX'Ot.l)  =  XlX,*),  , 


> 


(2-20) 

where  UH%,b)  *■  ,  a  x  %.*■  C  is  known,  but  ,  C£-%x.co 

is  still  unknown  and  leads  to  another  integral  equation  after 
inversion. 

To  solve  this  problem,  the  corresponding  homogeneous 
solution  in  the  complex  velocity  is  found,  i.e.,  for  the  problem 
where 

tin  /ijO+.b)  =  0,-<»  £.%*-0  " 

UTH  1% O  ,  ©k.%xc  ■ 

U.H  /*,*►»  »  O  ,  X^ooJ  # 


(2-21) 

This  homogeneous  solution,  Uo(^,h)  ,  say.  Is 

constructed  by  considering  the  behavior  required  of  the  inhomogeneous 
solution  at  the  edge  points,  i.e.,  the  leading  and  trailing  edges, 

Ot  *  O  and  Tt »  C  .  Here  the  usual 


\ 
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singularity  is  permitted  at  the  leading  edge  O  ,  while 

at  an  singularity  is  not  penaitted 

since  the  jet  emerges  here,  preventing  flow  around  the  trailing 
edge. 

In  the  classical  non- jet- flap  case,  flow  around  the 
trailing  edge  is  excluded  by  the  Kutta  condition.  With  the  jet-flap 
a  stronger  condition  is  imposed,  since  the  angle  of  deflection  of 
the  jet  relative  to  the  slope  of  the  airfoil  at  the  trailing  edge, 
X£4r)  >  may  also  be  prescribed.  Anticipating  later  results,  it 

is  necessary  to  remark  that  exclusi'm  of  the  1 /square- root  singularity 
does  not  prevent  a  logarithm! c  singularity  at  the  trailing  edge. 

Such  a  singularity  is  weaker  and  represents  not  flow  completely 
around  the  trailing  edge,  but  only  deflection  of  the  flow  by  some 
angle.  This  corresponds  to  a  discontinuity  in  the  downwash  at  the 
trailing  edge.  In  a  similar  fashion',  logarithmic  singularities  on 
the  airfoil  surface  sure  permitted  surd  will,  in  fact,  result  from 
prescribed  discontinuities  in  the  downwash.  The  classical  example 
of  such  a  singularity  is  that  at  the  kink  of  a  bent  flat  plate, 
first  pointed  out  by  Glauert  (1927)  -  cf.  also  Spence  (1958)  - 
in  his  representation  of  an  airfoil  with  flap. 

To  make  the  condition  at  the  trailing  edge  more  precise, 


write 


0. 

(2-22) 


Since  no  singularities  stronger  than  1  /square  root  are  permitted 

in  the  velocity  on  the  airfoil,  the  circulation  around  the  airfoil, 

1-'.,  lim  &<pl%,4r)  ,  la  regular,  and  since  its  time 

%+c 

derivative  must  be  also, 
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Prom  (2-22)  and  (2*23),  therefore, 

+u.  34^>]  =  o , 

(2-24) 

and  using  (2-8),  therefore 


•jt-e-c 


O  . 


She  homogeneous  solution  Is,  then,  from  Cheng  and 
Rott  (1954), 


(2-25) 


(2-26) 


The  complex  velocity  function,  F/j,  h)  ,  say,  which  Is 
the  quotient  of  the  Inhomogeneous  and!  homogeneous  solutions. 


(2-27) 


has  the  property  that  Its  Imaginary  part,  l.e.,  Its  downwash,  Is 
known  everywhere  along  the  •Jl-  axis.  Thus  the  problem  of 
finding  Is  the  so-called  direct,  or  thickness  problem 

of  thin- airfoil  theory,  whose  solution  can  be  written  down  immediately 
by  considering  the  velocity  induced  by  a  source  distribution  whose 
strength  Is  twice  the  downwash.  This  technique  will  become  clearer 
upon  working  through  the  problem. 
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Tor  reference,  write  out  the  complex  velocity 
components  of  the  three  Auctions, 

Hoff, ir)  *  u » lx, yd)  ti 


(2-28) 


In  order  to  remain  on  the  same  branch  of  the  function 
it  is  necessary  to  restrict  06  M3S*ir- 

Tor  -e»4)U.o  ,  »0+  ,  and  d^ig-cj  »7T,  d^/j|  *  7T, 


from  (2-26), 


2  1^1  ,  -oo*%<iO. 


(2-29) 


From  (2-20),  (2-27),  (2-28),  and  (2-29) 


F/»,o+,*W  (’(j-.o+M+iftwi)  -  -jy^rgnl^ 

(2-30) 

and  so  the  imaginary  part  of  PY%,o+,t)  is 

h)  *  O  ,  -  *  X.*  O  . 


(2-31) 


In  a  similar  fashion,  for  o^-)L4-  i  ,  , 

and  OA^lycl  »7T  ,  °  /  fran  (2-19)> 

-  -  (SFT.  . 


(2-32) 


From  (2-20),  (2-27),  (2-28),  and  (2-32), 


3  6 


Flx,0i,i)  «  J f  i gb.O+.h)  m  ^  O^JCAC, 


and  so  the  Imaginary  part  of  F(Vf,4) 

-(Af  “f*. 


(2-33) 


is 


(2-34) 


Finally,  for  c^-y.x.es^  ^sof  ,  and 
d^Jg-i  =  O  ,  from  (2-19), 

Ha^et^s  I  p?)4> 


(2-35) 

From  (2-20),  (2-27),  (2-28),  and  (2-35), 

s  J)(y.to+,i)  ¥i£(%to+,4s)  *  ~ ^  j  », 

(2-36) 

and  so  the  Imaginary  part  of  F/jt»o+,4i  is 

S'/a 


px.ot.b)  ,  d. ^y.^.00. 


(2-37) 


From  the  three  expressions  (2-31 ),  (2-34),  and  (2-37),  it 
Is  seen  that  the  imaginary  part  of  F/x.,<n-,<r)  ,  i.e.,  jj !%, 0 t) 
for  —  •<.  oo  is  known.  By  equation  (1)  of  Cheng  and  Rott 
(1 954),  which  Is  Just  the  velocity  induced  by  a  source  distribution 
of  strength 


> 
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or  using  the  definition  of  p(fiir) ,  (2-37), 

"ty*  '  ■ 

(2-38) 

With  (2-1),  (2-19),  (2-26),  (2-24),  (2-27),  (2-31),  end  (2-37), 

.  ,  ,x  i  f[JL^JZS£4*  i  l \(1  f* 

s  "  7f  (t/  I JJC'!'  *-*  *’8'  J 

(2-39) 

This  gives  the  complex  velocity,  'Wj/M  ,  everywhere  in  the 
field. 

To  get  the  desired  inversion  results,  in  addition  to 
recovering  the  boundary  values,  the  limit  ^  O  +  of  equation  (2-32) 
must  be  taken  on  each  segment  of  the  %-  axis.  The  results  are 


ATfr.,o+(i)  4  4r 


(2-40) 


04)t*  C, 


UhLiti  m  i. 


(2-41) 


(2-42) 


3B 


The  boundary  values,  (2-20)  are  apparent  in  (2-4o)-(2-42),  and  by 
using  (2-7)  and  (2-1),  the  desired  inversion  results  are 


(2-44) 

The  integral  equations,  (2-19),  for  the  downvash  have  been 


inverted  and  the  equations  (2-43)  end  (2-44)  are  the  resuJ-  3-  The 
airfoil  vortex  distribution  if ( %,-V)  ,  04|UC  ,  does  not 

appear  in  vhat  are  now  the  two  governing  equations  of  the  problem. 

The  problem,  now,  is  to  solve  (2-11)  and  (2-44)  for  the  downvash  on 
the  jet,  Ut"l%,4r\  ,  and  the  Jet  vortex  distribution,  $7%, 4}  , 
both  C  tfi  , using  the  boundary  conditions  to  be  discussed 
in  Sections  2.5  and  2.6.  Alternately  Y/%,t)  and  , 

iWly-.-b)  and  >  or  and  A$(y.,4i\  may 

be  considered  to  be  the  unknown  functions.  The  airfoil  vortex 
distribution,  /%,-(?)  ,  O  4  %.4-C  ,  may  then  be  evaluated 

from  (2-43),  as  could  the  donwash  upstream  of  the  leading  edge  from 

(2-2*0). 

For  future  reference,  equation  (2-43)  may  be  integrated 
with  respect  to  %.  ,  giving,  by  (2-7),  the  potential  difference 

across  the  airfoil.  Since  there  are  no  tt.-  velocity  perturbations 
on  the  >%.-  axis  upstream  of  the  airfoil  leading  edge,  (2-20), 

A4>(Oj-t)-=o  , 


(2-45) 
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so 


Ablrtf  =  J* Xit.4r)d3  t  o*  j 

V#  0  >  c 


0*~K*c. 

(2-46) 


To  simplify  this  expression,  interchange  the  order  of  integration. 
This  is  permissible  according  to  Section  D12  of  Heaslet  and  Lcmax 
(1954)#  since  the  behavior  of  Xlh-b)  at  the  trailing  edge, 
(2-22),  and  the  assumption  that  on  the  airfoil  can  be 

represented  by  a  power  series  in  %  eliminate  any  residual 
singularities  of  the  Integration.  This  gives 

0  6  c.  « 


O , 

(2-47) 


Using  the  results  of  Appendix  A,  (A-48)  and  (A-49),  the  potential 


difference  across  the  airfoil  is 

c 


d9 


a  Jc 

04.1LC.C. 

(2-48) 

In  particular,  if  %v  c -  ,  (2-48)  becomes  the  circulation  around 
the  airfoil, 


A4>fc-,±)  *  ru)  M  2.  ur/S^if  +  ■  ^j4-  . 


(2-49) 
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2.4  Identification  of  the  Airfoil- Quasi -Stsady  Terms 

The  first  tens  on  the  right-hand,  sides  of  (2-43),  (2-44), 
(2-49),  end  the  first  two  tens  on  the  right-hand,  side  of  (2-46)  will 
now  be  identified.  A  quasi-  steady  quantity  is  for  the  present  purposes 
defined  as  one  dependent  only  on  the  instantaneous  unsteady  notion, 
independent  of  the  previous  tine  history  of  the  flow.  Since  the 
tens  in  question  are  functions  only  of  the  airfoil  motion,  and  are 
Independent  of  the  Jet,  they  will  be  referred  to  here  as  the  airfoil- 
quasi- steady  terns.  They  are  calculated  such  that  the  instantaneous 
tine-dependent  boundary  condition  that  the  airfoil  surface  be  a  stream¬ 
line  is  satisfied,  but  neglecting  any  effects  due  to  shed  vortidty 
or  the  jet.  These  terms  are  results  of  the  classical  steady- state 
solution  for  a  thin  airfoil  with  tine-dependent  downsash  on  the  airfoil, 
i.e. ,  the  mixed  boundary- value  problem  treated  in  Section  2.3,  with 
U(%i0+,4)»  o  ,  C*-7-*-ae  f  or 

U  -  O  > 

ur(%l0%t)  -  0**.<c 

u  0+(t)sO  9  00  • 

(2-50) 

Since  the  same  requirements  are  made  on  this  solution  at  the 
leading  and  trailing  edges  (although  thqfclassical  Eutta  condition  is 
Invoked  at  the  trailing  edge  instead  of  the  argument  of  page  33), 
solution  follows  the  procedure  of  Section  2.3  directly,  aquations  (2-43) 
and(2-44)  with  fr/fc.t)  =0  ,  >  are 

then  the  airfoil- qiiasi- steady  values. 


9m  sixfoil-  janal-  steady  rortex  distribution  is. 


trm  (2-43), 

♦(suffer 


0<.%<C, 


(2-51) 


Trm  (2-48),  the  airfoil- quasi- steady  potential  difference  across 
the  airfoil  is 


and  the  airfoil-  quasi-  steady  circulation  is  then 

c 

Hli)~  A  h(c-, 4)  =2(1^)  Urff,4r)Jt  . 

9»e  airfoil-  quasi-  steady  dcwnwaeh  behind  the  airfoil 
is,  from  (2-44), 


04JUC, 


(2-5?) 


(2-53) 


£<£.%« OO 


(2-54) 

Since  the  airfoil  and  gust  Shapes  were  assuaed  in  Section  2.2  to  give 
power- series  representations  of  the  dcwnseah  on  the  airfoil,  it  can 


be  shown  that  this  downsash  is  continuous  at  the  trailing  edge  with 
the  airfoil- quasi-  steady  downsash  behind  the  airfoil.  9»  power  series 
Mgr  lie  written 


ar/*,4)  *  ,  o*-x*c 

n*o  ' 

(2-55) 

where  G*l^)  is  the  appropriate  function  of  tine.  Iraluating 
(2-54)  with  this,  using  (Ar-33) >  gives 


ur.lt, D  * 

(2-56) 

Therefore 

wr.^+,t)=  £  6h^6H  -  UIiC-,4,) . 

M*0 

(2-5T) 

Since  problems  with  the  airfoil  shape,  YlxM,  om-c  > 
prescribed  are  to  be  considered,  these  airfoil-  quasi-  steady  quantities 
may  be  calculated  Immediately  in  a  given  case. 

2.5  Boundary  Conditions  at  the  trailing  Edge 

There  are'  two  boundary  conditions  on  the  jet  ordinate  at 
the  trailing  edge  of  the  airfoil.  First,  the  Jet  must  be  continuous 
with  the  airfoil  there;  l.e. , 


Yfc+'h)  =  Y/d-i)  . 


(2-58) 

Secondly,  as  mentioned  in  Sections  2.1  and  2.3*  the  angle  of  deflection 
of  the  Jet  relative  to  the  slope  of  the  airfoil  at  the  t  railing  edge 
may  be  prescribed,  giving,  using  (2-13), 


2MiKm^+  r.F«). 


In  terns  of  the  dovnvash,  (£.58)  and  (2-59)  nay  be 
c cabined,  using  (2-1)  to  give  the  single  condition 
ur(c+,l)  *ur(c-,l)  +U*Ti,  Hi), 


(2-59) 


and  from  (2-57), 
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Urfc+,4)  s  Wo  fc+,4)  +  Wo  To  F(&) . 


(2-60) 

Finally,  these  conditions  can  he  expressed  in  terms  of  an 
integral  over  the  jet  vortex  distribution  or  potential  difference. 

From  (2-44),  (2-54)  and  (2-60), 


f-* 


-2  Wo  To  FUr)  3 


or,  using  (2-7), 

a*, + wffer 


(2-61) 


(2-62) 


In  the  steady  jet-flap  analysis  of  I  and  II,  the  boundary  condition 
at  the  trailing  edge  can  be  expressed  in  exactly  the  form  of  (2-6l) 
and  (2-62),  less  the  time  dependence.  Spence  has  shown  in  I  that  this 
condition  can  be  satisfied  if  there  is  a  logarithmic  singularity  in 
the  vortex  distribution  at  the  trailing  edge,  and  his  analytic 
solution  of  II  confirms  this.  These  singularities  have  been  discussed 
in  Section  2.3*  If  there  is  no  jet  deflection,  To  *  ®  > 

(2-61)  implies  that  the  vortex  distribution  is  regular  at  the  trailing 


edge- 


2.6  Proof  of  Constancy  of  Circulation 

Since  the  problem  is  linear,  solutions  for  various  boundary 
conditions  may  be  superimposed.  Therefore  a  transient  problem, 
e.g.,  one  described  by  (2-13),  (2-14),  or  (2-15)  and  (2-16),  may  be 


considered  Independently  of  any  other  steady  or  unsteady  configuration 
the  system  of  airfoil  plus  Jet  might  have,  provided  the  assumption 
of  a  Jet  fully  established  In  length  before  Initiation  of  the  transient 
motion  is  met.  The  airfoil  and  Jet  may  then  be  assumed  to  consist  of 
a  flat  plate  with  an  infinitely  long  Jet  tangent  at  the  trailing 
edge,  both  aligned  with  the  free  stream  before  initiation  of  the 
motion.  Then  the  circulation  around  the  following  simply  connected 
contour  is  initially  zero: 


to 


In  the  classical  unsteady  theory  of  airfoils  without  Jets, 
the  flow  about  the  airfoil  is  doubly  connected.  Therefore  the  circulation 
about  a  contour  drawn  around  the  airfoil  is  zero  before  Initiation  of 
the  transient  motion  and  must  remain  zero  for  all  times  after  the 
motion  has  begun.  This  Is  a  direct  conclusion  from  Kelvin's  Theorem, 
cf.  Sears  (19^4).  In  the  Jet- flap  case,  the  simple  connectivity 
of  the  flow  eliminates  this  argument.  Instead,  a  proof  of  the 
constancy  of  circulation  in  the  Jet-flap  case  given  by  Spence  In  III 
will  be  reproduced. 

This  proof  Is  based  on  the  physical  property  that  sufficiently 
far  downstream  at  a  given  instant  of  time  after  initiation  of  the 
motion,  the  Jet  must  return  to  its  Initial  undisturbed  position. 

Moreover,  its  slope  and  curvature  will  also  vanish;  l.e. , 


(2-63) 


IMs  property  may  be  visualized  in  the  following  way.  If  the 
velocity,  Vo  ,  of  the  flow  in  the  jet  is  much  greater  than  the 
velocity  of  the  jet  boundaries,  as  expressed  by  the  downw&sh  velocity, 
UT/V.-fc)  ,  the  jet  will  be  continuous  from  the  trailing  edge 

to  infinity  downstream.  If  the  reverse  inequality  held,  the  jet 
might  cease  to  be  a  continuous  flow  and  break  in  seme  fashion.  For 
the  continuous  jet,  the  assumption  of  inviscid  flow  Implies  that  the 
Jet  is  impenneable;  hence  a  particle  of  air  above  (or  below)  the 
jet  initially  must  always  remain  above  (or  below) .  If,  then, 
the  jet  were  displaced  frem  its  initial  position  at  infinity  downstream, 
an  infinite  amount  of  work  would  have  been  done  in  a  finite  time 
to  move  the  infinite  amount  of  air  above  and  below  the  jet.  The 
condition  that  Vo  »  is  no  restriction, 

since  in  Section  1 . 5  the  assumption  of  Vo  U«  has  been 

made,  and  in  Section  2.1  the  small- perturbation  assumption  of 

U0  »  ur/Y,t)  has  been  made.  Again,  as  in  Sections  1.3 

and  2.1,  it  must  be  mentioned  that  although  unit- step  functions  and 
their  derivatives  may  be  treated,  they  are  considered  as  the  limiting 
cases  of  flows  in  which  Vo  »  W lV-,b} 

To  find  the  circulation  in  the  system  at  any  time,  i.e., 
the  potential  difference  as  ■%-+oo  ,  the 

unsteady  Bernoulli  equation,  (2-12),  is  solved  for  A<Pl%,b) 

This  solution  is 
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A4>/%^)  * 


l 

po, 


+  ,  IstVU 


(2-eo 

vhere  the  motion  is  assumed  to  have  been  Initiated  at  ^  *  o  , 
and  vhere  there  is  no  pressure  difference  across  the  %-  axis 
upstream  of  the  leading  edge.  p (%-  is  an 

arbitrary  function  of  the  integration,  but  since  (2-64)  must  be  valid 
for  all  positive  and  negative  •%■  and  -i.  ,  the  vanishing  of 

for  o  and  1 4-0  require  p7%-U*i r) 

to  be  zero,  identically. 

An  alternate  integration  of  (2-12),  vhich  is  related  to 
(2-64)  by  the  transformation  ^  *  •L  +  ~  is,  since 

F/i-iU)  so  , 


i*V«o 

-  ^ ^(wpuol'y)  itj  ,  . 


(2-65) 


For  H-Uki  ,  with  the  pressure  difference  across 

the  Jet,  (2-5),  (2-64)  becomes 

-  > 

A$/x,^) »  -  i  -  lUL^is)  Jt  J  %>c+u»4. 

TCrV,i 
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Now,  if 

K/x,i)  *  /***  £*  %-0*t  £f  4% 

then 

|a^)|  ^  i  u.^(u#£)  ia%,t)  # 

and  for  a  fixed  4.  ,  by  (2-63), 

j  A$/M)|  *  O  , 

flvwl 

(2-66) 

therefore  the  circulation  in  the  system  is  zero.  Note  that  the 
limit  is  non-uniform  with  respect  to  4  ,  since  for  longer  times 

it  is  necessary  to  go  farther  downstream  to  take  the  limit. 

For  a  Jet- flapped  airfoil  in  unsteady  motion,  then,  it  is 
no  longer  Just  a  case  of  shedding  vorticity  into  the  wake  at  the 
trailing  edge  in  equal,  and  opposite  amounts  to  the  changes  in  circu¬ 
lation  around  the  airfoil,  as  in  the  classical  unsteady-airfoil 
theory.  Bather,  equal  and  opposite  amounts  of  vorticity  must  be 
shed  from  all  points  of  the  airfoil  and  Jet,  in  such  a  way  that  the 
total  circulation  in  the  system  vanishes.  Blmbaum's  concept  of 
bound  and  shed  vortices,  cf.  Cicala  (1941 ),  is  a  good  physical  way 
of  considering  the  problem.  Die  airfoil  and  Jet  are  thought  of  as 

'it, 

being  represented  by  bound  vortex  distributions  whose  strengths  at  a 
point  are  equal  to  $5^  •  ach  of  these  bound  vortices 

sheds,  as  its  strength  varies,  the  amount  of  vorticity  necessary 
to  satisfy  (2-12).  2Ms  shed  vorticity  is  convected  downstream  at 
the  free- stream  speed,  U» 


An  alternate  equation  expressing  the  constancy  of  circulation 
may  he  found  by  considering  (2-1*9)  along  with  (2-7)  and  (2-53)*  giving 

r.lb)  -4*NM)  +A^C4,i). 

(2-67) 

With  (2-66)*  then*  since  the  potential  difference  must  be  continuous 

at  the  trailing  edge* 
to 

-r.w 

(2-6B) 

expresses  the  constancy  of  circulation.  This  is  recognized  as  the 
Wagner  (1925)  Integral  condition  of  classical  unsteady-airfoil  theory. 
The  form  is  identical*  but  If  here  contains  the  effects 

of  the  jet*  in  addition  to  the  vorticity  shed  from  the  airfoil  due  to 
its  motion. 

The  above  proof  and  results  are  readily  extended  to  cover 
divergent  steady- state  oscillations,  O  .  At  i=-oo 

these  started  vith  zero  amplitude*  so  the  jet  must  return  to  its 
initial  position  at  infinity  downstream.  For  pure  oscillations* 

*  care  must  be  exercised*  since  the  oscillations 
began  at  finite  amplitude  at  4*  The  same  difficulty 

arises  In  the  classical  unsteady-airfoil  theory  and  will  be  discussed 
in  Section  5*3* 

2.7  Complete  Equations  for  the  System:  Some  Properties  of  Then 

The  complete  set  of  equations  required  to  state  the  problem 
have  now  been  derived  and  will  be  collected  here  for  convenience.  Since 
the  equations  can  be  written  in  a  variety  of  forms*  depending  on  Which 
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functions  are  taken  as  the  primary  unknowns,  several  alternatives 
will  be  given. 

Hie  dynamic  equation  coupling  the  Jet  to  the  main  stream 
is,  from  (2-8)  to  (2-11),  either 

(2-69) 

or 


D  . _ I  ii*/'/'- 

-  2  V*  <-c3  . 57?- 


(2-70) 

or 

_  _  _L  /•  /■— .  5^2^  .  c*.x.*.ao. 

T>fc  *  *  T  3%* 


or 


gg&ft  *  -iu.’-ccj 


(2-71) 


(2-72) 

Hie  kinematic  equation  giving  the  downwash  required  to 
make  the  Jet  be  a  streamline  at  each  instant  of  time  is,  using  (2-44), 
(2-54)  and  (2-7),  either 


ur/*.,t)  * 


Dy/MO 

TT- 


jmjdi 

t=% 


0  j 

(2-73) 


or 


50 


or 


TT 


JL 

I-* 


where 


(2-74) 


C*-K±ao  , 


(2-7?) 

The  boundary  conditions  at  the  trailing  edge  of  the  airfoil 
are  from  (2-58)  to  (2-62),  either 


yu+fi)  =  Y lc-M 


and 


(2-76) 


3YfcU) 


+  rc  Frt) , 


or 


(2-77) 


UT/c+.i)*  ttf.fc+.-fc)  +  UoT*  , 


(2-78) 


or 


jfom 

y«^c+ 


-2WTo  F^) 


(2-79) 
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The  conservation  of  circulation  in  the  system  is 
expressed,  using  (2-66)  to  (2-6B),  by  either 


(2-80) 


insuring  that 


*  O 


t  (W 


(2-8T) 


is  continuous  at  the  trailing  edge,  or 


J/i'jW-tWl  *  -  P./i)  , 

it 

(2-82) 


c* 

(2-83) 

By  explicit  use  of  the  Itagner  integral  condition  (2-82) 
or  (2-83),  the  downvash  equations  (2-73)  or  (2-74)  can  be  cast  in  a 
form  useful  in  exhibiting  certain  properties  of  the  solution  of  these 
equations.  Using  the  identity 


i  (f-c)  J_ 

*  (%rc)(t~%)  '  1 trC  , 

(2-84) 


equations  (2-73)  and  (2-75)  may  be  written 
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£frj£M& \J1 

i-t  t 


+  -  uTjkv&fy  H^SM H 

H:  c 


and  thrid 

urM).  5g 


4:;*%* 

condition,  (2-82),  the  second 


Fran  (2-53)  >nd  the  Wagner  integral - 

'  “  ' ‘  Integrals  cancel,  leaving  either 


(2-85) 


or 


(2-86) 


where 


- 


i  Hifc-jjf  pyilM&i 

^JLW-oJ  S-r 


(2-8J) 

Finally,  although  decoupled  by  the  inversion  of 
Section  2.3  front  the  equations  now  required  to  solve  the  problem, 
the  vortex  distribution  on  the  airfoil  is,  from  (2-43)  and  (2-51), 


llv.b)  «  fc 


t-ci 


«- * 


0*X**-C 


(2-8&) 


and  can  be  evaluated  once  the  Jet  vortex  distribution  is  found. 

It  is  of  interest  to  notice  that  the  effect  of  the  Jet 
appears  explicitly,  l.e.,  by  appearance  of  C.j  in  the  equations, 
only  in  (2-69)  to  (2-72).  These  are  dynamic  equations  involving  the 
pressure  difference  across  the  Jet,  which  is  proportional  to  Cj 
The  fundamental  downwash  equations  (2-19),  and  the  results  obtained 
therefrom  by  inversion  do  not  explicitly  contain  Cj  ,  although 
the  effect  of  Cy  is  Implicit.  Furthermore,  the  boundary 
conditions  at  the  trailing  edge  are  kinematic,  and  although  a 
dynamic  argument  was  used  to  prove  the  constancy  of  circulation, 


it  may  be  expressed  by  (2-81 )  to  (2-83)  >  which  are  found  from  the 
kinematical  equations. 

In  the  case  of  the  functions  and  , 

single  equations  may  be  found  by  eliminating  W i%,4rt  from  the 

fundamental  pair  of  equations.  For  the  vortex  distribution, 
differentiate  (2-73)  three  times  with  respect  to  %■  and 
eliminate  between  this  and  (2-72),  giving 


yWtttt  _  Mice,  >> 


C+%4-00  . 


(2-89) 

Using  the  transformation  %  *  C+  /%~c)y  in  the  integral, 

the  three  derivatives  may  be  taken,  giving  the  complete  result 


& 


2  >3^  > 


C<-K*#>  f 


(2-90) 


The  alternate  font  of  this  equation,  explicitly  incorporating  the 
Hngner  Integral  condition,  is,  from  (2-72)  and  (2-85), 


tHUK 

t>t*  ~  4i 


L 


^fJ(¥)VM  -  &  fef  2$>A  -  *#$$*#* 

— ■  >1  ®  ^S» 


,  3£ 

t 


-  UoCfr  gjg/Lfe) 

J  >y.#  > 


£  *1  £*100  . 


(2-91) 

Either  of  these  equations  can  then  be  solved,  using  the  conditions 
(2-79)  and  (2-82). 
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Following  the  ease  procedure,  (2-70)  and  (2-74)  yield 
for  the  potential  difference  across  the  Jet 

*  ifvf  v  * :  t  *#*  *} 


2 


(2-92) 


or.  In  the  font  explicitly  incorporating  the  Hhgner  integral  condition. 


tMliM  -  L &S& 


+  f  tijpHt} 

C.  J 


2.  , 


(2-93) 

Either  of  these  nay  be  solved,  using  the  conditions  (2-80)  and 
(2-81)  or  (2-83). 

ttiese  equations  in  the  Jet  vortex  distribution  and  potential 
difference  are  unwieldy.  However,  upon  approximation  for  small 
values  of  Cy  in  Chapter  6,  they  simplify  greatly  and  are  useful 


56 

forma.  Slagle  equation*  in  the  Jet  ordinate  or  dommaab  on  the  Jet 
cannot  be  found,  since  the  Jet  vortex  distribution  appearing  in  the 
integral  term  cannot  be  eliminated  completely. 

Ey  considering  some  special  limits  of  the  equations  of  the 
problem,  Important  asymptotic  properties  of  the  solution  near 
%  »  co  may  be  found. 

First,  if  the  downvash  equation,  (2-73),  is  multiplied 
by  %  ,  and  the  limit  as  %-*■  Co  is  taken  for  finite  fixed 

■fc.  ,  the  result  is 

e  r 

jti*  % 

(a-94) 

Using  (2- 53),  the  tihgner  integral  condition  (2-82)  makes  the  right-hand 
side  of  (2-9*0  vanish.  Therefore 


Jj*  *•  O  , 


or,  asymptotically  in  %  , 


(2-9"' 


US  ~ 


(2-96) 

where  D  >  I  ,  but  as  yet  undetermined.  Equation  (2-72)  for 

the  Jet  Ifcen  implies  ihat 


(2-97) 


or 
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**-3 


(2-96) 


These  latter  conditions  may  also  be  shown  to  follow  from  (2-90) 
by  multiplying  it  by  X*  ,  taking  the  limit  aa  X-*o«>  , 

and  Integrating  by  parts  the  Integrals  containing  ,  thus 

reducing  them  to  &  form  which  vanishes  by  the  Vagner  Integral 


condition. 

Next  treat  the  downwash  equation  in  its  form,  (2-89); 
which  incorporates  the  Vagner  integral  condition.  Multiplying  it 
by  'X-i  and  taking  the  limit  as  00  gives 

>99) 

the  integral  over  ft (9,-tr)  existing  because  of  (2-98).  The 

right-hand  side  of  this  equation  is  in  general  not  zero,  and  will 
be  shown  in  Chapter  3  to  be  related  to  the  lift  on  the  airfoil.  The 
downwash,  then,  behaves  asymptotically  in  X  like 


ur/x.-fc) 


(2-100) 

l.e.,  like  the  downwash  of  a  doublet,  which  follows  since  the  circu¬ 
lation  in  the  system  is  conserved.  Since  m  »  Z  in  (2-96), 

(2-97)  ®nd  (2-98)  give 
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$1M  5 

W  ■* 

(2-101) 

and 

Xlr,€)~  0TS  . 

(2-102) 

lhese  latter  results  any  also  be  shown  directly  from  (2-91),  where 
multipli cation  by  *i?  ,  taldng  the  limit  of  '%•*■ <*>  > 

and  several  Integrations  by  parts  yields 


Jfc*  VjjW)  J*  -jVw  _ 

(2-103) 

Similar  considerations  of  the  potential  difference  equations, 
(2-92)  and  (2-93),  show  that 

c  « 


indicating  that  asymptotically  in  %  > 


(2-104) 


(2-105) 


CHAPTER  3  -  CALCULATE  Of  QP  THE  LEFT  AID  PITCHTHO-MCMEET  COEFFICIENTS 


3-1  Calculation  of  the  Lift  Coefficient 


The  total  lift  on  the  airfoil  ia  the  integral  of  the  pressure 
difference  across  the  airfoil  plus  the  vertical  component  of  the  Jet 
momentum  flux.  This  Jet- momentum  flux  reaction  acts  upon  the  internal 
ducting  of  the  airfoil  and  is  not  an  external  pressure  force.  The  lift, 
positive  upwards,  is  written  in  the  linearized  approximation,  using 
(2-1),  as 


4  j  avfoi) 

*  DX 


(3-D 


Defining  the  usual  non-dimensional  lift  coefficient,  using 
Bernoulli  equation  (2-12),  and  using  the  definition  of  the 
momentum- flux  coefficient  (2-4),  (3-1)  may  be  rewritten 


the  unsteady 
Jet- 


**  • 
(3-2) 


von  Karmen  and  Sears  (1938)  argue  that  since  the  total 
circulation  in  the  system  of  an  airfoil  without  Jet  vanishes,  the 
vortices  occur  In  equal  and  opposite  pairs.  The  lift  may  then  be 
found  by  taking  the  negative  of  the  time  derivative  of  the  impulse 
of  the  flow.  Despite  the  fact  that  part  of  the  lift  in  the  Jet- flap 
case  is  given  by  the  internal  Jet- reaction  force,  the  total  lift  is 
again  given  by  this  result.  To  show  this,  consider  the  following 
integral,  using  (2-8), 
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(3-3) 

where  the  upper  limit  vanishes  by  (2-63).  Therefore,  although  the 
Jet-mcmentuBr flux  reaction  acts  upon  the  internal  ducting,  this  force 
is  represented  by  the  integral  of  the  pressure  difference  across  the 
jet.  Substituting  (3-3)  into  (3-2)  gives 


(3-4) 


The  total  lift  of  the  system,  then,  is  given  by  the  integral  of  the 
pressure  difference  across  both  the  airfoil  and  the  jet.  Integrating 
the  3  ^  term  immediately,  and  taking  the  time 

derivative  outside  the  Integral, 


cm  *  ^  * +  h  *H. 

the  second  term  vanishing  at  both  limits  by  (2-43)  nod  (2-103). 
Integrating  the  remaining  integral  by  parts  and  using  (2-7)  yields 


(3-*.) 


I 


61 


where  the  first  tern  vanishes,  again  by  (2-45)  and  (2-105)  •  Therefore, 

i 

(3-6) 

The  lift  Is  also  given  In  this  case  by  the  negative  of  the  tine 
derivative  of  the  Impulse  of  the  flow. 

To  get  the  lift  coefficient  in  forms  useful  for  calculation,  It 
is  more  convenient  to  treat  (3-2).  Integrating  the  ^ &&(%, £ 
term  by  parts, 

'i  • 

(3-7) 

the  first  term  vanishing  at  both  limits  of  integration  by  (2-22) 
and  (2-45).  Substituting  (3-7)  end  (2-7)  into  (3-2)  and  talcing  the 
time  derivative  outside  the  integral  gives 
c  c 

CJiV-  tc  +  Ct 

(3-8) 

a  general  expression  for  the  lift  coefficient  in  terms  of  the  airfoil 
vortex  distribution  alone. 


To  get  the  lift  coefficient  in  terns  of  the  vortex 
distribution  representing  the  Jet,  substitute  (2-43)  end  (2-31) 
into  (3-8) ,  giving  c» 

CLU)  *  £rc  4  f ~ 

#  (3-9) 

The  order  of  integration  nay  be  Interchanged  In  (3-9) >  using 
the  arguments  following  (2-46),  and  the  resulting  integrals  over  % 
evaluated  by  (A- 51)  and  (A- 53).  Using  these  results  and  the  classical 
lift  terms  discussed  below,  the  lift  coefficient  becomes 

e» 

cjto*  G,«nA« 


(3-10) 

This  is  a  completely  general  result  far  the  unsteady  lift  coefficient. 

The  constancy  of  circulation  in  the  system  is  not  incorporated  in 
this  result. 

The  integrals  over  the  airfoil- quasi- steady  vortex  distribution, 
tfolV-jt)  ,  have  been  identified  as  the  same  terms  which  arise 


in  the  classical  unsteady-airfoil  theory.  These  have  been  investigated 
by  von  Kerman  and  Sears  (1938),  whose  notation  ana  definitions  are  to 
be  used  here.  The  airfoil- quasi- steady  lift  is  that  circulatory  force 
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due  to  the  airfoil  chape  and  notion  alone,  neglecting  any  wake  or  Jet 
contributions  aa  dlacuaaed  In  Section  2.4.  Ita  coefficient  la  denoted 
by  «a4 1*.  u#inj  (2-53). 


CL,W 


(3-H) 


The  airfoil  apparent- mass  lift  Is  the  force  of  non- circulatory 
origin  and  depends  on  the  airfoil  shape  and  motion  only.  Its  coefficient 
Is  denoted  by  Cl,  (-ir)  and  Is  defined  by 

6-.  M *  *  vfC  ^  ' 

0 


(3-12) 

To  get  In  terns  of  the  downwash  on  the  airfoil,  substitute 

(2-51)  into  (3-12),  giving 

c 

Inversion  of  the  order  of  Integration  using  the  arguments  following 
(2-46),  and  evaluation  of  the  resulting  Integrals  over  %  by 
(A- 50)  and  (A- 52)  gives 


(3-13) 


Since  the  potential  difference  must  vanish  at  infinity 
downstream  as  discussed  in  Section  2.6,  (2-66),  this  result,  (3-10), 
can  be  simplified.  Identification  of  certain  of  the  Integral  terns 
in  (3-10)  by  the  Wagner  integral  condition,  (2-82),  and  by  the 
definition  of  the  airfoil- quasi- steady  lift  coefficient,  (3-1 1), 


substitution  of 


3A#*tyj 


substitution  of  (3-3)  for  Cj 
cancellation  of  terms, 

c 


for  Y ly.,4A  ,  and 

give,  upon 


(3-?4) 


Integrating'  the  second  integral  of  (3-14)  by  parts  yields 

\  c 

«0 


the  contributions  at  the  limits  of  integration  vanishing  by  ( 2-22 ) 
and  (2-105)*  The  integral  remaining  cancels  the  last  integral 
in  (3-14),  giving  the  compact  form  for  the  lift  coefficient 

r* 

CJi)=  Cl,U)  . 

(3-15) 

In  this  expression,  all  the  lift  of  a  circulatory  nature  is  included 
in  the  integral. 

The  lift  coefficient  can  be  found  in  terms  of  the  limiting 
behavior  of  the  dovnvash,  vortex  distribution  and  potential  difference 
at  infinity  downstream,  as  promised  in  Section  2.7*  Examination  of 
(3-13)  and  (3-’5)>  in  conjunction  with  (2-99),  (2-103)  and  (2-104) 
gives 


(3-16) 


(3-17) 


and 


CM>WcX,&ki?t^) . 

(3-18) 

For  the  sake  of  completeness,  an  expression  for  the  lift 
coefficient  can  be  found  which  is  of  the  form  of  the  von  K&rm&n  and 
Sears  (1938)  result.  To  show  this,  add  and  subtract  £&•  I4r) 
from  (3-15)  using  the  Wagner  Integral  condition,  (2-82),  take  the 
time  derivative  inside  the  integral  In  (3-1 5) >  and  substitute  (2-10)  ' 
for  inside  the  integral,  giving 

CM  *  ClM  +&..«)  +  fefeJWi* 

+ jjVwr«fc  *&“>  ZgFp*. 


(3-19) 


Integrating  the  Integral  over  by  parts  gives 
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since  the  contribution  at  the  Units  of  integration  vanish  by  (2-22) 
and  (2-102).  Substituting  this  into  (3-19)  sinplifies  it  to 


CLU)  *  Cl,  a)  +ClJ&  + 


(3-20) 


(3-21) 


In  the  classical  case,  Ci.^(b)  is  the  vake- effect  tens  but 

here  includes  the  Jet  effects  as  does  the  last  term.  This  tern 
vanishes  with  Cj  ,  leaving  the  von  Kantian  and  Sears  (1938) 
result. 


3.2  Calculation  of  the  Pitching- Moment  Coefficient 

The  total  pitching  moment,  positive  nose  downward  about 
sene  reference  axis,  is  the  integral  of  the  moments  of  the  individual 
lifting  elements  of  the  airfoil  about  this  axis,  plus  the  moment  of 
the  Jet- momentum- flux  reaction  about  this  axis.  To  determine  this 
direct- Jet- reaction  moment,  consider  the  following  arbitrary  airfoil 
at  some  instant  of  its  motion  . 

D 


•Inc* 


By  the  linearized  approximation, 

i^e±Vi-*)c 


AB  *  Xb 

Geometrically, 

-yrc+,t)+bc. 


The  moment  due  to  this  force  is 

mu)  >  j[«-<uc  *kJ , 

The  total  moment  is,  then, 

MU)*  -Jto-**) h  +■  jEMc  - W*.«  *le]. 

(3 

Defining  the  non- dimenai  onal  pitching-moment  coefficient,  and  using 
the  unsteady  Bernoulli  equation,  (2-12),  and  the  definition  of  the 
momentum- flux  coefficient,  (2-4),  (3-22)  becomes 


+k]  . 


(3- 

As  for  the  lift,  this  equation  may  be  shown  to  reduce  to  a 
form  analogous  to  that  which  von  K&rmin  and  Sears  (1 938)  derived 
considering  the  time  derivative  of  the  moment  of  Impulse.  To  show 
how  this  differs  in  this  case  due  to  the  moment  of  the  internal  Jet- 
reaction  force,  consider  the  following  integral,  using  (2-8), 


which  upon  integration  by  parts  *  m 

-(**!**] 


A  further  Integration  and  use  of  (25-63)  in  evaluating  the  Units  of 
the  Integrals  gives,  therefore, 

W-¥l. 

(3-24) 

Substituting  (3-24)  into  (3-23)  yields 


(3-25) 


Integrating  (3-25)  by  parts,  and  using  (2-7), 

CmU)* 

The  limits  of  the  first  and  third  tezms  vanish  by  (2-45)  and  (2-105)* 


so  the  moment  expression  becomes 

00  Q0 

Cm  (h)  -  -  ft?  +  Cr 

(3-26) 

This  corresponds  to  the  von  K&rman  and  Sears  (1938)  expression  except 
for  the  tens  Cy  b  .  This  term  is  the  moment  of  the  momentum 
flux  at  infinity  downstream  times  the  distance  of  the  reference  axis 
below  the  undisturbed  position  of  the  jet  at  infinity  downstream.  It 
is  not  surprising  that  this  additional  term  appears,  since  the  moment 
due  directly  to  the  Jet  reaction  acts  on  the  internal  ducting  of  the 
airfoil,  not  as  a  pressure  force  on  the  airfoil  surface. 
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To  get  the  moment  coefficient  in  forms  useful  for 
calculation,  t  treatment  to  that  node  In  the  lift-  coefficient 
case  will  be  made.  Integrating  the  term 

in  (3-23)  by  parts  gives 

c 

*  ui^V/^  +zac%-%x]  , 

0 

the  integral  vanishing  at  both  limits  of  integration  by  (2-22)  and  (2-45) . 
Substituting  this  result  and  (2-7)  into  (3-23)  and  taking  the  time 
derivative  outside  the  integral  gives  a  general  expression  for  the 
pitching-moment  coefficient  in  terms  of  the  airfoil  vortex  distribution 
alone, 

c  r 

0,U)*  ^  jipWl-iA) 

v0  ft 

(3-27) 

The  pitching-moment  coefficient  also  can  be  expressed  in 
terms  of  the  vortex  distribution  representing  the  Jet.  Substitute 
(2-43)  and  (2-51)  into  (3-27),  giving 
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Cm  U)  =  {j^v  ^|/lrA«)  fcfctl  k 

+tte- •  iT^?  ^ 

V  c 


(3-26) 

The  order  of  integration  again  nay  be  interchanged  by  the  argument 
following  (2-46)  and  the  resulting  integrals  over  %  evaluated  by 
(A- 51 ),  (A- 53)  and  (A- 55)  •  Using  this  and  the  classical  terms 
discussed  below,  the  pitching-moment  coefficient  reduces  to 

-  i it? 

Jc 

+ ajfr*-  *  4-*>  >JSP  ■ - ' ^  A . 

c 

(3-29) 

This  is  a  general  result,  the  constancy  of  circulation  not  having 
been  used. 


The  Integrals  over  the  airfoil-quasi- steady  vortex 
distribution  in  (3-29)  have  been  identified  In  terms  of  the  airfoil- 
quasi-  steady  and  apparent- mass  lift  coefficients  defined  in  (3-1 0  to 


I 
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(3-13)>  and  the  analogous  pitching- moment  coefficients  discussed 
by  von  Karmln  and  Sears  ( 1 938 ) .  The  airfoil-  quasi- steady  pitching  moment 
is  that  moment  of  circulatory  origin  due  to  the  airfoil  shape  and  motion 
alone,  neglecting  wake  and  Jet  contributions.  Denoted  by  C*c 
its  coefficient  is 

c 

C/nt  l &  *  J*  . 

O 

(3-30) 

By  the  same  argument  leading  from  (3-12)  to  (3-13)>  this  can  be  written, 
using  (A- 50)  and  (A- 52),  as 

C*eM)*  -^S&C'tf'vrtuUi  . 

•  (3-3i) 

By  convention,  this  is  a  moment  about  midchord  (a*V)  •  n* 

airfoil  apparent-mass  pitching  moment,  also  about  midchord,  is  of  non- 
clrculatory  origin  and  depends  on  the  airfoil  shape  and  motion  only. 

Denoted  by  lb)  ,  its  coefficient  is  defined  by 

c 

Cm, Ur)  * 

(3-32) 


By  the  argument  leading  from  (3-12)  to  (3— 1 3) >  this  also  can  be 
written,  using  (A-50),  (A- 52),  and  (A- 54),  as 


t 


(3-23) 


By  use  of  the  condition  that  the  potential  difference  at 
infinity  downstream  vanishes,  (2-66),  this  result,  (3-29),  can  be 
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simplified.  Identification  of  certain  of  the  integral  terns  in 
(3-29)  using  the  Hsgner  integral  condition,  (2-82),  and  the  definition 
of  the  airfoil- qpasl- steady  lift  coefficient,  (3-11),  substitution 
0f  f0r  ft !%,{;)  by  (2-7)  where  needed,  and  substitution 

of  (3-24)  for  CapK)^^'  -  give,  upon  cancellation 

of  terms, 

Gn(  4)  *  c*.ih)  +Cm,w  +  i  it-2A)Cita) 


c  C 


&c)  +  Cj  b  t 


(3-34) 


Integrating  the  second  Integral  in  this  expression  by  parts  yields 

it ,  ^pwMr«wK)Ae»)]c 

fir  ^ 

-  (fa 


,  el 
*  «??  35 


OO 


the  contributions  at  the  limits  of  integration  vanishing  by  (2-22) 
and  (2-105) •  The  remaining  integral  cancels  the  last  Integral  in 
(3-34),  which,  upon  factoring  of  the  first  integral  beccaes 
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CM*C»,M  +  Cm 


I 

~5*> 


fr/f>WS  +CtI. 


(3-35) 

To  calculate  the  retching- moment  coefficient  by  (3-35)  ,  only  one 
more  Integral  over  the  Jet  vortex  distribution  need  be  calculated 
besides  the  Integral  required  for  the  lift  coefficient  In  (3- 1 5) * 

An  expression  for  the  pitching-moment  coefficient  of  the  form  of 
the  von  Karmfin  and  Sears  (1938)  expression  is  found  from  (3-35)  by 
adding  and  subtracting  -za)  £/..  Cir)  from  (3-35)  using  the 
Vagner  integral  condition,  (2-82),  taking  the  time  derivative  Inside 
the  Integral  and  substituting  (2-10)  for 
This  gives 


CM*  CmoOt)  +  +iO-u.)a0fr) 


fri-cf tmus  fG b. 


(3-36) 


terns  by  parts  gives 


7* 


Integrating  the 


frW* mji'/i  +  kfc/t-J'igpdt  - 

't  * 

^  oo 


j.  [m&A t  - 

HfcJ  Wi-c)7*-  ~ 


a^L^it-ctmui 
**  1 


to  * 

-  l!z™±^£jU  Wf.iUi  -  > 


where  the  contributions  at  the  limits  of  integration  vanish  by  (2-22) 
and  (2-102).  Substituting  this  into  (3-36)  and  cancelling  terms  gives 

Cm  1^)*  +  Cm , Ur)  014;)  +  Cl, /^)J  +  ^  Ci%U?S 

+  itJt)  C^lt-ct  +  £rb# 

(3-37) 

This  form  corresponds  to  the  result  of  von  Earaan  and  Seams  (1938) • 

In  addition  to  their  terms,  there  are  two  additional  integrals, 


i 
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proportional  to  Cj  >  and  also  the  C-g  lo  tern  discussed  ahovo. 

If  Cy  SO  ,  the  classical  results  are  recovered,  and  the 
moment  contains  only  the  same  terns  which  were  evaluated  to  find 
the  lift  coefficient,  (3-20). 

There  are  no  expressions  for  the  pitching-moment  coefficient 
in  terms  of  the  behavior  of  the  Jet  at  Infinity  downstream  as  were 
found  for  the  lift  coefficient. 

Die  pair  of  equations  (3-15)  and  (3-35)  are  the  simplest 
to  use  for  the  calculation  of  the  lift  and  pitching-moment  coefficients, 
since  only  two  integrals  over  the  jet  vortex  distribution  need  be 
calculated.  For  convenience  they  are  written  here. 


CM)*  C..H)  -  QL 


if*"' 


cfirtividi , 


(3-38) 


and 


(3-39) 


CHAPTER  4  -  EQHAHCHB  FOR  PARTICULAR  PROBLEMS 


4,1  Jet- Deflection  Problem 

The  fundamental  Jet-flap  problem,  both  In  the  steady  and. 
unsteady  cases  is  that  of  "Jet  deflection,"  i.e.,  the  emergence 
of  the  Jet  from  the  trailing  edge  at  a  prescribed  angle  with 
respect  to  the  slope  of  the  airfoil  ordinate  there.  This  problem 
has  been  called  by  Malavard  (1957)  the  "singular-blowing"  problem  in 
recognition  of  the  logarithmic  singularity  in  the  vortex  distribution 
at  the  trailing  edge  required  to  satisfy  (2-61 ) .  Practically,  Jet 
deflection  might  be  achieved  by  a  very  small  flap  at  the  trailing 
edge,  or  by  internal  ducting,  and  is  assumed  to  be  adequately 
represented  by  the  present  model.  She  importance  of  the  problem 
resides  in  its  being  unique  to  the  Jet- flap,  for  if  Cx  *  O  , 
specification  of  the  Jet  deflection  angle  can  no  longer  be  made  and 
the  problem  is  trivial. 

As  a  model  for  the  Jet- deflection  case,  a  flat-plate  airfoil 
aligned  with  the  free- stream  direction  is  chosen,  the  Jet  having 
the  time- dependent  angle  of  deflection,  (2-13), 


T(i)  -  To 


(4-1) 

For  transient  inputs  this  problem  is  of  considerable 
practical  importance,  since  the  jet* flap  has  been  proposed  as  a  control 

*  i 

mechanism  and  the  response  of  the  airfoil  to  such  Inputs  is  basic  to 
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the  full  understanding  of  such  usage.  Oscillating  inputs  are  also 
important,  since  use  of  the  Jet  flap  for  cyclic  pitch  control  of 
helicopter  rotor  blades  has  been  suggested  by  Dorand  (1959) • 

Since  the  airfoil  is  aligned  with  the  free  stream. 


y/*,o  -  o , 

and  hence  from  (2- 1 ) , 


(4-2) 


0*.% . 


(4-3) 

All  the  airfoil- quasi- steady  and  airfoil  apparent-mass  terns  are 
identically  zero. 

The  basic  equations  sure,  in  terms  of  ISl't.i)  and 

ux/y.,i)  -  the  others  sure  omitted  here  and  in  the  remainder 

of  this  chapter  for  simplicity  -  from  (2-69)  to  (2-87), 


w  " 


y*C  ^  |U) 

i-  4  »  a 


and  either 


(4-4) 


<t 


(4-5) 
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with  either 


(4-6) 


urfc+.-t)  *  UoX#  FU), 


(4-7) 


2Uot,  F/4;. 


(4-8) 


00  '4 

j(£j  •  o  . 


(4-9) 


dace  this  set  has  been  solved,  the  other  properties  of 
interest,  e.g.,  the  lift  and  pitching-moment  coefficients,  may  be 
found  from  the  appropriate  results  of  earlier  sections. 

4.2  Problem  of  Airfoil  in  Plunging  Motion 

The  problem  of  a  flat-plate  airfoil  moving  in  a  purely 
plunging  motion  vith  a  Jet  emerging  tangential ly  from  the  trailing 
edge  may  be  described  by 
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V%,t)  *  FH)  ,  Q*-%cc.  , 


(4-10) 


where  h0  is  the  non- dimerisi oual  amplitude  of  the  plunging 
motion,  and  is  a  small  parameter  consistent  with  linearization  of 
the  problem.  This  problem  does  not  have  much  application  in  transient 
motion,  but  in  oscillating  motion  its  results  are  important  for 
determining  the  coupled  binary  bending- torsi  on  and  tertiary  bending- 
torsion,  control- surface- rotation  flutter  stability  criteria. 


From  (2-1)  and  (4-10),  the  dovnwash  on  the  airfoil  is 


urit.b)  -  ho C  F'**)  #  0* TL<-C . 


(4-11) 


With  (4-11),  the  important  airfoil- quasi- steady  terms  in  the  basic 
equations  are  from  (2-51),  (2-53),  (2-54)  and  (2-8 7),  using  (A-9), 
(A- 38),  and  (A-39) 


&/?.*>  *2 kcCffifFU)  ,  o^c, 

(4-12) 


uTo/r,t)  *  u[>  *  Pt)4]  , 


(4-13) 


Pc  Li)  *  nhc?  f'Li), 


(4-14) 


and 

«r,h*\  •  U[i-  frf- C*.  , 


(4-15) 
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The  basic  equations  are,  then,  from  (2-69)  to  (2-87), 

I Tth&  tic  , 

TN+*  *  "  T  , 


L<-  %*-<*> i 


(4-16) 


and  either 


(4-18) 


with  either 


*  Jio  c  fWj 


(4-19) 


or 


A* 


milht  -  o 


(4-20) 


and 


\l0m i*  *  -vU'rtt). 

(4-21) 


Upon  solution,  the  other  properties  nay  he  found.  The 
airfoil- quasi- steady  and  airfoil  apparent-nass  coefficients  are, 
given  by  (3-n),  (S-^),  (3-31)  «nd  (3-33),  using  (A-18)  and  (A*21), 
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and 


CP'(l) 

~ur  , 


Ci.,14:)*  Sj* 


Cx  F'U) 

~uT", 


C*M 


.ttL  cfV£) 
*  u«  > 


U(f)>  o . 


(4-22) 


(4-23) 


(4-24) 


(4-25) 


4.3  Problem  of  Airfoil  in  Pitching  Motion 

A  flat-plate  airfoil  moving  In  a  pitching  motion  about  an 
axis  a  distance  £C  behind  the  leading  edge  may  be  described 
by 


Y/M)=  «*(x-ec)fH) ,  o *  , 

(4-26) 

where  0t»  is  the  amplitude  ocT  the  Incidence  angle  and  Is 
nail  consistent  vlth  the  linearisation.  The  jet  emerges  from  the  air¬ 
foil  tangentially  at  the  trailing  edge,  like  the  plunging  case,  the 
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results  of  this  problem  for  oscillatory-  notions  are  Important  in 
the  binary  and  tertiary  flutter  analyses. 

The  downwasb  on  the  airfoil  is,  by  (2-1 )  and  (4-26), 


urfc.On  Ued.F(4)  4  oU(x-ec)F'te),  . 


(4-27) 

The  important  airfoil- quasi-  steady  quantities  in  the  basic  equations 
are,  from  (2-51),  (2-53),  (2-54)  and  (2-87),  using  (A-9),  (A- 11), 

(A-38),  (A-39),  (A-42)  and  (A-43), 

XJtM •  aud.(^f{ Fli)  4  [* 4rt-eW)  £^] , 

(4-28) 


UJoltM*  Uod.F(±)[l-/*£f]  +  do  fwfex-ec) , 


C  co  t 


(4-29) 


Po^»  TrUodoCFW  +1Td.cV4-«)F/^i, 


and 


(4-30) 


witti  ’  m.  F«)[|  -  (2*f-  ,7nra*.] 
4-tl.FW)f/t-«d-r*4«-<)c](l 'if  - 


(*-31) 
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B»e  basic  equations  are,  again  from  (2-69)  to  (2-87), 


im#  .  _  Ufe  rT 


A  ’  3-itJ  ,  CX%x«oJ 


and  either 


«D 


♦  «(.&.-«)-  &jj  f'ftl , 


or 


(4-32) 


(4-33) 


«*.«  •  • -if  (M Wntf-  £**!«*» 

+*[fv€d-&»  ^]fw, 


C4V14« 


J  (4-34) 


with  either 


or 


urfc+.tf  =  Uod.  Ftt)  +«(.McFW, 


t t^ntm  •  ° . 


(4-35) 


(4-36) 


and 
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&(lMi  *  -7T  O.eU  FUr)  -Ttfoc'Cj-t)  F'll). 


(4-37) 


The  airfoil- quasi- steady  and  airfoil  apparent-mass  lift 
and  pitching-moment  coefficients  are,  by  (3-11),  (3-13),  (3-31) 
and  (3-33),  using  (A-l8),  (A-21)  and  (A-23), 


Cl.  Li)  «  in  M.  F/«  + w ^  -d)  ^  > 


(4—38) 


Ct,  H=>  -  +  Sgs^-e) 

(4-39) 

Fft)  -  3g»  (i-e)  £££> 

(4-40) 

and 

(4-41) 


4.4  Problem  of  Blown  Flap  in  Unsteady  Motion 

Another  practical  jet- flap  system  which  has  been  considered 
is  the  so-called  "blown  flap,"  or  "Jet- augmented  flap."  In  this 
configuration  the  airfoil  has  an  ordinary  trailing- edge  flap,  the  Jet 
emerging  at  the  hinge  point,  following  the  flap,  and  leaving  tangentially 
at  the  trailing  edge.  This  has  been  represented  in  the  thin- airfoil 
model  as  a  bent  flat  plate  with  a  jet  emerging  at  the  trailing  edge. 
Spence  (1958)  has  discussed  the  physical  considerations  of  such  a 
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model.  Here  the  problem  concern*  the  time- dependent  motion  of 
the  flap  about  its  hinge  point,  which  is  located  a  distance  EC 
ahead  of  the  trailing  edge,  l.e.. 


y/t,^  *  r  o  ,  o*%*u-*\c 

,  U-*K 

(4-4z) 

Both  the  transient  and  oscillatory  responses  are  of  importance  for 
this  problem,  the  fomer  because  of  the  control  use  of  the  flap,  the 
latter  because  of  its  importance  in  the  tertiary  bending-torsion, 
control- surface- rotation  flutter  characteristics . 

Defining  an  angle  by 

£  *  , 

(4-43) 

the  downwash  on  the  airfoil  is,  from  (2-1)  and  (4-42), 

ur/tvfe'l  O  4  o*% 

|  U.p.  R4)  lx- cut  J.J  F7«  ,  CU4%±  *  % <c. 


(4-44) 

The  ai rf oil- quasi- steady  quantities  in  the  basic  equations  are, 

from  (2-51),  (2-53),  (2-54)  and  (2-87),  using  (A- 36),  (A- 37),  (A-40),  (A-4l), 

(A-8),  and  (A-10), 


"  -*]Jh 


(4-45) 


^  W)  *  ^[iW [g^f &*£*]  -  ^&)'*j  W 

■*■  * *f  f){*r[j^^  *]  ♦  [l  \<*  1  *» v  j  pw; 

(4-46) 


fo^)  -  ^  O  FTt)  +  ■*'W*»Ks/a-o*>s)Jp#/^), 


(4-47) 

arid 


fW^fJ  -  [\«-  l^v^-'O]  ttff 


C^%<-4 O  . 


(4-48) 
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The  basic  equations  are,  from  (2-69)  to  (2-87), 


C*%<- «o  f 


(4-49) 

and  either 

to 

+ r<) 
c 

+^-£3./lL-cc«>t'$.) -  tv* *  £6*f«*^“SfaO]  F'l4rl 


or 


vrirM  * 


(4-50) 


-  f  *  -S?[»r*- <**£]  i*-  [t~b(£T] 


-|>*-  -  f  m , 

C*.%±*o  f 

(*-50 


with  either 
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txr(c+,t\  *  U*p.Rt) 

(4-52) 

or 

feAOtf&W-o. 

(4-53) 

and 

4  SW.\,/x-ce«\.)3  F'^)  ,  (4-54) 

Airfoil- quasi- steady  and  airfoil  apparent-mass  lift 
and  pitching-manent  coefficients  are,  from  (3-1 1),  (3— 1 3) >  (3-31 )> 
and  (3-33) >  using  (A-17),  (A- 20)  and  (A-22), 

CLo  iVl '  ip.  Fli)  f ^  tKtiW  +*■  \fe-w-iO]  CJ^\ 

(4-55) 

Cl,#)* 

(4-56) 

C*ft)  -%  »v»0  W)  -  »£2^k] 

(4-57) 

and 

6 i.m*  -fe 

(*-58) 


If  the  limit 


,  or  V® 


, 


is  taken 


in  this  problem,  the  equations  reduce  to  the  pure  Jet- deflection 
equations  of  Section  4.1,  provided  care  is  taken  at  the  trailing  edge 
to  preserve  the  Jet- deflection  angle  in  the  limit  of  vanishing  flap 
chord.  Also,  in  the  limit  E  "*>  I  ,  or  K.  — •  7T  ,  the  equations 

reduce  to  those  of  an  airfoil  in  pitching  motion  about  its  leading 
edge,  i.e.,  the  limit  of  €.-*0  in  the  equations  of  Section  4.3* 


4.5  Problem  of  Airfoil  Bitering  Sharp- Edged  Oust 

As  seen  from  equation  (2-15),  there  are  many  possibilities 
of  gust  configurations  which  might  be  considered.  The  fundamental 
transient  gust  problem  which  has  been  considered  in  the  literature  is 
that  of  a  sharp-edged  gust  of  constant  upwash  amplitude  moving 
over  a  flat  plate  airfoil  with  relative  speed  ^°/\  » 

\H(xM  =  -  ^j). 

(4-59) 

For  A*  I  ,  this  is  the  well-known  Kfissner  (1936)  problem,  cf. 

von  K&xman  and  Sears  ( 1 938 ) .  Miles  (1 956)  generalized  it  for  all 
positive  and  negative  values  of  X  ,  positive  values 
corresponding  to  gusts  proceeding  over  the  airfoil  from  the  leading 
edge  to  the  trailing  edge,  and  negative  values  to  gusts  overtaking 
the  airfoil  from  behind.  The  corresponding  oscillatory  problem 
is  defined  by 

v  _  i  u)(4 


(4-60) 

This  is  the  Sears  (1941)  problem  of  a  sinusoidal  gust  of  constant 

amplitude  if  I  ,  and  was  extended  by  Kemp  (1952)  for  a  general, 
complex  X 
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For  the  Jet- flapped  airfoil,  the  presence  of  the  Jet 
extending  to  infinity  downstream,  along  tilth  the  assumption.  Implicit 
in  Section  2.3,  that  the  disturbances  of  the  Jet  must  die  out  at 
Infinity  downstream  indicate  that  the  transient,  sharp- edged  gust 
moving  over  a  flat  plate  airfoil  with  tangential  Jet  from  the 
leading  edge  towards  infinity  downstream,  \>0  ,  is  a 

reasonable  problem  to  treat. 

The  dcwnwash  distribution  of  the  gust,  (^59 ),  implies  that 
the  basic  downwash  equations  for  the  system,  (2-19),  rewritten  to 
express  the  additional  downwash  which  must  be  induced  by  the  vortex 
distribution  to  cancel  (4-59)  must  be 

W(%M)  tw  1/4-  iJ)  -  “  jir  ^  ,  0 , 

0 

(4-61) 

where  W/itt4r)  is  still  given  by  (2-1).  The  Inversion  of 
Section  2.3  again  holds  for  this  equation,  so  the  airfoil- quasi- steady 
expressions  may  again  be  calculated  in  the  usual  fashion,  replacing 
t by  Wd/4-\S)  in  (2-51),  (2-53), 

(2-54)  and  (2-87)  and  using  (A-7),  (A-34)  and  (A-35),  giving 


(4-62) 
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^  exU 

4  j 


(4-63) 

C  iwcfacfigf]  -  [Mlk-M)fj 

Cw4 

[^  t rwc  00  J 

(4-64) 


Wlt.-L)  sA 


v. 


Cd.yu£.t»  , 

*  [l "  '  ^kd.tot 


(4-65) 

The  fundamental  equations  are,  then,  using  the  inversion 
results  of  Section  2.3  with  (4-59),  and  (2-69)  to  (2-87) » 
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at*  TTCy  >  c*%**t 


rw 


(4-66) 


and  either 


-  id^^kr  w  ♦  vwg&uni 

“Mf  •“IWI ,  •-**« 


xw 

Tf 


t4«p( 


-  i(^f  |(AT i^f  -  w[i-(^fJi/«-^) 


(4-67) 


or 


urft.-fc)* 


± 

zrr 


+*<«*!* 

c 


C**-4CO, 


-iferw.  »»o-«r- 


-» ta- 


(4-6B) 
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with  either 


urfc+,fi*o , 


(4-69) 


(4-70) 


WMii  * 


-irBci^-^)  .  %**■■*">  . 


<4-71) 


The  airfoil-  quasi- steady  and  airfoil  apparent-mass  lift  and 
pitching-moment  coefficients  are,  calculated  in  the  same  fashion  as 
(4-62)  to  (4-65),  from  (3-11),  (3-13),  (3-31)  and  (3-33),  using  (A-l6) 


and  (A- 19) 


Cl.W« 


Cl.M' 


0  )  k £  f 


(4-72) 


(4-73) 


C*MV 


U.  FCX 


O, 


(^75) 


CHAPTER  5  -  LIKITI1I0  THEORIES  OP  THE  UX8TEADY  JEP-JLAP  THEORY 


5.1  Reduction  of  the  Equations  to  the  Classical  Has teady- Airfoil  Theory 

The  classical  theory  of  the  unsteady  motion  of  thin  airfoils 
without  Jets  is  contained  in  the  shove  formulation.  To  recover  the 
fundamental  equations  of  that  theory  it  is  sufficient  to  set  Cy  !  O 
formally,  provided  certain  remarks  about  the  equations  are  made.  The 
vortex  distribution  behind  the  airfoil, 

and  the  down  wash  behind  the  airfoil,  LXTlt.i')  ,  C^%£.«o> 
are  convenient  functions  to  treat. 

As  discussed  in  8ection  2.7,  the  Jet  effects  appear 
explicitly  only  in  the  dynamic  Jet-interaction  equations  (2-69) 
to  (2-72).  All  the  equations  derived  from  the  inversion  of  the 
kinematic  dcwnwash  equations,  (2-19),  do  not  change  if  CySO  , 
provided  the  l/square- root  singularity  at  the  trailing  edge  is 
understood  to  be  excluded  by  the  Kutta  condition.  Therefore  (2-73) 
and  (2-88)  remain  the  same.  In  the  absence  of  the  Jet,  it  is  no 
longer  possible  to  prescribe  the  slope  of  the  Jet  deflection  at  the 
trailing  edge.  Hence,  at  the  trailing  edge  the  downwash  is  always 
continuous,  (2-57)*  The  circulation  in  the  system  is  constant  by 
the  argument  of  Section  2.6  using  Kelvin's  Theorem,  so  the  Wagner 
Integral  condition,  (2-68),  which  follows  from  a  kinematic  equation, 
remains  the  same.  Therefore  the  only  equation  which  formally  chaiges 
is  the  Jet- interaction  equation,  (2-71),  say.  The  pressure 
difference  across  the  %  -  axis  behind  the  airfoil  is  proportional 
to  Cy  ,  so  in  the  absence  of  the  Jet, 
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(5-1) 


Ibis  equation  is  a  statement  of  Helmholtz'  Theorem  for  the  conservation 
of  vorticity  of  a  fluid  particle  convecting  with  the  stream,  cf .  Sears 
095*0*  Equation  (2-69)  becomes,  when  Cj  SO  ,  using  (2-7), 

(2-45)  and  (2-49), 


(5-2) 


and,  in  particular  for  , 


w,  ,x  j. 


(5-3) 

A  clear  picture  of  the  physical  nature  of  the  problem  can 
be  seen  from  (5-0  and  (5-3)*  As  the  circulation  around  the  airfoil 
changes  with  time,  vorticity  is  shed  off  the  trailing  edge  in  equal 
and  opposite  amounts  according  to  (5-3)*  This  vorticity  then  convects 
with  the  free  stream,  its  strength  remaining  constant  by  (5-0* 

Solution  of  the  problem  has  become  considerably  simplified 
by  the  disappearance  of  the  downwash  from  (5-1 ),  because  .it  can 
be  solved  immediately  for  the  vortex  distribution,  without  simultaneously 
solving  for  the  downwash.  According  to  Chapter  II  of  Webster  ( 1 95 5 ) > 
the  solution  of  the  first-order,  homogeneous,  partial  differential 
equation,  (5-1);  Is  that  the  vortex  distribution  must  be  a  function 

l  T-t 

of  the  characteristic  variable  *t  -  ■=^* 


; 


namely, 
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Ift, t)« 


(5-4) 

The  functional  form  of  ^71-  -3^)  is  then  determined  by 

considering  the  Wagner  Integral  condition,  (2-82),  as  an  integral 
equation  for  If  , 

ttk-tgjiiz-nti). 

(5-5) 

This  is  the  problem  considered  by  von  Harman  and  Sears  (1938),  and 
many  others,  although  the  equations  were  derived  there  in  a  more 
physically  direct  manner. 

Once  the  vortex  distribution  in  the  wake  behind  the 


airfoil  is  known,  the  dovnwash  behind  the  airfoil,  (2-73),  the 
vortex  distribution  on  the  airfoil,  (2-88),  and  the  other  properties 
of  the  airfoil  and  wake  may  be  found.  The  lift  and  pitching  moment 
coefficients  may  be  calculated  from  (3-20)  and  (3-37)  as  discussed 
in  Chapter  3- 

A  sharp  contrast  may  be  drawn  now  between  the  classical 
problem  and  the  Jet-flap  problem.  Classically  the  vortex  strength 
in  the  wake  behind  the  airfoil  depends,  cf.  (5-1)  and  (5-3);  on  the 
history  of  the  circulation  around  the  airfoil,  or  as  expressed 
through  the  Vagner  integral  equation  (5-5);  on  the  airfoil- quasi¬ 
steady  circulation.  Thus  is  the  important  input 

to  the  problem,  the  actual  distribution  of  down  wash  on  the  airfoil 
having  been  Integrated.  Many  different  chordwise  downwash  distributions 
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would  give  the  same  vortex  distribution  la  the  woke  if  t hey  led 
to  R(«  of  the  sane  tine  dependence.  In  the  jet- flap  problem 
the  vortex  distribution  representing  the  Jet  depends  not  only  on  the 
history  of  the  circulation  around  the  airfoil,  but  also  on  the  details 
of  the  downwash  distribution  at  every  point  of  the  Jet,  through  the 
dynamic  coupling  expressed  by  (2-69)  to  (2-72) . 

5.2  Properties  of  the  Classical  Transient  Solutions 

The  two  types  of  transient  problems  are  those  in  which 
the  motion  of  the  airfoil  is  initiated  at  isO  ,  in  accordance 
with  (2-1 4)  and  (2-1 6),  and  those  in  which  the  leading  edge  of  the 
airfoil  enters  a  gust  at  b  M0  ,  in  accordance  with  (2-15)  and 
(2-16).  In  these  linearized  problems  the  shed  vortex  wake  behind 
the  airfoil  grows  in  length  with  rate  (Jo  after  initiation  of 
motion,  so,  at  any  time,  b  ,  it  extends  downstream  to  »b 

Hie  fundamental  equation  to  be  solved,  then,  is  the  Wagner  integral 

equation,  (5-5);  written  as 

C4Uoi 

fcfw*- • 

c  (5-6) 

For  airfoil  motion,  the  fundamental  problem  is  the  so-called 
Wagner  (1925)  problem,  where  the  airfoil- quasi- steady  circulation  has 
step  function  time  dependence, 


r./ti*  u.t  ±ft> . 


(5-7) 

Once  the  solution  of  this  problem  has  been  found,  solutions  for  other 
transient  R  ft)  may  be  found  by  Duhamel  superposition.  The 


application  of  tbe  Vaguer  problem  to  actual  airfoil  motions  of  the 
types  considered  In  Sections  4.2  to  4.4  has  been  misunderstood  at 
tines  in  the  literature.  Bsference  to  the  particular  problem*  of 
the  above  Sections  will  clarify  this  assertion.  Tor  an  airfoil  in 
plunging  motion  with  F(t)*  ilt)  ,  l.e.,  plunging  at  a 
constant  speed,  (4-l4)  gives 


Glib 

l  f 


(5-8) 


a  direct  application  of  the  Wagner  problem.  On  the  other  hand, 
for  an  airfoil  in  pitching  motion  with  Ftt) *  3.^4)  }  i.e. ,  the 
airfoil  being  snapped  up  instantaneously  to  an  angle  of  incidence, 
(4-30)  gives 


r.M.  mUU»ci{4)  -*)  £&), 

(5-9) 

■where  the  Dirac  delta  function,  S  ,  is 

««•  , 

and  has  the  properties  that 


ill)--o  ,  l  #o' 

but  such  that  ao 

JSIAdfc  «  i  f 

— oo 


(5-10) 


It  is  important  to  recognize  that  the  term  must  not  be 

neglected,  as  It  Is,  for  Instance,  by  Boblnson  and  Lauxmann  (1958) 
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on  p.  505.  Physically  this  term  represents  the  limit  of  the  motion, 
which,  although  occurring  in  a  vanishingly  will  time  Interval,  occurs 
at  such  a  rapid  rate  that  the  result  is  a  finite  rotation  through 
the  angle  do  Considering  this  term  as  such  a  limit,  it 

is  clear  that  it  will  have  an  important  effect,  particularly  for 
small  times.  Therefore,  only  for  an  airfoil  pivoted  about  the 
three-quarter  chord  point,  6  *  Vf  ,  is  the  Wagner  problem 
directly  applicable  to  the  airfoil  snapped  to  angle  of  incidence, 
likewise,  the  trailing- edge  flap  snapped  to  an  angle  of  deflection 
with  F^)  *  d-lty  gives  an  airfoil-  quasi-  steady  circulation, 

(4-47),  'Which  contains  both  IF)  and  £  ('t)  , 

(5-11) 

The  fundamental  gust  problem,  that  of  a  sharp- edged 
gust  of  constant  amplitude  fl  moving  with  a  speed 
relative  to  the  airfoil,  X >°  ,  is  formulated  in  Section  4-5« 

The  airfoil- quasi- steady  circulation  is,  from  (4-64), 

-  tBcfwggy ]  -  H 

r.w  •  • 

(5-12) 

Sharp-edged  gusts  overtaking  the  airfoil  from  behind,  O  , 

can  also  be  treated,  cf.  Miles  ( 1 956) ,  but  cannot  be  extended 
readily  to  the  Jet- flap  case  and  will  not  be  treated^iere . 
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A  method  of  solution,  valid  for  all  transient  |#( 4)  , 
vill  now  be  outlined  and  certain  results  noted.  Upon  the 
transformation 


1  U>  » 


(5-13) 


(5-6)  becomes 


(5-14) 


In  a  similar  fashion  the  lift  coefficient  due  to  the  wake,  ClxU^  , 
is,  from  (3-21), 


-  m  .  X  f, 


(5-15) 


f 


and  finally,  from  (2-67),  (2-68),  (2-7)  and  (5-13), 

*  '  '  ' 

(5-16) 

Sears  (1940)  observed  that  the  integrals  in  (5-14)  to  (5-16)  are 
of  the  convolution  type,  and  used  Laplace  transforms  as  a  convenient 
technique  for  treating  them.  Defining  the  Laplace  transform  of  a 
function  by 


00 


(5-17) 


A  method  of  solution,  valid  for  all  transient  P*/ 1)  , 
will  now  be  outlined  and  certain  results  noted.  Upon  the 


transformation 


i  _  _  T 

t  U,  ~  **  * 


(5-13) 


(5-6)  becomes 


(5-14) 

In  a  similar  fashion  the  lift  coefficient  due  to  the  wake,  , 


is,  from  (3-21), 


t 

5  ttj , 


(5-15) 


and  finally,  from  (2-67),  (2-68),  (2-7)  and  (5-13), 


JVlTldt  «  -jj; 


(5-16) 


Sears  (1940)  observed  that  the  Integrals  in  (5-14)  to  (5-16)  are 


of  the  convolution  type,  and  used  Laplace  transforms  as  a  convenient 
technique  for  treating  them.  Defining  the  Laplace  transform  of  a 


function  by 


00 


(5-17) 


Sears  (19**0)  showed  that 


102 


2 


2. Up  fy?) 

cpP^i>®  , 


and 


(5-18) 


(5-19) 


-  tc 


)(«(IEi)  r*(  p) 

jMuy^isSi 


j 


(5-20) 

where  He/§f0)  and  are  the  Modified  Bessel 

Functions  of  the  Second  Kind.  From  the  definition  of  the 
airfoil- quasi- steady  lift  coefficient,  (3-11),  and  using  the 
airfoil  apparent-mass  lift  coefficient  appropriate  to  the  particular 
motion  considered,  (3-13);  the  lift  coefficient,  (3-20),  is  upon 
transformation, 


CUf)  --  Cl, tp )  4 


Dw&yw&Vj , 


For  the  Whgner  problem,  (5-7); 


(5-2T) 


(5-22) 
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and  the  well- known  results  emerge.  Upon  inversion  of  the  Laplace 
transforms,  for  instance,  the  nail-  and  large- tlae  lift  Units  are 


and 


jEL 

£-►0+ 


CM) 

Cu<*>' 


.  jL  *  aw  +  o«) 

l+o+  CL(*  1  2  * 


(5-23) 


+ 

(5-2k) 

It  is  of  interest  to  note  in  the  application  of  (5-23)  to 
the  plunging  case  that  the  lift  coefficient  as  is, 

using  (4-23), 


h*  ttL  *  i  tIVl  +  i  Hi)  +  oft) . 

4^04  Cu»l  ^  * 

(5-25) 

Therefore,  for  naan  time,  airfoil  apparent-mass  lift  is  dominant. 

This  follows  for  all  types  of  airfoil  motion,  since  , 

and  hence  Cl0Us)  and  Cl^Us\  have,  for  small  time,  the 

same  time  dependence  as  the  airfoil  dcsnwash,  cf.  (5-20)  and  (3-11), 
while  Cl,1$  has  time  dependence  like  the  time  derivative  of  the 
airfoil  dowrrwash,  cf.  (3-13)*  That  is,  CutUr)  and 
are,  for  small  tlae,  proportional  to  the  velocity  of  the  airfoil, 
while  Cl,  LC)  Is  proportion*!  to  its  acceleration. 

In  view  of  the  nature  of  the  Jet- flap  equations,  sene  of 
the  local  properties  of  the  flow  are  of  interest.  From  (5-18)  and 
(5-22),  it  can  he  shown  that,  upon  taking  the  inverse  Laplace  transforms, 
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(5-26) 

This  1/square-root  singularity  la  the  vortex  distribution  at 
the  trailing  edge  of  the  wake  vortex  sheet  is  the  Mathematical 
representation  of  the  "starting  vortex"  observed  by  Fraadtl  in  his 
striking  pictures  of  these  flows,  cf.  Prandtl  and  Tdetjene  (1934). 

The  downwash  behind  the  airfoil  can  be  calculated  from  (2-73)  or  (2-85), 
in  principle,  once  (i~  has  been  found.  Hobbs  (1957) 

made  such  calculations  in  order  to  estimate  the  loading  an  the 
tail  of  an  aircraft  whose  wing  is  flying  through  a  gust  field. 

Without  details  of  the  calculations,  Hobbs  gave  curves  for  the 
Wagner  problem  which  Indicate  that  the  downwash  is  finite  as 

%-+■  C4-Uoi-  on  the  wake  vortex  sheet,  has  a  square- 

root  singularity  as  %-r  c+U*4  +  off  the  sheet,  and  then 

dies  off  at  infinity  like  •/V*  In  analogy  with  the  behavior 
bf  the  downwash  induced  by  the  singularity  of  the  vortex  distributions 
near  the  leading  edge  of  a  flat-plate  airfoil,  this  behavior  at 
C+U, ck  is' consistent  with  (5-26). 

As  -fe  -►  o  ,  the  end  of  the  shed  vortex  sheet  moves 
in  towards  the  trailing  edge.  Also,  inversion  of  (5-19)  with  (5-22) 
shows  that  P(k)  **  'kf*  as  4  o  .In  the 

absence  of  circulation  around  the  airfoil,  there  must  be  flow  around 

* 

the  trailing  edge  in  the  first  Instant,  4*0*  ,  after 

initiation  of  the  motion,  in  order  to  satisfy  the  downwash  boundary 
condition  on  the  airfoil.  Presence  of  the  jet  for  Cy  >  O 


> 
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vith  the  requirement  that  it  always  remain  tangential  to  the  trailing 
edge  of  the  airfoil,  (2-57),  will  therefore  greatly  modify  the  flow 
pattern  near  %  »  C  +Uoi  by  preventing  flow  around  the 

trailing  edge  in  the  first  instant,  and  thus  for  all  later  time. 

If  the  response  to  fl  ti)  *  ^  WO  is  desired  in 
order  to  treat  the  problems  of  en  airfoil  snapped  up  to  Incidence 
and  the  impulsive  deflection  of  a  trailing- edge  flap,  (5*9)  and  (5-11), 
ft/l*)  *  ^  may  be  substituted  formally  into  (5-18)  to 

(5-21).  However,  consideration  of  ?(p)  for  this  input  Indicates 
that,  for  •  Inversion 

of  this  transform  cannot  be  carried  out  because  it  is  too  singular, 
corresponding,  as  seen  by  differentiation  of  (5-26),  to  Xlx)  71 

for  small  ~C  This  implies  .  that  although  Cf)  * 

can  be  substituted  formally  into  (5-18)  to  (5-21),  the  Integrals 
(5-14)  to  (5-16)  and  these  transforms  fall  to  exist.  Therefore  the 
cases  involving  Ptlb)~  S  ^)  are  too  singular  to  permit  a 
solution.  However,  the  step- function  inputs  sure  themselves  idealiza¬ 
tions  of  physically  realizable  inputs,  so  these  physical  inputs  can 
be  treated  by  direct  application  of  the  results  (5-18)  to  (5-21) 
with  their 

For  the  sharp-edged  gust  problem,  as  defined  by  (5-12), 
the  transform  of  flit)  may  be  found  to  be,  cf.  Sears  (1940), 
or  Miles  (1956), 

(5-27) 


where 


I.m) 


and 


z,m 


are  the  Modified 
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Bessel  Functions  of  the  First  Kind.  Substitution  of  (5-27)  into 
(5-18)  to  (5-21)  gives  among  its  results  the  fascinating  one 
noticed  by  Sears  ( 1 9^*0)  that,  for  the  Xflssner  problem,  A9! 

U(f>-  -c?e  +*,&)]  > 

(5-26) 

which  has  the  same  dependence  as  flip)  for  the  Wagner 

problem,  using  (5-19)  and  (5-22).  The  lift  coefficient  for  an 
airfoil  flying  through  a  stationary  gust  depends  on  4;  in  the 
seme  fashion  as  the  circulation  in  the  Ifegner  problem. 

The  vortex  distribution  at  the  trailing  edge  of  the  shed 
vortex  sheet  is,  upon  inversion  of  (5— 18)  with  (5-27)/ 

(5-29) 

vanishing  as  %-*»  C+Uo4 '  on  the  sheet,  hence 

continuous  at  '%»£+  Ub4  .  Hobbs  (1957)  also 

gave  curves  of  downv&sh  for  the  gust  cause  with  several  X>0 
He  found  that,  consistent  with  the  continuity  in  the  vortex 
strength  at  the  end  of  the  shed  vortex  sheet,  the  downv&sh  is 
continuous  there.  However  a  discontinuity  in  dotnv&sh  appears 
at  4L*  G+U*  (4-  ,  which  is  the  point  in  the  wake 

corresponding  to  the  arrival  of  the  gust  at  the  trailing  edge  of  the 
airfoil.  By  the  arguments  of  p.33 ,  this  discontinuity  in  downwash 
implies  a  logarithmic  singularity  in  the  vortex  distribution  at 
that  point.  That  the  discontinuity  arises  at  this  point  rather  than 
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at  Ve-t'  ,  which  corresponds  to  the  arrival  of  the  gust 

at  the  leading  edge  of  the  airfoil,  can  he  seen  by  considering 

>  (5-12)  •  function  of  tine  is  continuous 

and  has  a  continuous  slope  at  \z  *  O  ,  hut  is  continuous  with 

a  discontinuous  slope  at  ^  ^  ,  the  discontinuity  in 

slope  leading  to  the  downwash  discontinuity. 

The  presence  of  the  Jet  for  Cj  >0  in  this  problem 

will  modify  the  flow  pattern  near  't,»c4-Uo  {4-  >  since  the 

Jet  is  required  to  remain  tangential  to  the  airfoil  at  the  trailing 
edge,  (4- 69) ,  preventing  a  velocity  discontinuity  when  the  gust 
arrives  there. 

5-3  Properties  of  the  Classical  Solutions  for  Steady- State  Oscillations 
The  appropriate  functions  of  time  for  steady- state 
oscillations  of  airfoils  have  been  discussed  in  Section  2.2  Airfoils 
In  infinite  sinusoidal  gusts  have  been  treated  by  Sears  (19*M)  and 
Kemp  (1952),  but  will  be  emitted  hers.  For  oscillations  beginning 
at  t  «  -  os  ,  the  shed  vortex  wake  behind  the  airfoil  will 
be  infinite  in  length  at  the  finite  times  under  consideration. 

The  fundamental  problem  for  airfoil  motion  is  the  Theodorsen 
(1934)  problem,  in  which  the  airfoil- quasi- steady  circulation  is 
exponential  in  time;  i.e., 

u.ceto‘, 

(5-30) 


so  the  governing  Wagner  Integral  equation  is 
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Cf  *  \Vi  iwi 

J(i^)  ir/u^§  *  -Uoc  *  . 


-c  (5-31) 

This  problem  was  shown  by  Sears  (1940)  to  be  related  to  the  tfegner 
problem  of  Section  5-2  in  the  fashion  typical  of  linear  systems; 
i.e. ,  the  steady- state  response  to  an  exponential  input  like  (5-30) 
is  found  from  the  Laplace  transform  of  "the  unit- step- function 

by  fi«0)  , 

.  It  is  of 

certain  interest  to  go  through  the  solution  in  the  maimer  of  von  Kanafn 
and  Sears  (1938).  The  application  of  the  Theodore  sen  problem  to  the 
airfoil  motions  considered  in  Chapter  4  presents  no  difficulties  of 
the  type  encountered  in  the  Wagner  problem.  Considering  the  airfoil 


response  by  replacing  the  transform  variable  p 
and  multiplying  the  response  by 


m  4* 

in  plunging  motion,  for  Ft 4:)  *  £ 

Fill?)  *  iiUc’oie^f 


rut) 


Is 


(5-32) 

and  is  a  direct  application  of  the  Theodorsen  problem.  The  airfoil 
oscillating  in  pitch  and  the  airfoil  with  oscillating  trailing- edge 
flap,  both  with  F/4)  *  are  also  direct  applications 
of  the  Theodorsen  problem,  since,  from  (4-30)  and  (4-47), 


4 u|4  »  v  W 

*  7tO*(Lc  e  +  i7r«uc  u>  1%-e)  c  , 


itiit 


and 


(5-33) 


imi 


(5-34) 


respectively. 
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For  steady- state  oscillations,  all  the  functions  have 
time  dependence  a/  g  ,  therefore 


yiv  *£) .  i  iwi  eiuJli 

(5-35) 

and  the  Wagner  integral  equation,  (5-30),  becomes 

•  U«C  ,  blrtso, 

(5-3 6) 

The  integral  in  (5-36)  exists  for  S<»t  *  o  ,  but  for  pure 

oscillations,  §fo>V*o  ,  the  Integral  is  finitely  oscillatory. 

A  cannon  approach  to  this,  cf . ,  Robinson  and  Lauxoann  (1956),  has 
been  to  solve  the  problem  for  §lui\<0  and  then  argue  by 

analytic  continuation  that  the  result  is  valid  for  , 


too.  von  Kerman  and  Sears  (1938)  treated  this  in  another  way.  They 
considered  as  the  basic  equation  not  the  Wagner  Integral  equation, 
(5-5),  but  the  equation  (5-3),  fraa  which  (5-5)  can  be  derived  if 
vanishes.  Substituting  of  (2-67) 

into  (5-3X  and  urlng  (5-30)  and  (5-35)  gives 

(5-37) 

The  integral  on  the  lefthand  side  of  this  exists  for  §(<*0  £  o 
To  evaluate  that  integral,  von  Karmen  and  Sears  (1938)  used  the 
integral  representation  of  the  Modified  Bessel  Function  of  the 
Second  Kind,  cf.  Jahnke  and  Bade  (19^5),  with  identically  equal 
integrals  added  and  subtracted,  l.e., 
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(5-38) 


Ed. f f  erentiating  this  with  respect  to  ^utc^U0  ,  and  using  the 
identity 


(5-39) 

the  integral  in  (5-37)  can  he  evaluated  to  give 

Jiurt  *  -  iwc  t TO 

(5-1*0) 

where 


£1%)*  .  S'“’U0' 

(5-M) 

is  the  Sears  Function,  and  was  tabulated  by  Kemp  (1952).  The 
circulation,  also  from  (2-67),  la 


m *  U.c e““W’ 

(5-W) 

The  circulatory  lift  coefficient  may  be  calculated  from  (3-20), 
using  (5-35),  there  being  no  difficulties  with  convergence,  giving 

Ct.Uh  Cujt)  «  2  (U%,)  eluH  9M  *0 , 


(5-43) 


Ill 


where 


C(%) 


§(w)*o , 


(5-44) 

is  the  Theodorsen  Function,  tabulated  in  detail  by  Luke  end  Dengler 
(1951)- 

Of  particular  interest  for  the  jet-flap  analysis  is  the 
shed  vortex  distribution  in  the  wake,  given  by  (5-35)  and  (5-40), 
i.e. , 

^\srw)6o, 

(5-45) 

which  has  the  above-mentioned  oscillatory  properties  at  infinity  if 
§(u>)  -  o  .  The  dcwnvash  behind  an  airfoil  performing 

steady- state, purely  oscillatory  t  §  (uj)  " O  ,  plunging  motion 

was  calculated  by  Lapin,  Crook  shanks  and  Hunter  (1 952) .  Their 
results  were 


]> 

ijwr  c 

"*  [*./«>*.($£)]  4  Me  L 2£i>‘  T5NJ( 


(5-46) 


where 


rf  i4 


(5-47) 
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is  a  function  defined,  discussed  and  tabulated  by  Schwarz  (19^3)  • 

A  few  Halting  cases  of  (5-46),  for  which  siaple  egressions  may  be 
found,  were  given  in  their  paper.  For  (A) -wo  ,  (5-46) 
reduces  to  the  airfoil- quasi- steady  downwash  distribution, 
of  (4-13)  •  In  the  limit  of  tl)  oo 


r 


+  oM  t  wc+Uf, 


toK.t)  J 


|  -  [i  |i  (^)  d  **  +  /it-c)*4  *  ->  W*  %**■  t+g. 


(5-48) 

The  latter  expression  of  (5-48),  although  not  given  by  them,  may 
be  found  using  equation  (30*)  of  Schwarz  (1943).  Away  from  the 
trailing  edge,  the  downwash  increases  like  as 

10 -wee  .  finally,  for  arbitrary  CAJ  ,  but 

%-*eo  ,  Lapin,  Crook  shanks  and  Hunter  (1952) 

gave  -  they  have  neglected  a  minus  sign  - 


ur/%,^)  *- 


irk*  «oc  e 


MS&+  M»f£\ 


ILwco 

> 


(5-49) 


This  result  may  be  shown  to  be  the  downwash  induced  at  a  point  due 
to  a  dqUbly  infinite  distribution  of  vortices  of  strength 


a.7rtk*«uce 


i«U>- 


«-*«) 


(5-50) 
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This  is  the  vortex  distribution  shed  by  an  airfoil  in  plunging 
oscillations  as  seen  by  application  of  (5-45)  to  this  case, 
using  (5-32) .  The  dowuwash  and  vortex  distributions  in  the  wake 
oscillate  both  in  space  and  tine  far  downstream.  For  o  , 

these  results  damp  to  zero  as  %-*"  ,  and  it  is  also 

clear  to  see  from  these  results  uhy  ^  (to)  ">o  oust  be 
excluded  from  consideration. 


5.4  Reduction  of  the  Equations  to  the  Steady  Jet- Flap  Theory 
The  steady  jet-flap  theory  is  also  included  in  the 
above  formulation.  It  will  be  the  limiting  solution  of  the  equations 
for  £  -rcP  in  the  transient  case  if,  PIW-I 

For  steady  oscillations,  it  is  the  limiting  solution  as  W  O 
Formally,  however,  the  steady  equations  can  be  recovered  simply  by 
setting 


and 


p/4)  *0 . 


(5-51) 


Using  (2-1)  and  (2-7),  these  reduce  (2-69),  (2-73),  (2-78)  and  (2-79), 
for  example,  to  the  set 


Jfht)  «  -  or'/*.)  ,  c.4.  , 


(5-52) 
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with  either 


C*  tx-ao, 

(5-53) 


Uric*)  *  UToU'i+UoZ* 


(5-54) 


or 


MU< 


t-% 


•  -2U«T# 


%-+Ci 

(5-55). 

The  assumption  (5-51)  imply  that  b  '►co  much  faster  than 
OL-rw  ,  therefore  the  circulation  in  the  system  is  no  longer 
zero,  and  equations  (2-81 )  to  (2-83)  sod  (2-85)  to  (2-88)  are  not 
valid.  Likewise  the  lift  and  pitching-moment  coefficients  are  not 
given  by  (3-30)  and  (3-39)-  The  lift  coefficient,  using  (5-5) 
in  (3-10),  (3- 11)  and  (3-13),  and  using  (5-52)  in  the  integral  over 
til)  ,  becomes 


Cl-Cl.  - 

C 

The  pitching-moment  coefficient  in  similar  fashion  reduces  from 
(3-29)  to 


(5-56) 
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Cm  *  **£<.#  -  (• 4 C,l 

(5-57) 

Eliminating  %l%)  between  (5-52)  and  (5-53) >  the  single 


equation  treated  by  Spence  in  I  and  II 


ur/%) 


w*i%)  * 


(5-58) 


with  either 


or 


ur(c+)  ~  UTcCc)  4-UoZm  t 


(5-59) 


is  found. 


4JJcZg 
CCy  , 


(5-60) 


In  I,  Spence  assumed  in  the  Jet- deflection  problem 
(Section  4.1  above)  that 
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AtBST-'J KM] 


(5-61) 

where 


(5-62) 

The  first  tern  was  introduced  to  satisfy  (5-60),  the  series  being 
regular  at  %  *  C  •  Having  substituted  this  expression  into 
(5- 53),  the  resulting  equation  was  approximately  solved  by  a  collocation 
scheme,  i.e.,  the  first  H  coefficients,  Aw  ,  were  found 
by  satisfying  the  equation  exactly  at  H  points.  This  approximation 
converged  rapidly,  and  detailed  results  for  N-  9  were  found. 

Spence  also  solved,  in  I,  the  incidence  problem  (Section  4.3  above 
with  assumptions  (5-51))  in  the  same  fashion,  having  assumed 

(5-63) 

Finally,  Spence  (1953)  solved  the  Jet- augmented- flap  problem 
(Section  4.4  above  with  assumptions  (5-51))  by  collocation,  where 

UJI%)  .  - 

^ar  ***©  ^ 


(5-64) 
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Later,  In  II,  Qpence  solved  (5-58)  to  (5-60)  analytically. 
In  the  Jet- deflection  and  incidence  problems,  for  jX*- 1  , 

where 


This  was  done  by  transforming  the  downwash  distribution  by 


(5-65) 


orW  -  Ooz.fefrjp)  +  iM.[l  +  ijgf  p. (Ji)]  . 

(5-66) 

A  detailed  discussion  of  this  transformation  and  its  implications 
will  be  made  in  Section  6.1  for  the  full  equations  of  the  unsteady 
problem.  It  suffices  here  to  say  that  substitution  of  (5 -66)  into 
(5-58),  with  Ufc lx)  given  by  (4-29)  fander  the  assumptions 

(5-51))#  to  an  iterative  expansion  of  £*■(')£)  and 
of  the  forms 


(5-67) 

and 

fd  M  *  [j  +■  A 


(5-68) 
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where  &  is  Euler's  constant.  Equations  for  the  functions 
,  ftl'yt)  ,  P%.lt!)  and  ^Ix1) 

were  found  to  be 

L  ®o  ,  o*%'**> 

pm  » •, 


(5-69) 


M-d&fW. 


6*% U  «o 


f , lo)  -O  , 


(5-70) 


where 


L P%1%) '  ~ 


:)-  -  4%’^  ,  «**.'*•  1 

P%,lo) »  J 

t  6*1 tf* 


(5-71) 


J 


(5-72) 


w .  no  -  Jfrf  £0 ' 


(5-73) 


119 


These  equations  were  then  solved  by  a  Mellin- Transform  procedure, 
a  technique  to  be  discussed  in  Section  6.2.  The  lift  coefficients 
found  from  this  solution  agreed  very  well  with  the  numerical 
solution  of  I,  even  up  to  IK  -  ) 


CHAPTER  6  -  "BOUNDARY- LAYER"  NATURE  OP  THE  PROBLEM:  TRANSFORMATION 
OF  THE  EQUATIONS  TO  "BOUNDARY- LAYER  "COORDINATES 

There  are  many  problems  in  applied  mechanics  which  may 
be  typified  as  "boundary- layer"  problems.  The  name  arises  from 
the  mathematical  similarity  of  these  problems  to  Prandtl's  (1904) 
"grenzschicht, "  or,  as  it  has  been  translated,  "boundary- layer," 
in  the  flow  of  viscous  fluids  around  bodies.  Properties  and 
examples  of  "boundary- layer"  problems  have  been  discussed  by 
Carrier  (1953);  Friedrichs  (1955);  sod  others. 

Carrier  (1953)  states  three  requisite  criteria  for 
consideration  of  a  problem  as  of  "boundarylayer"  type.  First, 
the  coefficient  of  the  highest- order  derivative  appearing  in  the 
equations  of  the  problem  should  be  very  small  compared  to  0(1)  . 

Second,  the  other  important  terms  in  the  equations  should  have 
coefficients  of  0(0  •  Finally,  the  domain  of  the  problem 

in  the  coordinate  system  chosen  must  be  characterized  by  lengths 
of  0(0 

Spence  first  realized  in  II  that  the  equations  for  a 
jet- flapped  airofil  in  steady  motion  were  of  "boundary- layer"  type, 
provided  the  Jet  momentum  coefficient,  Cy  ,  was  sufficiently 
small.  In  most  practical  applications  proposed  for  the  Jet- flap 
the  required  Cj  are  within  the  limits  of  the  restriction. 
Spence  then  made  the  transfQzmation  to  "boundary- layer"  coordinates  - 
equation  (5-66)  of  Section  5.4  -  and  solved  the  steady  problem  within 
that  framework.  For  the  unsteady  problem,  Spence  made  a  similar 
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approach  in  III.  It  Is  important  to  discuss  the  unsteady  problem 
in  detail  with  respect  to  Carrier's  criteria  in  order  to  determine 
the  validity  of  this  approach. 

First,  the  highest- order  derivative  appearing  in  the 
problem  is  always  on  the  right-hand  side  of  the  particular  dynamic 
interaction  equation,  (2-69)- (2-72),  used.  This  derivative  is 
multiplied  by  JX  ,  (5-65),  so  if  the 

highest- order  derivative  is  o(l)  .  As  mentioned  in  Section  2.7 
this  is  the  only  explicit  appearance  of  A  in  the  equations 
of  the  problem,  hence  Carrier's  second  criterion  is  also  met. 
Considering  the  problem  in  the  form  of  the  single  equation 
in  jfl'hk.S  or  ,  (2-90)- (2-93),  these  criteria 

are  again  seen  to  be  met.  Therefore  these  two  criteria  apply  in 
the  unsteady  problem  exactly  as  in  the  steady  one. 

The  characteristic- length  considerations  of  the  domain  in 
the  unsteady  problem  are  somewhat  different  from  those  of  the  steady 
one.  In  the  steady  problem,  Spence's  numerical  solution  in  I 
clearly  indicated  that  the  downwash  and  vortex  distributions  die 
off  rapidly  within  a  few  chord  lengths  behind  the  airfoil.  Thus 
the  chord  length,  c  ,  is  a  significant  characteristic  length. 

In  the  discussion  of  the  classical  unsteady  transient  solutions 
in  Section  5.2,  however,  discontinuities  in  the  downwash 
distributions  auad  the  associated  singularities  in  the  shed  vortex 
distributions  were  ..oted.  In  particular,  U»4.  and  Uo (-L - 
sire  significant  lengths  for  adrfoil  motion  auad  gust  penetration, 
respectively,  in  addition  to  the  edrfoil  chord,  C 
Similarly,  equation  (5-^5)  indicates  that  for  steady  oscillations 
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the  wave  length, 


ajrU 

a/c 


,  of  the  shed  vortex  distribution 


Is  an  Important  length.  Therefore,  for  times,  Jam 


,  of 


ujc 
X7TV* 

the  airfoil  chord,  C 


0(1)  or  smaller  In  the  transient  case,  or  for  frequencies, 

,  of  OO)  or  larger  in  the  oscillating  case, 

,  is  again  a  sufficient  characteristic 
length  since  U.fc  ,  UcU-  an d  2^1 

are  of  the  same  order  as  C  •  Furthermore,  (2-100)  and  (2-102)  show 
that  the  dovnvasb  and  vortex  distributions  die  off  at  infinity  even 
more  rapidly  than  for  the  steady  case,  for  X>  .  For 

times  so  large,  or  frequencies  so  small,  that  Uot,  * 

are  much  greater  than  C  ,  a  "boundary- 


and 


■ZirUo 

u>c 


layer"  based  on  C  alone  is  inadequate,  since  important  effects 
would  be  well  outside  the  layer.  Spence,  in  III,  treated  the 
region  near  X  *  C  +  Uot  as  of  primary  importance  for  large 

times  in  the  transient  case,  along  with  the  region  near  the  trailing 
edge. 

Realizing  the  time  or  frequency  limitations  just  discussed, 
the  appropriate  transformation  of  the  coordinates  for  small  A 
is,  following  III, 


(6-1) 


4-  . 


(6-2) 


and 


tt> 


(6-3) 
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The  leading  edge,  %aO  ,  transforms  to  , 

so  for  vanishing  A  ,  the  airfoil  becomes  seal- infinite  in 
length.  In  conjunction  with  the  coordinate  transformation,  it  is 
convenient  to  follow  Spence's  procedure  of  IZ  and  III  and  transform 
the  dependent  variables  representing  the  Jet  by 


MM  * yu  ZU.c A  f  KA‘) , 

r  2U0A 

y/*,±1  *  B/%,1)  yxcA  Ji/*U' ) , 

and 

«.  u.A  l%fkh!,V) . 

The  time  dependence  transforms  by 


(6-4) 


(6-5) 


(6-6) 


(6-7) 


(6-8) 

Determination  of  the  non-dimensional  constant  A  ,  and  the 
function  will  be  discussed  after  transformation 

of  the  equations. 

The  equations  relating  down wash  and  Jet  ordinate,  (2-1), 
and  potential  difference  and  vortex  distribution,  (2-7),  transform  to 
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and 


MtU)  ■ 


ah/*',*)  .  >W?)  -  a 
’>«  ” /  a^j^ar)  , 


(6-9) 


(6-10) 


Trnasforaing  the  various  fonus  of  the  equations  of  the  Jet, 


main -stream  dynamic  interaction,  (2-69)  to  (2-72),  gives 


3t£&g  +  MteLji  -  ..  aHfcfl 

>*•  /  ajt'* 


+/A1_  M&Ljf  tkhUL 

/  I*?* 


»  a  O0j 

(6-n) 


££(%£) 


•42 


aw 


+ 


an-*  J 


a*-Jl 


3  MU*) 

jtiyfifr-j 


oo. 


(6-12) 


/sii/ViO 

+/  TTf^yP  , 


O*  CO  J 

(6-13) 


and 
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1)C,S  t  u\y*) 


ZktU') 


+,,*  iskh^e) 

/  iSyw*  j 

o £■  iid.  oo, 

(6-14) 


provided 


9 xWt.t)  _  „ 

_ —  s  O 


(6-15) 

The  two  forms  of  the  downwash  equation,  (2-73)  and  (2-85) > 
transform  to 


+  £a$Tu,,iw'‘  •***•* 


and 


(6-16) 


fc/*!,*')  *  S5,  (if  o^». 

T) 


(6-17) 


The  various  forms  of  the  trailing- edge  boundary  conditions, 


(2-7 6)  to  (2-80),  become 


K(o+,k')  * 


yfa<a-Bfo.i> 

^*-cA 


(6-18) 


%c 


A 


9 


(6-19) 
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=  uIa  + 


Dt 


A  » 


(6-20) 


and 


while  the  Wagner  integral  condition,  (2-82),  becomes 


(6-21) 


00 

f  . .  r.'tf) 

J  f,ya  V&hp* , 


(6-22) 


where 


p.w.  r.'ti'). 

The  lift  coefficient,  (3-1 5) >  and  the  pitching-moment 
coefficient,  (3-35) >  become, under  this  transformation, 

00 

Cut'**  c.,w) 


(6-23) 


(6-24) 


and 


C* It')  *  -K*,/e)  +  kgClJe) 


(6-25) 


80 


The  function  3(%>^)  can  be  determined  in  each  case 
that  the  factor  ($r»  (it/O4  no  longer  appears  in 

the  terms  &[£?[“* (%ti) '  *&?] 

of  (6-1 6)  and  (6-17)*  Uhls  gives  a  first-order  partial  differential 
equation  for  Bf'h-L)  which  can  be  Integrated  subject  to  the 
conditions  that 

and 


ZYle-.l)  -ZBfai)  _ 

5> 


o, 


or 


u*oCt+M  -  ^'.0, 

which  make  the  boundary  conditions,  (6-18)  to  (6-21),  dependent 
only  upon  the  Jet  deflection,  To  fe)  .  lhe  constant  A 
is  chosen  to  simplify  the  form  of  the  final  equations. 

The  equations  given  above  are  exact.  For  small  yU>  , 
however,  neglecting  terms  of  0(^~)  enables  the  equations  to  be 
greatly  simplified.  The  chief  simplifications  are  in  the  kernels 
of  the  integrals  of  the  downvash  equations,  (6-1 6)  and  (6-17),  and, 
as  will  be  seen  below  for  the  particular  problene  of  Chapter  4, 
in  the  inhomogeneous  terms  on  the  right-hand  sides  of  these 
equations.  These  simplifications  are  equivalent,  in  the  first 
approximation,  to  neglecting  JX  X  with  respect  to  I  ,  which 
is  -valid  only  near  %■'  -O  ,  i.e.,  within  a  "boundary  layer" 

near  the  trailing  edge,  and  to  neglecting  with  respect  to  I 

valid  only  near  4!  ~  O  ,  i.e.,  within  a  "boundary  layer" 


near  the  time  origin.  Higher  approximations  in 
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may  be  obtained  in  the  unsteady  problem  by  expanding 
the  solution  in  terns  of  JX  as  Spence  did  for  the  steady 
problem  in  II,  as  mentioned  in  Section  5.4. 

The  equations  of  the  various  problems  of  Chapter  4  will 
now  be  written  in  the  transformed  coordinates. 

Jet-Deflection  Problem:  The  input  to  this  problem  is 
the  boundary  condition  on  the  Jet  slope  at  the  trailing  edge. 
Therefore,  choosing 


A  =  To 


and 


(6-26) 


6  *o , 


(6-27) 

the  important  boundary  conditions  (6-19),  (6-20)  or  (6-21)  become 

^  5*1. 


(6-26) 


This  is  the  problem  formulated  by  Spence  in  III. 

Plunging- Alf roil  Problem:  Here,  A  and  Q/%,4^ 
are  found  to  be 


A* 


(6-29) 


and 
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ho c  Flh) . 

The  important  terms  are,  then,  in  the  dovnvash  equations  (6-1 6) 
and  (6-17), 


(6-30) 


&  *  y**'  *£r , 0 


4.  v.‘e.«o 


(6-31) 


and 


’"‘"i 


(6-32) 


in  the  boundary  conditions,  (6-18)  to  (6-21), 


*  kM)  =  -  0  ; 


(6-33) 


and  in  the  Wagner  integral,  (6-22), 


(6-34) 


The  leading  inhomogeneous  term  in 


r 


in  the  equations  is 


the  right-hand  side  of  the  Wagner  integral  condition  if  the  dovnvash 
equation,  (6-1 6)  vlth  (6-31 ) ,  is  used,  or  ^^)^<»/*.  if  the 
form  of  the  dovnvash  equation  incorporating  the  Whgner  condition, 
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(6-17)  with  (6-32),  is  used.  The  second  term  in  (6-32),  i.e., 
the  leading  term  in  (6-31 ),  is  proportional  to  jX  because  it 
is  a  higher  power  of  .  It  would  be  necessary,  then,  to 

carry  the  solution  to  the  second  approximation  in  JX  to  account 
completely  for  the  airfoil-  quasi-  steady  do®  wash.  The  first 
approximation  will  give,  however,  the  leading  term  in 
the  lift  and  pitching-moment  coefficients. 

Pitching- Airfoil  Problem:  Here  it  is  found  that 

r 


*  < 


55  .*** 


-  d* 


for 


(6-35) 


and 


The  important  domwash  terms  are,  fran  (6-1 6)  and  (6-17)> 


(6-36) 


rSfc 


did') 


e.**M 


3 

T 


A$W> yi x [x.* M  +  V* } ]  , 

(6-37) 


and 
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U»A 


♦**«»*#] 

+f'gZ)  fp**>  <■**  <^rJ  ,  e.*y* 


&  *%!  *40  * 


yt,lfr&klt*dp]  >e-* 


(6-38) 


the  boundary  conditions,  (6-18)  to  (6-21),  are 


w.^1  -  -  *<*«) .  -£jf(&  **&  •  °  i 


(6-39) 


and  the  Wagner  integral  condition,  (6-22),  is 

' e** 

1 T&W)  ,  . 


(6-40) 


The  equations  of  the  first  approximation  in  are  thus 

Identical  to  those  for  the  plunging  airfoil,  if  6. 

The  higher  approximations  in  JX  satisfy  different  equations, 
however,  including  the  £. » case  which  is  really 
a  second  approximation.  Here  the  analysis  must  be  carried  to  the 
third  approximation  in  jx-  to  account  for  all  the  airfoil- 


quasi- steady  terms.  The  first  approximation  is  now  clearly  restricted 
to  very  small  b  ,  since  certain  terms  in  (6-37)  and  (6-38) 
are  higher-order  in  jf-  because  flh)  is  higher-order  in 
yU  than 

Blown- Flap  Problem:  As  the  down  wash  equations  (4-50) 
and  (4-51 )  stand  for  this  case,  they  are  inconvenient  for  calculation; 
however  the  arctangent  may  be  expanded  like 

U-[Ub  (&fj  o 


valid  If  Cet  \  ijh'p*  .  Provided  1b  not  too 

small,  use  of  this  expansion  is  consistent  with  the  neglect  of 


(6-41) 


compared  to  | 


in  the  first  approximation.  Therefore 


a  first  approximation  should  be  obtainable  in  this  manner.  Furthermore, 


n  \=  TT 

e*o 


the  equations  agree  with  the  pitching  problem  for 
.  In  this  approximation, 


As  - 


JjF  w 


(6-42) 


B/x/fc)  *  p*  cc*S$r)  F lb) . 


(6-43) 


The  downwash  terms  including 
and  (6-17), 


Ol/A 


only  are,  from  (6-1 6) 


(ifTi* it,*)-  -  m  4^ 


(6-44) 
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and 


*12 •*#♦/!! 


,i{^ 


♦  «•»**«*  ; 

(6-45) 


the  boundary  conditions,  (6-18)  to  (6-21),  sure 


U,t) .  ^  .  k.Hi'1  -  -SL^lgifp  .  o  ; 


(6-46) 


and  the  Wagner  integral  condition,  (6-22),  becomes 


oo 


f  4(S',£)<U  _  jj- «f jrAtevofta) 

J-^vT  -W~JF  V  JBSS&SsP^  *  *}P  • 

0  (6-47) 

The  equations  for  the  first  approximation  in  JX  are  again 
identical  vlth  those  of  the  plunging  and  pitching  airfoils.  For 
the  first  approximation  alone,  the  expansion  (6-4l)  is  not  used, 
since  the  leading  tern  in  (6-45)  comes  from  the  airfoil- quasi- 
steady  circulation  which  is  exact.  Again  the  approximations  involved 
clearly  limit  the  validity  of  the  first  approximation  to  very  small 

i!  . 


Sharp- Edged- Gust  Problem:  The  dcwnwash  equations  (4-67) 

and  (4-69)  are  also  not  in  convenient  form  as  they  stand.  An 
expansion  for  small  -fc  ,  and  small  %-c  ,  may  be  made. 

From  such  an  expansion,  it  is  found  that 


A* 


(6-48) 


and 


6 W»o  . 


(6-49) 


The  dovnvash  terms  including 
and  (6-17), 


only  are,  from  (6-1 6) 


^j^]*  o**uoo , 


(6-50) 

and 


±  /*' 
U*A 


n/ofl 

Tit  J 


4*  -  **  *'*]  *  ^ 


the  boundary  conditions  are,  (6-18)- (6-21), 


(6-51) 


and  the  Wagner  integral  condition,  (6-22),  is 


(6-52) 


*  tt^  -y*- 


(6-53) 
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5he  equations  for  the  first  approximation  here  are  the  same  as  for 
the  above  three  cases  if  their  tine  dependence  is  explicitly  made  to 
be  .  The  approximations  have  the  sane 

implications  as  to  their  validity  as  above. 

From  consideration  of  the  above  five  problems,  it  is 
clear  that  in  the  first  approximation  for  small  yti  and  t/  , 

there  are  tvo  fundamental  unsteady  jet-flap  problems  to  be  solved. 

The  first  of  these  is  the  Jet- deflection  problem  first 
treated  by  Spence  in  III.  The  problem,  in  terms  of  the  downvash, 

Rli,^  ,  and  the  jet  vortex  distribution,  ^  hf-W)  , 

for  example,  is  governed,  using  (6-4)  to  (6-23)  and  dropping  the 
primes  for  convenience,  by 

(6-54) 

and 

(6-55) 

or 

-  if 

(6-56) 

with  either 


W<K,  t)  ■  ?(i), 


(6-57) 
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(6-58) 

and 


(6-59) 


The  problem  can  also  be  formulated  in  terms  of  the  jet  ordinate, 

,  and  the  potential  difference, 
using  the  first  approximations  to  (6-9)  and  (6-10), 


and 


(6-60) 


**  , 


(6-61) 


and  the  appropriate  equations  from  (6-1 1 )- (6-23), 

To  the  first  approximation  in  /*  ,  then,  the  lift  and 

pitching-moment  coefficients,  (6-24)  and  (6-25),  for  the  Jet- deflection 
problem  become 


(6-62) 
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and 


U  (M  »  -3(i-  jfWJi . 


(6-63) 

It  immediately  follows  that;  to  this  order;  the  center  of  pressure 
is  at  the  three- quarter- chord  point;  a.»/f  .  For  steady 

flow  Spence  found  numerically  in  I  (and  unpublished  calculations 
using  the  results  of  II  confirm  it)  that  the  center  of  pressure, 
to  first  order  in  ,  is  at  the  half- chord  point.  The 

validity  of  the  "boundary- layer"  transformation  for  small  time  only 
is  pointed  up  here;  since  the  center  of  pressure  remains  at  the 
three-quarter  chord  point  for  all  time  in  this  approximation. 

The  second  fundamental  problem  in  the  first  approximation; 
corresponding  to  the  airfoil-motion  problems  and  the  sharp-edged 
gust  problem  (if  ^76)®  )  is,  in  terms  of  j b  (%,<=) 

and  , 


V 


C*-X.*V3J 


(6-64) 


and 


c*-y.**j 


(6-6  5) 


or 


O 


(6-66) 
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with  either 


\ 


kKD  9 cj 


(6-67) 


or 


and 


(6-68) 


The  lift  and  pitching-moment  coefficients,  (6- 2k)  and 

oo 

Cm-  $[&$**** , 


(6-69) 

(6-25),  are 


(6-70) 


and 

CmU)  -  Ujtb  Cm,  IB  +  ^  j&tyw1*, 

(6-71) 

where  Cl,  l D  ,  £  Atoll)  ,  and  CM,l-b)  are  appropriate 

to  the  particular  problem. 

The  solutions  of  the  two  fundamental,  sznaU  l  , 

problems  are  closely  related.  To  explore  this  relationship  fully 
in  the  transient  case,  it  is  necessary  to  consider  Laplace  transforms 
taken  on  the  time  variable.  Using  the  definition  of  the  Laplace 
transform,  (5-17)  >  and  the  initial  conditions  that  the  flow  is 
undisturbed  prior  to  initiation  of  the  transient  motion,  expressed  by 
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(6-72) 

the  Laplace- transformed  equations  of  the  jet- deflection  and 
airfoil- motion  problems,  (6-54)  to  (6->59)  and  (6-64)  to  (6-69) ,  are 


and 


(6-79) 
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4  (6-80) 
where  ^X(Ptp)  and  ktyp)  are  the  transformed  vortex 
distribution  and  downwash  for  jet  deflection; 

rt;,) 

and  )£■  (x.ff>)  those  for  airfoil  motion,  and 

the  transform  of 

Differentiating  (6-77)  and  (6-78)  with  respect  to  % 

gives 


(6-81) 


(6-82) 

where  the  identity  &  9  irPh  +  +  has  been  used. 

By  inspection  of  (6-8l)  and  (6-82),  the  deflection  solution  of  (6-73) > 
(6*74)  and  (6-76)  suggests  that  the  airfoil-motion  solution  may  be 
written  as 


(6-83) 
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and 


a£]*a) 


Mp\lxh.)p  t 


(6-84) 


Integrating  these  with  respect  to  7L>  ,  and  using  the  t railing-edge 

boundary  condition,  (6-79) >  and  the  condition  that  the  vortex 
distribution  vanishes  at  infinity  downstream  gives 


(6-85) 


where  Alp1)  has  been  evaluated  from  the  Wagner  integral  condition, 
(6-80),  upon  integration  by  parts  and  use  of  (6-83)*  To  complete 
the  solution,  the  integral  required  for  the  lift  and  pitching-mcment 
coefficients,  (6-70)  and  (6-71),  is,  upon  integration  by  ports. 


(6-87) 


Therefore,  once  the  Jet-deflection  solution  has  been  found,  the 
airfoil-motion  solution  follows. 

The  Laplace- transformed  equations  of  the  transient  problems 
also  hold  for  the  problems  of  steady-state  oscillations,  i.e. ,  where 


4  lit)  *  J  * 


«?4 


(6-88) 
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and 


(6-89) 

provided  the  transfora  variable  p  is  replaced  by  (if) 
and  ?Cp)  by  I  .  Thus  the  relations  (6-85)  to  (6-87) 
with  the  proper  modifications,  also  hold  for  steady- state  oscillations 
in  jet  deflection  and  airfoil  motion. 


CHAPTER  7  -  CRITIQUE  OF  ATTEMPTED  8QLUTI0RB  IH  "BOUNDARY- LATER"  COORDINATES 


7-1  Critique  of  Spence's  Solution  of  the  Jet- Deflection  Problem  for 
Snail  Time. 

Spence,  in  Section  3  of  III,  approximated  to  the  equations 
of  the  transient  Jet-deflection  problem  for  small  time  after  a  unit- 
step-  function  input,  and  then  obtained  a  "solution”  of  the  simplified 
equations.  In  the  present  section  it  will  be  shown  that  this 
"solution”  is  incorrect,  and  that  no  valid  solution  of  Spence's 
approximate  equations  can  be  found. 

The  approximation  of  III  is  to  neglect,  for  «it»n  time, 
the  derivatives  with  respect  to  'JL-  compared  to  the  derivatives 
with  respect  to  i.  in  the  convective  derivatives  which  appear 

in  the  problem;  i.e.,  it  is  assumed  that 


2. 


which  is  valid  for 


(7-1) 


'L  *  oft) 


(7-2) 

This  holds,  for  small  time,  everywhere  except  so  near  the  trailing 
edge,  O  ,  that  and  ^,*01^) 

Since  the  important  boundary  condition  must  be  evaluated  at  the 
trailing  edge,  it  must  be  satisfied  there  in  a  sense  consistent  with 
(7-2) .  The  small-time  solution,  if  found,  clearly  will  be 
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non- uniformly  valid  near  the  trailing  edge.  Very  near  the  trailing 
edge,  on  the  other  band,  if  and  hence  2^  »  *  ( ^ 

the  equations  reduce  to  those  for  the  steady  Jet- deflection  problem 
treated  in  II  (and  discussed  briefly  above  in  Section  5*4).  The 
steady  solution  is  valid  near  the  trailing  edge  in  a  region  growing 
with  time. 

Instead  of  the  downvash,  Spence  treated  the  Jet  ordinate, 

,  as  the  unknown  function,  along  with 
Cross-differentiation  of  the  small-time  approximation  of  the 
appropriate  equations,  (6-13)  and  (6-54),  and  elimination  of 
CJ^  gives  the  equation  in  alone, 

denoting  by  subscripts  (  ^  and  (  }  partial 

differentiation  with  respect  to  those  variables, 

0  (7-3) 

The  form  of  (7-3)  implies  that 


9 


i.e.,  that 

■x*0 U*)' 


which  from  (7-2)  implies  that 


(7-4) 


(7-5) 


(7-6) 
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Instead  of  solving  (7-3)  for  k  1%,^  directly 
In  terms  of  a  similarity  variable,  ,  as  in 

III,  it  is  more  instructive  to  retain  %  and  k  explicitly, 
and  to  solve  simultaneously  for  and  4)  while 

retaining  the  similarity  approach.  Improper  treatment  of  the 
unit- step  function  appearing  in  the  problem  is  one  reason  for  the 
incorrect  result  of  III. 

If  ,  the  equations. (6-13),  (6-55), 

(6-18),  (6-19),  (6-58)  and  (6-59)  are,  in  this  approximation, 

=  -  htt*  Iv-M  i  o<-%*-co 

(7-7) 

and 

*  (7-8) 

with  either 


h  »  0 


(7-9) 


and 


or 


*  1  , 


(7-10) 


«- 1  it) 


9 


(7-.l») 
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and 


00 

jr*^1  m< 


Following  III,  but  In  terms  of  'JL  and  t 

of  ^  ,  define  the  Mellln- transform  pairs 

C4i0. 


c-i«* 


and 


oo 


C4»9» 

c-i*> 

Ml  ^'^Gls)  -  At, 


(7-12) 

instead 


(7-13) 


(7-14) 


(7-15) 


(7-16) 


where  C  is,  as  yet,  an  undetermined  real  constant. 

From  the  form  of  (7-13)  and  (7-15),  the  assumption  of 
an  *Xy^Vs  similarity  is  evident.  For  ***/■£**  <■  I  , 
the  functions  h  !%£}  and  ^  /%,k\  can  be  evaluated  in 

series  in  increasing  powers  of  ^/f^s  ,  [ 

etc.  by  moving  the  line  of  integration,  5»C  ,  to  the  left 

past  the  singularities  of  MlS)  and  &  ts)  ,  respectively,  in 
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the  S -  plane-  Because  of  (7-2) ,  such  series  would  be  valid 

In  the  United  range  where  4^*%  ,  i.e.,  where 


4;*  ol%)  *o[oU*/t)]  m 


(7-17) 


In  a  similar  manner  for  >  I  >  h  l%jb)  and  ^ />,,-/>) 

can  be  evaluated  In  series  in  Increasing  powers  of 

,  etc.  by  moving  the  line  of 

integration,  S  »  C.  >  to  the  right  past  the  singularities 

of  /-//S')  and  6>/&)  >  respectively,  in  the  g-  plane. 

These  expansions  are  valid  for 


’■*  o(%)  . 


(7-18) 


By  successive  differentiation  of  (7-13)  with  respect  to  *%■  , 

cuV o 

=  -  iFt  [ i  <«  s  H/s)  Js , 

and,  provided  $/S+i)/Sf»)  W/s)  is  regular  in  the  infinite 
strip,  C-3<£  fc/S)  *  C  ,  the  line  of  integration  can  be  moved 

to  the  left  to  S  -  C-3  .  Writing  g  =-2-3  and 

dropping  the  bars  then  gives 

C+Ho 

hwh.i)  *  -  is-!)  Hls-1)  Js . 

CriM 

(7-19) 

Substituting  (7-15)  into  the  Integral  operator  in  (7-8), 
interchanging  the  order  of  integration,  and  using  the  integral 
given  in  Section  3.1  of  III  (where  the  factor 
advertently  omitted),  namely. 


was  in- 


(7-20) 


x 

tr 


gives 


ifliffl't-'P  ~  iw  ea)Jtj 

°  (7-21) 


if 


W“i 


Next,  differentiating  (7-13)  with  respect  to  X  gives 
C+Uo 

kM  -  s/i>]  #«**, 

CrVa  (7-22) 

vhe:  ■*.  is  the  Dirac  delta  function,  (5-10).  For  m  <■  c  , 

1>y  (7-17)  and  t.  is  necessarily 
greater  lhan  zero.  Therefore  the  term  may  be  neglected 

since  it  is  zero  for  O  .  For  fu%)  >  t  a 

toy  (7-18)  and  *L  is  not  necessarily  greater  than  zero.  Thus 
the  Hi)  tern  must  be  considered.  By  the  definition  of  s/«  , 
however,  it  follovs  that  if  €  'S'O  ,  /i)  •  O 

Therefore  the  %ll)  tern  may  be  neglected  provided  the  first 

pole  of  His)  to  the  right  of  fc/sWc  occurs  for 

(Ji/si  >  -  I  *  This  is  clearly  satisfied  here  since 

|c|<.^  ,  so 

CtlSD 

J,*M  *  H/s)Js. 
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Differentiating  (7-15)  with  respect  to  gives 

WM  •  ♦***»*»]««*. 

L-iM 

(7-24) 

By  the  same  arguments  given  for  )l* /*,<:)  ,  the  $IL) 

term  may  be  dropped  if  the  first  pole  of  Cj  l S')  to  the  right 

of  *C  occurs  for  {Lli\>  ^  Assuming 

this  to  be  satisfied  gives 

cttto 

*  ^jic*4WM*i<«  [->¥)  is . 

‘•to  (7-25) 

Finally,  differentiating  (7-13)  and  (7-15)  with  respect 
to  4C-  gives 


CfUP 

c-Ue 


(7-26) 


and 


C.+J* 


c*b>» 


(7-27) 


Substituting  (7-19)>  (7-21),  (7-23)  and  (7-25)  into  the 
equations  (7-7)  and  (7-8)  gives 


04  Wo 


s  *  o 


C-l* 


(7-28) 
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and 

c4W 

J -  f  4T* 6&)  +[^s)  H/$$  Js  -  o . 

C.in 

(7-29) 

These  will  be  satisfied  for  arbitrary  <%  and  £  provided 

<SM  =(2$f)otwm. 

(7-30) 

Prom  (7- 13)#  it  follows  that  the  trailing- edge  boundary 
condition,  (7-9)  >  will  be  satisfied  if  and  WM 

is  regular  at  5  *  0  .  The  input  trailing- edge  boundary- 

condition,  (7-10),  will  be  satisfied,  using  (7-2 6),  if 
has  a  simple  pole  with  unit  residue  at  5  *  - I 

By  inspection  of  (7-16),  the  Itagner  integral  condition, 
(7-12),  will  be  satisfied  provided  <sr±i=o  . 

Summarizing,  (7-30)  must  be  solved  for  H&)  and  (y(s) 
subject  to  the  conditions  that 

S  (s+Ms-iz)  Hfs)  is  regular  in  the  infinite  strip, 

-3* (Us)* 

(7-31) 

is  regular  in  the  infinite  strip, 

(7-32) 


is  regular  at  S  -  O  , 


(7-33) 


sVl 


near  S  «  - 1 
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i 


i 

(7-3*0 

and 


G(i\  -  o. 

Uie  difference  equation,  (7-30),  for 
be  satisfied  if 


(7-35) 

W/S)  will 


Hfe)  - 


f^rM 


Af/$)  Vis) 

> 


which  implies  that 


(7-36) 


Gtsv  -V'lay  a/tvsMtiiyis'i  t 

(7-37) 

where  AVs)  is  some  function  of  period  three,  and  Vl$) 
is  the  function  introduced  and  discussed  by  Spence  in  III. 

The  Spence  function,  V(S)  ,  satisfies  the 

difference  equation 

V(s)  tWns  *  o  # 


and  may  be  vzltten  as 


(7-38) 


(7-39) 

where  Qj%  (s)  is  the  Lighthlll  function  discussed  by  Spence 
in  II.  This  latter  function  satisfies 


152 


G,'(s)  -  Wits  =  o 


and  is  represented  by  the  infinite  product; 


(7-to) 


Cs0 Is)  *Go(~ i~s) 


■  U  t  <<%<<-£)} , 


(7A1) 


S  s  kl*  j 


i.e. ,  (yp  l S ) 
and  at 

and  at  S  *  -  w  -  ^ 
and  S  *  -4^ 
interest,  say 
at  Sr  i,  4,  ^ 

S* -1 3 

Using  these  properties, 
to  be  regular  in  -  3  ^  ft/s)  <-  4 


has  poles  of  order  h  at  S' «  -  W 

,  zeroes  of  order  H  at  S » n 

,  and  equals  unity  at  S  *  O 
Therefore,  in  the  region  of  primary 

,  ijJ/s)  has  simple  poles 
and  simple  zeroes  at 


fr-/)(  w/cv 

/*+*!){  T  ' 

,  except  for  the  required 


is  seen 


simple  pole  at  S  •  -  1  ,  and  a  double  pole  at  S  *  -  Z 

where  only  a  simple  pole  is  permitted.  Also,  Cfkws  tyls) 

is  regular  in  0*-%(s)*-l  ,  but  6/4)  ^  O 

Therefore  the  function  M  /s)  must  have  a  zero  at  S  *  “  2. 
to  satisfy  the  regularity  condition,  (7-31  ),  and  a  zero  at  S!  4 
to  satisfy  the  Vagner  Integral  condition,  (7-35) • 

To  Insure  the  existence  of  the  Mellin  transforms,  (7-13) 
and  (7-1 5),  it  is  necessary  that  ( '*'/■£**)  H/s)  and 

('*’/■£**)  SQ>&  be  integrable  as  f  $&))  — ►  Co 

along  .  In  this  limit  ( %/&* ) 

is  bounded,  as  is  C6&  ITS 


rS 


.  Using  the  appropriate 
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form  of  Stirling's  formula,  (S~/)j  +*  |9/s)| 

(a^J)/C  /M'1**'*  e**"* 

imf'**-* e1**' 

in  II  and  III  that  <flS)  ~ 


««•*  -IN 

,  and 


Spence  found 
,  vbile 


~  .  Therefore, 

(*/&*)'*  H(s)  ~llL/&*)*&(s)  -v  ,  which  is 


integrable . 

The  requirement  that  AMs')  introduce  zeroes  at 
S  -  -2.  ,  S  -  ^  ,  could  be  met  in  only  two 

ways.  In  the  first  of  these  it  would  be  at  the  expense  of  introducing 

corresponding  poles  into  H  (s)  and  Cr(S)  In  violation 

either  of  the  regularity  condition,  (7-31 ),  or  of  the  regularity 
condition,  (7-32) .  The  second  would  be  at  the  expense  of  introducing 

exponential  factors  making  )-)  IS)  and  Qls)  fail  to  be 

integrable  along  (R(s)  *■  C  .  Therefore,  no  function  Mis) 

with  the  desired  properties  can  be  found.  Thus  no  solution  of  the 
approximate  set  of  equations,  (7-7)  to  (7-12),  can  be  found.  Discussion 
of  reasons  for  this  will  follow  the  details  of  Spence's  "solution" 
of  III. 


Spence  obtained  his  "solution"  in  III  by  choosing  M/s) 


as 


Mis)  - 


4 

3  (-1)6.  Hi 


fr*TT*r 


(7-42) 

where  the  numerical  factor  is  chosen  to  give  the  pole  of  WlS) 
at  S  *  -  I  unit  residue.  The  zero  required  at  X  *  -  2. 

is  introduced  at  the  expense  of  a  simple  pole  of  Mis) 
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at 


S’  i 


S-'l  ,  etc.  While  the  simple  pole  at 
does  not  violate  the  regularity  of 

In  , 

has  a  simple  zero 


His) 

(7-31),  since  V/S> 

there,  the  condition  on  the  regularity  of  G/S)  In 

0  JUS)  6:  2.  >  (7-32),  is  violated  hy  the  simple  pole 

of  G(S)  introduced  at  S  *  2.  by  Mis) 

Therefore,  the  neglect  of  the  sw  tezm  in  (7-24)  was  not 
justified.  Furthermore,  of  course,  the  Wagner  integral  condition, 
(7-35) >  is  violated;  in  fact  the  integral,  (7-12),  does  not  even 
exist  with  the  G<  s)  found  here.  On  these  two  grounds,  the 

"solution"  is  clearly  invalid. 

Investigation  of  the  properties  of  this  invalid  "solution" 
helps  indicate  the  failure  of  this  attempt  to  approximate  the  full 
equations  for  small  time,  and  is  thus  worthwhile .  Since  the 
behavior  of  Hfs'J  G(S)  in  the  whole 

S  -  plane,  in  particular  their  poles,  is  known  from  (7-3 6), 
(7-37)  and  (7-42),  series  expansions  for  h  (t,h)  and 
iua/  be  found  to  be 

f  *1IU  ‘  ,  **** 


*  < 


and 


(7-43) 
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&/*/&•)  +  °^*),  *** 


(7-44) 

(The  third  term  in  the  small  -  '*/i£*4  expansion  for 

|r|  /%£)  has  the  sign  corrected  from  equation  (38)  of 

III.)  Similarly,  the  first  terms  in  the  expansions  for  the  important 
derivatives  and  integrals  are,  from  (7-19),  (7-21),  (7-23) ,(7-24), 

(7-26)  and  (7-27), 


*•(*$&).  *>t*. 


(7-^5) 


f  4  oil)  , 


(7-46) 
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*  ^ 


(7-47) 


4 

JH)!^ 


♦  0/&)  ,  x^* 


3gj4) 
uK  %** 


(7-1*8) 


r 


g. 

3irt 


|  -  +  O £ X'vH*i*  /*/£**)]  a  xyi***% 


and 


(7-^9) 


-  j^eij^.  Aa/x/^*)  +  ^V|. 


(7-50) 


157 


From  (7-1*5)  and  (7-50)  it  is  seen  that  the  approximate 
equation,  (7-8);  1b  satisfied;  from  (7-1*7)  and  (7-1*8)  it  is  seen  that 
the  approximate  equation,  (7-7);  is  not  satisfied  for  large 

,  the  extraneous  %lis)  tern  in 

(7-1*8)  arising  since  it  was  incorrectly  neglected  in  (7-24) .  The 
boundary  conditions  at  the  trailing  edge,  (7-9)  and  (7-10),  are 
satisfied,  considering  (7-43)  and  (7-46' ,  however,  the  alternate 
statement  of  these  conditions  in  terms  of  the  vortex  distribution,  (7-1 1); 
is  not  satisfied,  considering  (7-50) • 

The  lift  coefficient  was  evaluated  in  III  using  (6-62), 

which  is 

CM  *  - jUj f'pmH . 

(7-51) 

It  is  important  that  in  the  derivation  of  this  form  of  the  lift 
coefficient,  satisfaction  of  the  Wagner  integral  condition  was 
required.  Since  this  condition  is  violated  here,  (7-51)  should 
not  have  been  used.  Furthermore,  (7-44)  shows  that  the  integral  in 
(7-51)  does  not  even  exist.  In  III,  the  calculation  was  actually 
carried  out  by  taking  inside  the  integral  in  (7-51)  and 

using  (7-7),  the  ensuing  integral  existing  only  because,  as  seen 
above,  (7-7)  is  not  satisfied  by  the  "solution"  of  III. 

As  a  final  property  of  the  "solution,"  the  assumptions 
that  and 

will  be  checked.  For  / 

valid  in  the  range  b  ^ 


,  the  expansions  are 
,  (7-17).  If 
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4*  oorO 

o  fir*) 

r  % 


(7-52) 

where  fC  »l  is  some  constant,  this  Inequality  Is  met. 

Furthermore,  the  Inequality  ,  (7-6),  Is  also 

satisfied.  Substituting  (7-52)  into  (7-^5),  (7-^),  (7-W)  and  (7-^9) 
indicates  that 


0(f) 

(7-53) 


pM)  *  ocf) 

(T.*, 

Therefore  as  assumed,  but  hv/*,4) 

is  of  the  same  crder  as  kt  /y.,t)  ,  inconsistent  with  the 

assumption  that  kt /*,£)  ■»  .  In  a  similar 

fashion  for  >/  ,  the  expansions  are  valid 

for  %» ,  (7-18),  such  that  the  inequalities  (7-2) 
and  (7-6)  are  also  satisfied.  These  are  if 


i-  Q  (K'%) 

%*  O  (fc*0 


(7-55) 
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where  j£>7(  is  again  some  constant.  Substituting  (7-55) 
into  (7-1*5)  >  (7-1*6),  (7-W3)  end  (7-1*9)  indicates  that 


ht/%£)  --  ol&hti 

> 

k%  1%M  *  o(K  ***) 


and 


(7-56) 


OlK?x)  ,  4=o 
O  1&JM,  t>« 


*  O(^) 

(7-57) 

Therefore  /%,-t)  ■»  as  assumed.  At  4  "  O  , 

as  assumed;  however  the  solution 

should  also  be  valid  for  4  >0  provided  (7-6)  is  satisfied, 

but  ^t/jt-4)  ,  and  the  assumption  breaks  down. 

The  question  arises  whether  the  unit- step- function  time 
dependence  assumed  here,  being  discontinuous  at  4=0  ,  is  too 

singular  to  permit  a  solution.  However,  if  the  general  input, 


'f(i)  *  / 


.'4*1/4) 


K  >0 


,  is  treated  in 


precisely  the  above  fashion,  again  no  solution  can  be  found  that 
satisfies  the  regularity  condition  on  S  (S+i)(s*t)  H(s) 

satisfies  the  regularity  condition  on  the  new  G/S)  and 

satisfies  the  Wagner  integral  condition. 
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This  approach  to  the  solution  of  the  small-time  Jet- deflection 
problem  clearly  is  inadequate.  Consideration  of  the  result;  (7-53); 
indicated  that  came  out  to  be  of  the  same 

order  as  kt  (%tl)  >  indicating  that  it  should  not  have 

been  neglected.  This  does  not  seem  surprising  in  retrospect;  since 
the  important  input  to  the  problem  is  the  boundary  condition;  (7-10); 
on  h%/r  ,D  •  Neglecting  with 

respect  to  In  M)  in  the  downwash  equation;  yet 

requiring  satisfaction  of  the  boundary  condition  on  k*  M 
seems  inconsistent.  Ibis  suggests;  as  an  alternative;  formulation  in 
terms  of  the  downwash;  k  fat-.t)  ,  and  the  vortex  distribution; 

jj  If-iVl  ,  still  making  the  approximation  ^  ^ 

for  small  time.  The  next  section  will  deal  with  this  formulation 
and  its  subsequent  failure  to  give  a  valid  approximation  to  the 
complete  set  of  equations  for  small  time. 

7.2  Further  Critique  of  Small-Time  Approach:  Jet- Deflection  and 
Airfoil-Motion  Problems. 

Formulation  in  terms  of  the  downwash  on  the  Jet,  k  (%,±\ 
and  the  vortex  distribution  representing  the  Jet;  ^ /%,i)  , 

by  equations  (6-54)  to  (6-5 9)  in  the  Jet-deflection  problem;  and  in 
the  airfoil-motion  problem  by  (6-64)  to  (6-69)  eliminates  the 
contradiction  of  satisfying  the  input  boundary  condition  on 
while  neglecting  it  with  respect  to  kt  l%,4)  in  the 

equations.  The  small-time  approximation  of  neglecting 

•%-  derivatives  with  respect  to  t  -  derivatives; 

wherever  convective  derivatives  appear;  will  be  retained.  Therefore; 
the  various  order-of-magnitude  restrictions;  (7-1 )>  (7-2)  and  (7-4) 
to  (7-6)  still  hold. 
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The  approximate  small-time  Jet- deflection  equations  are 
then;  from  (6-54)  to  (6-59);  denoting  this  approximation  by  a 
subscript  (  )  , 


and 


with 


f 


(7-58) 


j£  <%t)  ’ 


0  f 


(7-59) 


(7-60) 


*0. 


(7-61) 


Equation  (7-59)  is  exact,  so  only  in  (7-^)  bas  it  been  necessary 
to  approximate  for  small  time. 

To  treat  the  problem,  again  define  the  Mellin- transform  pairs 


(7-62) 


(7-63) 
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and 


t+lm 

»  jjJV'i*  iU)  ^‘lS)  1,5 

1M  K  l s)  *  ^ ** 


(7-61*) 


(7-65) 


where 


give 


C  is,  as  yet,  cm  undetermined  real  constant. 

Successive  differentiations  of  (7-62)  with  respect  to  dt 

CtUfi 

fej  /^,i)  *  "  £irl  (V’$"*  ^  -L^)  6Cf+0(s*S'z)^* /s)  » 

w '  J 

c-i* 


(7  66) 


Substituting  (7-64)  into  the  integral  operator  of  (7-59)  and  evaluating 
it  as  in  Section  7.1  gives 

c«k* 

c-t» 


(7-67) 

if  /  cl  *.  4^  .  Differentiating  (7-64)"  twice  with  respect 

to  4  gives 

c.«t# 

fs/t-v  * 

c-i* 
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but,  provided  Is  regular  in  the  Infinite  strip, 

C  Qli)  <■  C*3  ,  the  arguments  of  Section  7*1  Justify 

neglect  of  the  S(l)  and  J  lb  terms,  so 

t ♦*» 

« ii]Vw  ±*  1 «)  fc  ♦¥)  0 *■  H)  £tsn)i$ . 

C-K* 

Substitution  of  (7-62),  (7-66),  (7-67)  and  (7-68)  into  equations 

(7-58)  and  (7-59)  gives  ‘ 
ctu* 

(7-69) 


(7-68) 


C-*> 


and 


C4|A0 

-  Uwstfts)}  Js  *  o  , 

c-U* 

(7-70) 

which  will  be  satisfied  if 

*.Tfe)  --  tiftH&'a)  *  ■ 

(7-71) 

Fran  (7-67)  it  is  seen  that  the  trailing- edge  boundary- 
condition,  (7-60),  will  be  satisfied  if  O  C  *  £  and 

has  a  double  pole  at  the  origin  with 
coefficient  ^  .  The  Wagner  integral  condition,  (7-61), 

will  be  satisfied,  using  (7-65),  if 

Summarizing,  (7-71 )  may  be  considered  a  difference 
equation  for  Qcf  ($)  which  must  be  solved  subject 


to  the  conditions  that 
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(s)  is  regular  in  the  Infinite  strip,  o  fc/s)  j  , 


(7-72) 


and 


&.Z(S)  -  near  S  *  O  , 


4>.TH)  =  o. 

The  difference  equation,  (7-71 ),  will  be  satisfied  if 

fS-aj  SWiirV 


(7-73) 


(7-74) 


Ao(s)  =  ^ 


($)! 


yv/s)V(s-3), 


(7-75) 


which  implies  that 


Z 


r^1  m  v,3)  > 


(7-76) 


where  Nl$)  is  again  seme  function  of  period  three,  and  tyls) 
is  Spence's  function,  (7-39)*  For  ,  the 

conditions  (7-72)  and  (7-73)  are  both  satisfied.  However,  the 
requirement,  (7-74),  to  satisfy  the  Wagner  integral  condition  is 
not  met.  As  in  Section  7*1>  if 

and  is  integrable  as  /  9/S))  along  fe/s)  *  C 

Introduction  of  a  zero  at  S>^  by  a  function  fills') 

of  period  three  would  either  Introduce  a  pole  into  the  strip 

0<d  (Us)  <■  ^  in  violation  of  the  regularity  condition 

on  Alls)  ,  or  would  introduce  an  exponential  factor 
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making  tike  Mellin  transform  of  ($)  fall  to  converge 

as  /  0  (s)  I  along  Oils)  »  C  •  Therefore,  a  solution 

satisfying  all  the  conditions  cannot  be  found.  Evaluating  the 
integral,  (7-61),  for  AMS)  2 1  ,  using  (7-65),  (7-75), 

(7-40)  and  (7-41)  gives 


jV'p.iW, .  ^  #im , 


(7-77) 

clearly  greater  than  zero.  If  the  order- of-magnltude  considerations 
of  the  previous  section  are  observed,  the  assumptions  of  neglecting 
and  vith  re6Pe°t  to 

is  satisfied.  Therefore,  failure  to 
obtain  a  solution  here  results  solely  fran  the  inability  to  satisfy 
the  Wagner  integral  condition  in  this  approximation. 

The  approximations  f  (%,t)  and  k»Z 
may  be  considered  as  the  first  terms  in  an  expansion  of  the  full 
equations,  (6-54)  to  (6-59).  This  expansion  may  be  written  in  terms 
of  the  Mellin  transforms  as 
C4k» 

--  lit)  A, , 


(7-78) 


and 


k'M  .  f  tWlWXW 

c-|e 


(7-79) 


The  second  approximations,  A,Tw 
are  the  transforms  of  (%,-i) 

say,  which  satisfy 


and  %X(S) 
and 


and 


•0 


with 


-  o 


and 


e» 


*  O. 


for 

and 


3 
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(7-80) 


(7-81) 


(7-82) 


(7-83) 


Operating  on  the  Mellin  transforms  (7-78)  and  (7-79) 
and  fei *  l%,l)  as  for  (%,i) 

ko  ,  the  inhomogeneous  difference  equation, 


%T(s)  =  UirsKu  -  (SgU>a ' 

(7-84) 

must  he  solved,  with  Q£.C&  ^  ,  Qff  l S)  regular  in 

o±Us)£  \  ,  #,T  /s)  having  no  stronger 


a  singularity  at  the  origin  than  a  simple  pole  in  order  to  satisfy 
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the  downwash  condition,  (7-82),  and  -  O  to  satisfy 

the  Wagner  integral  condition,  Die  solution,  which  again  satisfies 
everything  but  the  Wagner  integral  condition,  is 

*4  H. 

(7-85) 

and  the  corresponding  W(s)  .  Evaluating  the  integral, 

(7-83),  using  (7-78)  and  (7-85)  gives 


*  -§2p  i*  1«« 

*  ‘  (7-86) 

Therefore  the  second  approximation  in  £  not  only  fails  to  correct 
the  error  in  the  first  approximation  of  failing  to  satisfy. the  Wagner 
Integral  condition,  but  introduces  further  error. 

To  complete  the  picture,  the  small-time  approximation  will 
be  applied  to  the  airfoil-motion  problem,  as  defined  by  (6-64)  to 
(6-69j^  and  the  solution  attempted.  The  time  dependence  here  will  be 
assumed  to  be  to  avoid  the  difficulties  with 

the  Dirac  delta  function  as  discussed  in  Section  5.2.  With  the  same 
small-time  approximations  of  (7-1),  (7-2)  and  (7-4)  to  (7-6),  the 
equations  simplify  to 

-  kZh,i)  , 

0  w- 


and 


(7-87) 


(7-88) 
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with 


fo+,4)  • 


5  O 


(7-89) 


and 


Jf *  7^1^) 


(7-90) 

The  important  input  to  the  problem  is  the  Wagner  Integral  condition, 
(7-90),  as  contrasted  to  the  trailing- edge  boundary  condition, 

(7-60),  in  the  jet-deflection  problem. 

In  view  of  the  relation  of  the  airfoil-motion  solution  to 
the  jet- deflection  solution  as  pointed  out  in  Chapter  6,  it  is 
impossible  that  an  airfoil-motion  solution  of  this  type  will  be  found. 
However,  it  is  instructive  to  carry  through  ohe  Mellin- transform 
approach  to  see  exactly  how  it  fails. 

Assuming  a  similarity  solution  in 
Mellin- transform  pairs  are 


the 


c-ieo 
CO 


(7-91) 


(7-92) 


and 


(7-93) 
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/ 


***•% Ai  Kis) *  Jt . 


(7-94) 


Three  %-  derivatives  of  (7-93)  give 
c+iv 

ktM  *  - ;&j . 


C-l® 


(7-95) 


Substitution  of  (7-91)  into  the  integral  operator  in  (7-88)  and 
evaluating  it  as  usual  yields 

*»  c«U» 


c-Wp 


(7-96) 


if  Id^-  2,  •  Differentiating  (7-91)  twice  with  respect  to 

■fc-  ^  and  dropping  the  £(■£)  and  ^  ^ t  terms  by 

the  arguments  of  Section  7-1 ,  provided  Is)  is  regular 

in  the  infinite  strip,  C4-3  ,  gives 

ctwo 

$*,*) *  &  JV~*  i a)  (**■;  i i?)  . 


6-Ko 


(7-97) 


The  equations  (7-87)  and  (7-88)  are,  upon  substitution  of 
(7-93)  and  (7.95)  to  (7-97), 


C4l« 


%'s'3  s  /st/)/s**)%f  /s)  +(^)  s 


and 


£-iuo 


C«i« 


0 

(7-96) 


4* 


C-  iso 


(7-99) 
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which  will  be  satisfied  if 

tfffl  -  u  Vm  *  TOS  if  *•>. 

(7-100) 

Hie  trailing- edge  boundary  condition,  (7-89) ,  will  be 
satisfied  if  O and  if  Is)  has  a  singularity 
at  the  origin  no  stronger  than  a  simple  pole.  The  Wagner  integral 
condition,  (7-90),  is  satisfied,  using  (7-92),  if  &,*«.)*  TT 
In  summary,  the  difference  equation,  (7-100),  must  be 
solved  for  [s)  subject  to  the  conditions  that 

i8  regular  in  the  infinite  strip, 


(7-101) 


and 

8£li)*TT. 


(7-102) 

The  difference  equation,  (7-100),  will  be  satisfied  if 


1-fiP 


and 


(7-103) 


P/3)  v/s) , 


(7-104) 
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vhere  P(s)  is  s  function  of  period  three  sod  Is 

Spence's  function,  (7-39).  Kith  Pf$)  S )  ,  S  ftj/s) 

is  regular  in  the  strip,  Ot-t-li)*-  J  ,  snd  B  ^ 

However,  A.a  (s)  has  a  double  pole  at  the  origin,  removal  of 
which  by  a  suitable  function,  Pis)  ,  would  again  either 
introduce  a  pole  into  the  strip  0  £fs)£  £  ,  violating 

the  regularity  of  4tt(s)  ,  or  would  add  an  exponential  factor 

ng  b'  Li)  fail  to  converge  as  /  §ls)j  -vao  along 

(jt(i)  s  C  •  Thus,  as  expected,  a  solution  satisfying 

all  the  conditions  of  the  problem  cannot  be  found.  The  double  pole 
at  the  origin  of  gives,  in  terms  of  fee  ll.i)  , 

from  (7-93)  and  (7-104) 


(7-105) 


which  is  not  zero,  but  in  fact  singular,  at  -fc  =  0 

Introduction  of  a  continuous  time  dependence  at  *  O 
by  =  J)?  £*  d./i)  for  both  the  Jet- deflection  and  airfoil- 

motion  problems  lead  to  the  same  failure  as  above  to  yield  solutions. 
Therefore,  it  is  not  a  question  of  the  motion  being  too  singular 
at  O 


It  has  been  seen  that  in  both  the  Jet- deflect! on  and 
airfoil-motion  problems,  no  solution  using  the  present  small-time 
approximation  can  be  found.  For  the  Jet-deflection  problem,  the 
input  boundary  condition  at  the  trailing  edge  was  satisfied,  but 
not  the  Wagner  integral  condition.  On  the  other  hand  for  the 
airfoil-motion  problem,  the  input  through  the  Wagner  Integral 
condition  was  satisfied,  but  not  the  t railing-edge  boundary  condition. 
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The  condition  left  unsatisfied  In  each  Instance  is  an  important 
physical  condition  on  the  problem,  and  a  "solution"  which  falls 
to  satisfy  it  is  meaningless.  Moreover,  the  approximation  does  not 
change  the  order  of  the  equations,  since  the  highest  derivative, 

(or  if  Mfctt  is  eliminated,  M-*'  ) 

is  retained,  so  there  is  no  Justification  for  failure  to 
satisfy  one  of  the  conditions.  Failure  to  get  a  solution  in  this 
small-time  approximation,  i.e.,  the  failure  to  get  a  similarity 
solution  to  the  full  equations  for  small  time,  must  arise  frcm  the 
inability  of  this  approach  to  give  a  valid  result  near  the  trailing 
edge.  As  discussed  in  the  beginning  of  this  section,  any  solution 
foundvould  have  been  non- uniformly  valid  in  time  near  the  trailing 
edge.  Furthermore,  as  expressed  by  (7-17),  a  similarity  solution 
for  small  *4  would  be  limited  in  validity,  because  of  the 
small-time  approximation,  to  the  range  where 

This  restriction  apparently  prohibits  satisfaction  of  the  equations 
at  X*0  ,  and  results  in  the  fed. lure  of  the  present  approach. 

A  recent  private  communication  from  hr.  h.  A.  Spence  indicates 
that  he  has  solved  the  small  Jet-deflection  problem  exactly, 

using  Laplace  transforms  on  both  the  %  and  -4  variables.  The 
exact  leading  lift  term  for  small  time  has  been  found.  His  results 
also  indicate  -chat  the  solution  near  the  trailing  edge  is  not  expressible 
in  terms  of  the  similarity.  These  results  will  be 

reported  shortly  by  Dr.  Spence. 
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7-3  High-Frequency  Steady-State  Oscillations;  Jet  Deflection  and 
Airfoil  Motion  Problems. 

No  similarity  solution  in  could  be  found 

for  small  time  In  the  previous  sections,  yet  the  same  approach  for 
high  frequencies  of  steady- state  oscillations  appears  to  lead  to  an 
similarity  solution  satisfying  the  equations  and 
all  the  conditions,  still  approximating  by  neglecting  X- 
derivatives  with  respect  to  -t-  derivatives.  In  Section  5.2 
of  III,  Spence  found  high-frequency  "solutions"  to  the  Jet- deflection 
problem  in  terns  of  the  Jet  ordinate  and  vortex  strength.  Although 
these  "solutions"  do  not,  of  course,  encounter  the  step- function 
difficulty,  they  do  fail  to  satisfy  the  Wagner  integral  condition, 
hence  are  incorrect. 

The  equations  for  the  Jet- deflect! on  problem  in 
steady  oscillations,  with  !)(?)  ^  O  ,  are,  from  (6-73)  to 
(6-76),  using  (6-88)  and  (6-89)  and  the  remarks  following  them, 


(itf  'f  (vV  +2  ft?)  pW  +  J 


(7-106) 


and 


04  00, 


(7-107) 


with 


»  i 


t 


e 


(7-108) 


17^ 


and 


and 


f)  il  x  o 


(7-109) 


In  full  generality,  Mellin  transforms  may  be  defined  by 

c«-W» 


ft;))  *  jrJV*  K'(s;»ds' 

r°  \ 

kz(s;))~  \j%s"kx(rJ)<h- 


otieo 


d-*D» 

6xis\1\  ty) 

6 


Using  the  techniques  of  the  preceding  sections,  the 

various  derivatives  and  integrals  are 

c-U  co 

tjp®.  -jSijV*'1  sfc+i)tf«)K*fs;Ws  ; 

C-ioO 

A 

and,  provided  G* [$]?)  is  a  regular  function  of  ^ 

in  the  infinite  strip,  C*H(s)4C+3  , 

c*i*e 

oT/*;J)  *  lj|V’'3^Vnj;J)  Js  , 

c-Uo 


(7-110) 


(7-m) 


(7-H2) 


(7-113) 


(7-1 14) 


(7-H5) 
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6Vs»; M dJ, 


and 


c-i» 


CHS 


(7-116) 


•  ^JV5*3  (s+/)  Ish)C,x(sh ;?)<}$  ; 

fc-V» 


and,  if  jc/A^  , 


(7-117) 


-fef  _ l  L- 

TTJt  SJ  u  f-v  -  5Ti  ^ 

0 

£-tt* 


*+l* 

1^$fciKirsGl^V)<is . 


(7-H8) 


Substituting  these  into  the  governing  equations,  (7-106) 
and  (7- 107),  gives 
C4 too 

-  i(i1)lsix)Gx(s+rJ)  +  ts+i)ls+z)C,T{su;?) 

-  s(s+i)/s+*)  fcT($;?)J  ds  9  o 


C-ioo 


(7-119) 


and 


C4l« 

S[KZIS;?)  ~  'binfi.'lTS  GXl*j})*h  ds  a  O  , 


c-lfl0 


■»rtiioh  will  be  satisfied,  if 


(7-120) 


•  4*»sSVs;f)  .  — -|^6Vs4»;J)  4  6>;». 


(7-121) 
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The  trailing- edge  boundary  condition,  (7-108),  is  satisfied 
if  OC.C *■  2,  and  G  tS’jl)  has  a  double  pole 

with  coefficient  Jj,  at  S  =  O  The  Wagner  integral 

condition,  (7-109),  is  satisfied  if  =  O 

Summarizing,  the  difference  equation,  (7-121),  must  be 

A 

solved  for  G  (S'r)  subject  to  the  conditions  that 

Gzls;  is  regular  in  the  infinite  strip,  j 

(7-122) 


Gj1  (Sj?)  =^,.near  S»0, 


(7-123) 


and 


SVi;P) »  o . 

Examination  of  (7-121)  indicates  that  an  approximation 
scheme  in  ?  ,  for  large  J  ,  may  be  made  .by  expanding 

ICV  (S'?)  and  £Vsj  ?)  in  terns  of  ^  by 


(7-124) 


(7-125) 


and 


6T(t;n»  ****)]. 


(r-126) 


This  approximation  is  precisely  the  same  as  that  of  the  previous 


- 
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sections;  i.e.,  it  is  assumed  that  (7-1)  is*  here, 

(7-127) 

The  same  limitations  in  the  regions  of  validity  of  the  solution  hold 
here,  (7-2  ),  (7-5),  (7-&  (7-17)  and  (7-18),  vith  1/f  replacing 
.  Equations  (7-121)  to  (7-124)  may  thus  be  written 

Ko  (S)  IS)  +o(?'v»)  *  iu^Trs  6*  (s)  +i?'</*6k'iTs6]c/s) 

=  -  Gl  (in)  -  i^[6^sn)  6* /sfri]  f  otf-**), 

(7-128) 

with 

Goli)  +OfP"^>)  regular  in 

(7-129) 

-±-  jjS?ls)  *5/^]  *  ^  near  S  *  O , 

(7-130) 

and 

fir}-**  [6?  (i)  *  if*  +  Of}-*'*)]  =  o . 


Equating  like  powers  of  then  gives  the  iteration  scheme 

for  successive  approximations  in  I? 


(7-131) 
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The  first  approximation  requires  solution  of 

jCo  -  ($)  *  td*.TW  Go  ($)  -  -  6*  (S+3)  , 


(7-132) 


subject  to  the  conditions  that 

(S)  is  regular  in  the  infinite  strip,  O*-  \  f 


near  S*Oj 

and 

GM)*  o  ; 

the  second  approximation  requires  solution  of 

fc,T(S)  *  t^TTS  -  -  <ZUS+*)  -  Z  fis’/S+l)  , 

subject  to  the  conditions  that 


(7-133) 

(7-13*0 


(7-135) 


(7-136) 


5  (p*  (s^  is  regular  in  the  infinite  strip,  Oi  *•  \  , 


(7-137) 

and 


O 


(7-138) 
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The  difference  equation,  (7*132),  of  the  first  approximation 
is  satisfied  if 


G-M'-*zkui)  ***•**"«>• 


(7-139) 


Furthermore,  all  the  conditions,  (7-133)  to  (7-135) >  are 
satisfied,  including  the  Wagner  integral  condition.  Using  the 
results  of  Section  7* 1>  fS rl)!&»ts)  as  |9^S))**«o 

along  ljl($)-C  .  The  complex  may  he  written 

as  |&fe))-*oo  along  ,  GX(Sj/)"l9^l  ^ 

as  j  VS)|  — ►  ce  along  02.1S)  -  C  .  Therefore,  since  the 
Mellin  transforms,  (7-1 10)  and  (7-1 12),  must  he  integrahle  along 
{LfS»)sC  for  0£-£*"z.  ,  it  is  necessary  to  restrict 

.  This  is  a  mild  restriction  for  the  case  of 
steady- state  oscillations,  eliminating  only  purely  divergent  motion. 

In  a  like  manner,  the  inhomogeneous  difference  equation, 
(7-136),  of  the  second  approximation  is  satisfied  if 


— -  Y>/s-i) , 


■**  9GoHW|) 

(7- 1^*0) 

where  the  first  term  is  the  particular  solution  to  the  Inhomogeneous 
equation,  (7-136),  and  the  second  term  is  the  general  solution  of 
the  corresponding  homogeneous  equation.  Both  conditions,  (7-137) 
and  (7-138),  are  satisfied  and  the  Mellin  transform  is  again  integrahle 
if  *  5 
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Higher  approximations  GJ(s)  ,  etc.  could  also 
be  found,  but  are  not  treatedhere.  Also,  higher  approximations 
In  could  be  found  following  the  technique  of  Spence  in  II 

for  the  steady  problem,  but  these,  too,  are  emitted  here. 

The  lift  and  pitching-moment  coefficients  can  be  evaluated, 
using  (6-62),  (6-63),  (7-H3),  (7-126),  (7-139),  (7-140),  (7-39) 
and  (7-1+0),  giving 


Cl  *  2  r°  { 1  +  +  0 

Cm  *  1%-ACu . 


(7-141) 


(7-142) 

As  derived  in  Chapter  6,  the  solution  for  airfoil  motion 
is  related  to  that  for  Jet  deflection  by  (6-85)  to  (6-87).  Here, 
in  particular,  using  (7-1 13), (7-126),  (7-139)/ (7-1 40),  (7-39)  and  (7-40), 

(7-1>t3) 


and 


ir,))  <>***$£ V,U,‘ 


f 

(7-144) 


while  the  important  integral  in  the  lift  and  pitchlng-mcment 


coefficients  is 


(7-145) 


These  results  check;  of  course;  with  those  found  by  a  direct 
Mellin-  transform  solution  of  the  airfoil-motion  problem. 


For  the  plunging  airfoil;  from  (6-70);  (6-71),  (4-23) ;  (4-24), 
(4-25)  and  (7-145), 


(7-146) 


and 


4  . 

(7-147) 

The  leading  tern  in  both  lift  and  pitching-moment  coefficients  is 
the  apparent-mass  contribution  and  is  lowest- order  in  both  ^ 
and  yu»  •  The  next  tens  is  circulatory  in  nature  and  is  the 
leading  jet  effect.  The  last  term,  also  circulatory,  may  be 
identified  as  precisely  the  leading  circulatory  tens  for  the 
classical,  ^-0  ,  limit.  That  is,  the  leading  jet- induced 

term,  proportional  to  ,  is  lower- order  than  the  leading 

classical  circulatory  tens,  which  is  recovered,  however,  in  this 
approximate  solution. 

With  like  interpretation,  the  lift  and  pitching-moment 
coefficients  for  the  pitching  and  blown-flap  airfoils  may  be  written. 

For  the  pitching  airfoil  they  are 


*2-Tra.(A 


l-l« 


f  /!-$)* 
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* 


Cm  - 

.5  +  . 


(7-149) 


for  the  blown- flap  airfoil, 


4  *  %$)'{[-*?*  ®Vv  * 

+ 1  [\0-2CoxJ  +w\/i-u»o]^57i  -  i[l  ftr  ***  v»0  v 


+  Oi 


(7-150) 


and 


+  jS^-^[xs//-xo*\) +  +  $t.V+Ss^ 

+  o^'4/’)} 

(7-151) 

It  appears  then,  that  in  the  sense  that  (it)) 
and  are  neglected  in  (7-106),  a  solution 

of  the  problem  of  steady- state  oscillations  can  be  found  satisfying 
all  the  equations  of  the  problem,  Including  both  the  trailing- edge 
boundary  condition  and  the  Vagner  Integral  condition.  Whereas 
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there  was  failure  to  solve  the  transient  problem  for  small  time, 
solution  to  the  oscillating  problem  has  been  found.  The  restriction, 


,  imposed  to  insure  existence  of  the  Mellin 


transforms,  is  interesting  in  view  of  the  Laplace- transform  approach 
to  the  transient  problem.  As  mentioned  above,  replacement  of  (W) 
by  a  real  p  ,  and  multiplication  of  the  result  by  ¥<r> 
should  give  the  Laplace  transform  of  the  transient  solution.  Such 
an  identification  is  excluded  by  the  restriction 

This  inability  to  extend  the  solutions  to  real  p  ,  hence  failure 
to  get  a  Laplace- transform  solution  of  the  transient  problem,  is 
consistent  with  the  failure  in  Section  7*2  to  get  a  similarity 
solution  in  ,  for  small  time.  However,  these 

circumstances  lead  to  the  remarkable  conclusion  that,  contrary 
to  the  usual  experience  in  such  problems  (in  particular  the 
classical  unsteady  solution  discussed  in  Sections  5-2  and  5.3),  the 
transient  response  for  small  time  is  not  related  through  its  Laplace 
transform  to  the  high-frequency  response  to  the  corresponding 
problem  of  steady- state  oscillations.  Finally,  it  should  be 
remarked  that  Spence's  latest  approach  to  the  transient  problem, 
if  correct,  must  be  made  to  tie  in  consistently  with  the  above 
solution  for  steady  oscillations.  It  must  be  concluded,  then, 
that  the  solutions  given  in  this  section  for  steady  oscillations 


must  be  regarded  as  tentative,  subject  to  further  study,  both  of 
them  and  Spence's  new  results. 


CHAPTER  8  -  CGHGLUBXCKS 


By  extension  of  the  existing  steady,  Jet- flapped.,  thin- 
airfoil  and  classical,  unsteady,  thin- airfoil  theories,  a  model 
for  unsteady  motions  of  Jet- flapped  thin  airfoils  has  been 
formulated  in  the  first  four  chapters.  Study  of  these  existing 
theories  has  clarified  certain  features  of  the  unsteady  Jet-  flapped- 
airfoil  behavior.  Such  studies  have  suggested  the  feasibility  of  a 
"boundary- layer"  transformation,  which  amplifies  the  region  near 
the  trailing  edge  of  the  airfoil  for  small  Jet-momentum  strengths, 

,  and  for  either  small  times  after  initiation  of 
transient  motions  or  for  high-frequency  steady- state  oscillations. 

The  invalidity  of  the  small-time  "solution"  found  by 
Spence  in  III  for  the  Jet-deflection  problem  has  been  pointed 
out  and  discussed  in  detail.  No  correct  solutions  using  Spcncc ' a 
approach  could  be  found,  however.  Likewise  Spence's  error  in  III 
for  the  high-frequency  response  to  steady- state  oscillations 
in  Jet  deflection  has  been  pointed  out.  For  this  problem  a 
tentative  solution,  as  yet  not  fully  understood,  has  been  put 
forth.  It  satisfies  the  equations  and  all  conditions  of  the 
problem  in  an  apparently  consistent  sense,  and  can  be  extended 
to  the  corresponding  airfoil-motion  problem. 

Spence's  long-time  solution  in  III,  i.e.,  the  approach 
to  the  steady  solution  after  a  transient  Jet  deflection  or  airfoil 
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motion,  is  the  only  established  result  remaining  for  this  problem. 
Although  not  discussed  in  this  report,  this  solution  is  found  from 
the  full  equations  by  examining  the  neighborhood  %  *  C+Utb 
and  showing  that  in  the  vicinity  of  this  point  for  very  large  times 
there  is  concentrated  an  amount  of  circulation  equal  in  magnitude 
but  opposite  in  sign  to  the  total  circulation  of  airfoil  plus  jet 
in  the  steady  solutions  of  I  and  II.  Considering  the  interaction 
of  this  "starting  circulation"  and  the  "steady- flow  circulation," 
the  lift  coefficient  for  approach  to  the  steady  solution  was  found 
to  be 


- c*]  (£i)  W&) 


(8-1) 

where  is  the  lift-curve  slope  in  the  steady, 

Incidence,  Jet-flap  solution,  say  equation  (65)  of  II,  and  -fc 
is  the  physical  time. 

The  surface  has  only  been  scratched  in  finding  the  lift 
and  pitching-moment  responses  to  unsteady  motions  of  jet- flapped 
airfoils.  The  approach  through  the  "boundary- layer"  transformation, 
if  solved  would  only  give  results  in  a  limited  range  of  and  b 

or  U)  .  Nevertheless  these  equations  have  the  strength  of 
their  relative  simplicity  and  give  some  hope  for  further  analytic 
attempts  to  solve  them.  It  remains,  as  mentioned  in  the  previous 
chapter,  to  investigate  fully  Spence's  new,  as  yet  unpublished, 
solution  for  the  transient  case,  as  well  as  the  tentative  high- 
frequency  solutions  for  steady  oscillations  given  in  that  chapter. 
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With  the  resolution  and  understanding  of  these  limiting 
results,  digital  computation  of  the  full  equations  would  probably 
be  necessary  to  give  solutions  for  all  intermediate  times  between 
small  and  large,  and  for  all  frequencies  up  to  the  high  ones.  As 
briefly  discussed  in  Section  5.4,  Spence  in  I  obtained  rapid 
convergence  to  the  steady- state  solution  using  a  numerical  collocation 
scheme,  since  the  functions  being  sought  had  a  monotonic  behavior 
in  %  .  A  collocation  scheme  should  also  be  applicable  in  the 

unsteady  problems,  although  convergence  would  by  no  means  be  as  rapid, 
due  to  the  much  more  ccmplicated  behavior,  in  <90-  and  b  or  W 
of  the  unsteady  solution.  Furthermore,  the  collocation  points  would 
have  to  be  chosen  in  a  fashion  to  adequately  handle  the  important 
effects  in  the  immediate  vicinity  of  the  trailing  edge,  'TL'*  C  , 
and  those  in  the  vicinity  of  »  d-'tUo'fc  or  , 

say. 

Finally,  and  probably  most  important,  is  the  need  for 
definitive  experiments  to  test  the  validity  of  the  model  formulated 
here,  and  any  solutions  which  might  be  found.  Only  in’ this  way  can 
the  ultimate  value  of  this  theory  be  established. 

It  is  strongly  felt  that  further  research  along  the  above 
lines  is  worthwhile,  not  only  to  obtain  results  of  use  to  the 
design  engineer,  but  also  to  understand  the  interesting  mathematical 
and  physical  points  raised  by  the  present  model  and  its  equations. 
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APFEHZZIX  A 


Evaluation  of  Certain  Integrals 

Although  many  of  the  integrals  to  be  evaluated  in  this 
Appendix  may  be  found  elsewhere,  it  is  convenient  to  treat  then 
in  general  and  collect  their  results  here  for  the  particular 
applications  required  in  the  text. 

The  first  type  of  integral  to  be  treated  is  defined  by 


(A-1) 


for  £C  ,  with  integer  values  of  Yi  , 

such  that  0  ^ -I  .  The  superscript,  O  ,  refers  to 

this  type  of  Integral  and  the  subscript,  H  ,  to  the  exponent 
of  %  .  The  integral  is  readily  evaluated  by  making  the 

transformation. 


f  *  c  Sin1 6 

dix 


(A-2) 


which  gives 


(A-3) 
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Using  the  trigonometric  Identity  given  by  extension  at  equation  4o4 
of  Dwight  09^7);  namely 

*»« 

*  (T«)  + 


Ml 


(A-4) 


(A- 3)  nay  be  inmediately  integrated  to  give 


(«.f) = -  «4#]} 

4 


(A- 5) 


A  particular  case  of  this  is 


T>.c)  *  (Cf-7). 


(A-6) 


The  examples  of  (A- 5)  and  (A-6)  which  are  required  in  the  text  are 

x:  *  t  e^oser]  -  ± 


(*-7) 


X*f«.d  *  £  , 


(A-8) 


(Ar9) 
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I 


(A- 10) 


and 


I>.0  , 


(A-tl) 


The  second  type  of  integral  is 

x;m,p).x| Lrfyj*. 

V  (A- 12) 

where  0&c(*^£C  with  integer  values  of  V)  ,  such  that 
Y\  •  Although  this  is  readily  evaluable 

using  the  transformation  (A-2) ,  it  nay  be  related  directly  to  the 
previous  integral,  I*  l&,$)  .  With  the  identity 


(A- 13) 


it  is  seen  from  (A-1)  that 


xk/o*  -  iu*,?). 


(A-14) 


A  particular  cause  is 


(A- 15) 


The  required  examples  are 

/iuffo-u,tA 

(A- 16) 

I,'  Uuf^'C)  a  c.  twvwO, 

(A- 17) 

I'fo.e)  •  £  . 

(A- 18) 

I>, 

<1  i  (Ucif’/M-U.i)*- 

"  /Mr  X1,  *"  A*  > 

(A- 19) 

c)  -- 


un 


2*n 


un 


(A-20) 


and 


a  i*,c)  * 
t;  (c<4*±td  * 


(A-21) 

^+^1-  D*'  904  1 tUf  iC»f\l 

MTT  3M7T  x  J> 

(A- 22) 


r,'  fo,  o 


/»? 


(Ar23) 
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The  third  type  of  integral  is 
/* 


/tf  Ml 
(~sj  «->  > 


(A- 24) 


where  o 6 ±C  and  h  is  an  integer  such 

that  ft  ^  /  .  The  integral  exists  for  O* 

^  4.%  +  c  and  ^7C  ,  and  in  the 

sense  of  the  Cauchy  Principal .  Value  for  (X  *■  OC  p 
Substitution  of  (A-2)  makes  it 

•  w-fftr-] 


Sot  &_a0 
Sk*©-  £ 


The  integrand  may  be  rewritten,  using  the  Identity 


(A^25) 


•  i/ft*) 


Jtf  if  Iff**'  '  , 


(A- 26) 


and  the  definition  of  X*/«,p)  ,  (A-1) 


I.  Kp.*>  -  % 


as 


^Riti 

r 


lu® 


16 


?-*>*»  . 


(A-27) 


The  remaining  integral  may  be  evaluated  using  equation  436-7  of 
Dwight  (1947),  as 
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f_i£_  = 

J  i-*» 


c*x**». 

(a- 29) 


The  final  result  is,  in  general. 


9 


o±%<  c 

(A-30) 


c  <■%<■*>. 


(A-31 ) 


A  particular  case  is 
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(A- 32) 


(A-33) 


The  particular  integrals  necessary  to  the  text  cure 


-m-wn 


FgPf 


►Wf 


(A-34) 


(A-35) 


o*  x*c 


(A- 36) 


c*-%*.oc 


(A- 37) 
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where  O  S  «<►  ^  SC  and  is  an  integer  such  that 

M  O  .As  for  i:  iAf.%)  ,  it  exists 

for  o*-%*’*  ,  and  , 

and  in  the  sense  of  the  Chubby  Principal  value  for  (X  *  %*■  $ 

Use  of  the  identity,  (A-13'/>  reduces  it  to 

I*  -  c.  I*.,  Id f,%)  -  . 

(*-1*5) 

A  particular  case  of  interest  is 


r 


o  *%.<-<- 


(A-W) 


a*» 

(A-47) 


The  integrals  required  in  the  text  are 


(ArW). 


-4  s* 


(Ar^9) 
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1 1  t‘,c,x) 


-I  ♦  (Jtf, 


(Ap50) 


(a- 51) 


-%  +  4 


ce%*«o  , 


and 


(A- 52) 


(A- 53) 


X l(0,C& 


r  «A%*C 


(A-5^) 


(Ar55) 


